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We present a new general algorithm for calculating arbitrary jet cross sections 
in arbitrary scattering processes to next-to-leading accuracy in perturbative 
QCD. The algorithm is based on the subtraction method. The key ingredi- 
ents are new factorization formulae, called dipole formulae, which implement 
in a Lorentz covariant way both the usual soft and collinear approximations, 
smoothly interpolating the two. The corresponding dipole phase space obeys 
exact factorization, so that the dipole contributions to the cross section can be 
exactly integrated analytically over the whole of phase space. We obtain explicit 
analytic results for any jet observable in any scattering or fragmentation process 
in lepton, lepton-hadron or hadron-hadron collisions. All the analytical formu- 
lae necessary to construct a numerical program for next-to- leading order QCD 
calculations are provided. The algorithm is straightforwardly implementable in 
general purpose Monte Carlo programs. 
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1 Introduction 



Most of the recent progress in the understanding of strong interaction physics at large 
momentum transfer has been due to the comparison between precise experimental data 
and very accurate QCD calculations to higher perturbative orders 0. 

The perturbative QCD approach for computing hadronic cross sections is based on the 
parton model picture. According to this picture, the cross section for any hard-scattering 
process (i.e. any process involving a large transferred momentum Q) can be written as 
a convolution of structure (f a (%,Q 2 )) and fragmentation (d a (x,Q 2 )) functions of partons 
(quarks and gluons) and a hard-cross section factor. The former are non-perturbative 
quantities but are universal, that is, they are process independent. The latter is instead 
dominated by momentum regions of the order of Q and hence, provided that Q ^> A (where 
A is the QCD scale), it can be computed in QCD perturbation theory to the lowest order 
in the 'small' (due to asymptotic freedom) running coupling «s(Q) ~ (A) In Q 2 /A 2 )~ 1 . 

This naive parton model approach corresponds to the so-called leading-order (LO) ap- 
proximation. It is justified by the high-momentum behaviour of the running coupling 
«s(Q). However, just because of its perturbative nature, the running of the QCD cou- 
pling can be hidden in higher-order corrections by the replacement as{Q) = aig [1 + 
K(Q) as(Q) + ...], being the values of ag at a fixed (and arbitrary) momentum scale. 
It follows that a LO calculation predicts only the order of magnitude of a given cross section 
and the rough features of a certain observable^ The accuracy of the perturbative QCD 
expansion is instead controlled by the size of the higher-order contributions. Any defi- 
nite perturbative QCD prediction thus requires (at least) a next-to-leading order (NLO) 
calculation, and NLO definitions of as, f a (x,Q 2 ) and d a (x,Q 2 ). 

This is the reason why the results of the higher-order QCD calculations which are 
available at present have been proved to be of vital importance to assess the progress 
mentioned above. 

These higher-order computations have been carried out over a period of about fifteen 
years, often long after the accuracy of experimental data has made them necessary, be- 
cause of the difficulties in setting up a general and straightforward calculational procedure. 
The physical origin of these difficulties is in the necessity of factorizing the long- and 
short-distance components of the scattering processes and is reflected in the perturbative 
calculation by the presence of divergences. 

In general, when evaluating higher-order QCD cross sections, one has to consider real- 
emission contributions and virtual-loop corrections and one has to deal with different kind of 
singularities. The customary ultraviolet singularities, present in the virtual contributions, 
are removed by renormalization. The low-momentum (soft) and small-angle (collinear) 
regions instead produce singularities both in the real and in the virtual contributions. In 
order to handle these divergences, the observable one is interested in has to be properly 
defined. It has to be a jet quantity, that is, a hadronic observable that turns out to be 

^ Typically, this poor predictivity is quantitatively signalled by a strong dependence on the (unphysical) 
renormalization and factorization scales. 
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infrared safe and either collinear safe or collinear factorizable: its actual value has to be 
independent of the number of soft and collinear particles in the final state (see Sect. [7| 
for a formal definition). In the case of jet quantities, the coherent sum over different 
(real and virtual) soft and collinear partonic configurations in the final state leads to the 
cancellation of soft singularities. The left-over collinear singularities are then factorized 
into the process-independent structure and fragmentation functions of partons (parton 
distributions), leading to predictable scaling violations. As a result, jet cross sections 
are finite (calculable) at the partonic level order by order in perturbation theory. All the 
dependence on long-distance physics is either included in the parton distributions or in non- 
perturbative corrections that are suppressed by inverse powers of the (large) transferred 
momentum Q that controls the scattering process. 

Because of this complicated pattern of singularities, the simplest quantities that can be 
computed in QCD perturbation theory are fully inclusive. In this case one considers all 
possible final states and integrates the QCD matrix elements over the whole available final- 
state phase space. Thus one can add real and virtual contributions before performing the 
relevant momentum integrations in such a way that only ultraviolet singularities appear at 
the intermediate steps of the calculation. Owing to this simplification, powerful techniques 
have been set up |§ to perform analytic calculations up to next-to-next-to-leading order 
(NNLO), i.e. to relative accuracy with respect to the lowest-order approximation. 

In the case of less inclusive cross sections, QCD calculations beyond leading order 
(LO) are much more involved. Owing to the complicated phase space for multi-parton 
configurations, analytic calculations are in practice impossible for all but the simplest 
quantities. However, the use of numerical methods is far from trivial because real and 
virtual contributions have a different number of final-state partons and thus have to be 
integrated separately over different phase space regions. Unlike the case of fully inclusive 
observables, one cannot take advantage of the cancellation of soft and collinear divergences 
at the integrand level. Soft and collinear singularities, present in the intermediate steps, 
have to be first regularized, generally by analytic continuation in a number of space-time 
dimensions d = 4 — 2e different from four. Then the real and virtual contributions should be 
calculated independently, yielding equal-and-opposite poles in e. Great progress has been 
made in recent years in the analytical techniques for calculating the virtual processes , but 
the analytic continuation greatly complicates the Lorentz algebra in the evaluation of the 
matrix elements and prevents a straightforward implementation of numerical integration 
techniques. Despite these difficulties, efficient computational techniques have been set up, 
at least to NLO, during the last few years. 

There are, broadly speaking, two types of algorithm used for NLO calculations: one 
based on the phase-space slicing method and the other based on the subtraction method^. 
The main difference between these algorithms and the standard procedures of analytic cal- 
culations is that only a minimal part of the full calculation is treated analytically, namely 
only those contributions giving rise to the singularities. Moreover, for any given pro- 
cess, these contributions are computed in a manner that is independent of the particular 
jet observable considered. Once every singular term has been isolated and the cancella- 
tion/factorization of divergences achieved, one can perform the remaining part of the cal- 

*We refer the reader to the Introduction of Ref. jl| for an elementary description of the basic difference 
between the two methods. 
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culation in four space-time dimensions. Although, when possible, one still has the freedom 
of completing the calculation analytically, at this point the use of numerical integration 
techniques (typically, Monte Carlo methods) is certainly more convenient. First of all, 
the numerical approach allows one to calculate any number and any type of observable 
simultaneously by simply histogramming the appropriate quantities, rather than having to 
make a separate analytic calculation for each observable. Furthermore, using the numerical 
approach, it is easy to implement different experimental conditions, for example, detector 
acceptances and experimental cuts. In other words, the phase-space slicing and subtrac- 
tion algorithms provide the basis for setting up a general-purpose Monte Carlo program 
for carrying out arbitrary NLO QCD calculations in a given process. 

Both the slicing |5j and the subtraction || methods were first used in the context of 
NLO calculations of three-jet cross sections in e + e~ annihilation. Then they have been 
applied to other cross sections adapting the method each time to the particular process. 
Only recently has it become clear that both algorithms are generalizable in a process- 
independent manner. The key observation is that the singular parts of the QCD matrix 
elements for real emission can be singled out in a general way by using the factorization 
properties of soft and collinear radiation [[7]. 

At present, a general version of the slicing algorithm is available for calculating NLO 
cross sections for production of any number of jets both in lepton || and hadron 
collisions. To our knowledge, fragmentation processes have been considered only in the 
particular cases of prompt-photon production QTUJ and single- and double-hadron inclusive 
distributions |11,12||. The complete generalization of this method to include fragmentation 



functions and heavy flavours is in progress [13]. 



As for the subtraction algorithm, a general NLO formalism has been set up for comput- 
ing three-jet observables in e + e~ annihilation [§,0 and cross sections up to two final-state 
jets in hadron collisions^] |4],|]1|]. Also the treatment of massive partons has been considered 
in the particular case of heavy-quark correlations in hadron collisions p~ 



In this paper, we present a completely general version of the subtraction algorithm. This 
generality is obtained by fully exploiting the factorization properties of soft and collinear 
emission and, thus, deriving new improved factorization formulae, called dipole factorization 
formulae. They allow us to introduce a set of universal counterterms that can be used for 
any NLO QCD calculation. Therefore, our version of the subtraction method can be 
compared with those used for three-jet observables in e + e~ annihilation |J and two-jet 
quantities in hadron collisions [[| (although, in these known cases our treatment turns out 
to differ in many respects from the previous ones). Moreover, we are able to consider 
the production of any number of jets in lepton and hadron cross sections and to provide a 
general treatment of fragmentation processes and multi-particle correlations. The inclusion 
of heavy quarks in the algorithm can also be performed in a completely general and process- 



independent manner | 19|[ . The extension of our method to polarized scattering is not 



considered here but it is straightforward. 

§ After completion of the present work, the method of Ref. Q] has been modified to deal with three-jet 
cross sections |l6j| . The formalism presented in Ref. jl6| can be extended to n-jet production both in lepton 
and hadron collisions. A similar method has been presented in Ref. fVE. 
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Besides discussing in detail our general formalism, we explicitly carry out the <i-dimen- 
sional analytical part of the NLO calculation for all the (unpolarized) scattering processes 
involving massless quarks and gluons. Knowing the relevant QCD matrix elements, the re- 
sults of our algorithm can be straightforwardly implemented in NLO numerical codes with- 



out any additional calculation. Detailed numerical applications to e + e annihilation pO 



and deep inelastic lepton-hadron scattering |21j are presented elsewhere. 



We begin in Sect. [| by giving a brief overview of the general method, describing the 
subtraction procedure and how our dipole formulae are used to implement it. In Sect. |3| we 
establish the notation used throughout the paper. In Sect. |4] we review the factorization 
properties of QCD matrix elements in the soft and collinear limits before presenting, in 
Sect. |5|, our dipole factorization formulae, which smoothly interpolate these two limiting 
regions. After briefly recalling, in Sect. ||, the precise definitions of QCD cross sections 
at NLO, we go on to describe in detail our subtraction method for evaluating these cross 
sections, in Sects. 0-11. In Sect. 12 we summarize and discuss our results. Appendix A 



gives more details, and some examples, of the necessary colour algebra. In Appendix B we 
explicitly perform the only difficult integral we encounter. In Appendix C we collect to- 
gether the main formulae needed to implement our method in specific calculations. Finally 
in Appendix D we work through a few simple examples of applying our method to specific 
cross sections. 

Since the paper is quite long, readers mainly interested in understanding the general 
method or in some particular application are advised to first read Sects. 0, |, [|, [5J] and [7[ 



Here we discuss in detail our general formalism and its use for processes with no initial-state 
hadrons like, for instance, e + e~ annihilation (in this case, a brief description of our method 
has already appeared p0|). Sections |8Hlll and the other Subsections in Sect. ^| can then be 



read quite independently from one another. The final formulae that are necessary for the 
actual numerical implementation of our algorithm in each different scattering process are 
summarized in a Subsection at the end of each of Sects. 0-4TI • 
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2 The general method 



In this Section we explain the general idea behind our version of the subtraction method 
by describing the subtraction procedure (Sect. p.l|) and considering mainly the simplified 
case of jet cross sections in processes with no initial-state hadrons, for instance, e + e~ 
annihilation (Sect. |2.2j ). A brief description of our method for more complicated scattering 
processes is sketched in Sect. p73[ 



2.1 The subtraction procedure 



Suppose we want to compute a jet cross section a to NLO, namely 

_LO | „nlo 
a = a + a 



(2.1) 



Here the LO cross section a LO is obtained by integrating the fully exclusive cross section 
da B in the Born approximation over the phase space for the corresponding jet quantity. 
Suppose also that this LO calculation involves m partons in the final state. Thus, we write 



a LO = da B , (2.2) 

Jm 



where, in general, all the quantities (QCD matrix elements and phase space) are evaluated 
in d = 4 — 2e space-time dimensions. However, by definition, at this LO the phase space 
integration in Eq. ( |2.2| ) is finite so that the whole calculation can be carried out (analytically 
or numerically) in four dimensions. 

Now we go to NLO. We have to consider the exclusive cross section da R with m + 1 
partons in the final-state and the one-loop correction da v to the process with m partons 
in the final state: 

a NLO _ I da NLO _ I da R + I da V (2 3) 

J J m+1 Jm 

The two integrals on the right-hand side of Eq. ( |2.3[ ) are separately divergent if d = 4, 
although their sum is finite. Therefore, before any numerical calculation can be attempted, 
the separate pieces have to be regularized. Using dimensional regularization, the diver- 
gences (arising out of the integration) are replaced by double (soft and collinear) poles 1/e 2 
and single (soft, collinear or ultraviolet) poles 1/e. Suppose that one has already carried out 
the renormalization procedure in da v so that all its ultraviolet poles have been removed. 

The general idea of the subtraction method for writing a general-purpose Monte Carlo 
program is to use the identity 



da 



NLO 



da R - da' 



+ da A + da v , (2.4) 



where da A is a proper approximation of da R such as to have the same pointwise singular 
behaviour (in d dimensions) as da R itself. Thus, da A acts as a local counterterm for da R 
and, introducing the phase space integration, 



a NLO 



m+1 



da R - da J 



+ / da A + / da v , (2.5) 

J m+1 Jm 
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one can safely perform the limit e — > under the integral sign in the first term on the 
right-hand side of Eq. fl2.5| ). Hence, this first term can be integrated numerically in four 
dimensions. 

All the singularities are now associated to the last two terms on the right-hand side 
of Eq. (|2.5| ). If one is able to carry out analytically the integration of do A over the one- 
parton subspace leading to the e poles, one can combine these poles with those in da v , thus 
cancelling all the divergences, performing the limit e — > and carrying out numerically the 
remaining integration over the m-parton phase space. The final structure of the calculation 
is as follows 

(2.6) 

e=0 

and can be easily implemented in a 'partonic Monte Carlo' program, which generates 
appropriately weighted partonic events with m + 1 final-state partons and events with m 
partons. 

Note that the subtracted term [ da R — da A ] in Eq. ( j2.6|) is integrable in four dimensions 
by definition. The fact that all the divergences cancel in the second term on the right- 
hand side of Eq. (|2.6j ) is instead not a general feature of all hadronic cross section^]. The 
cancellation of divergences is guaranteed only for the hadronic observables that we are 
considering in this paper, namely jet observables. 

These quantities have to be experimentally (theoretically) defined in such a way that 
their actual value is independent of the number of soft and collinear hadrons (partons) 
produced in the final state. In particular, this value has to be the same in a given m-parton 
configuration and in all m + 1-parton configurations that are kinematically degenerate with 
it (i.e. that are obtained from the m-parton configuration by adding a soft parton or 
replacing a parton with a pair of collinear partons carrying the same total momentum). 
This property can be simply restated in a formal way. If the function Fj 1 ^ gives the value 
of a certain jet observable in terms of the momenta of the n final-state partons, we should 
have 

F (m+l) _^ F (m) ^ (2 7) 

in any case where the m + 1-parton and the m-parton configurations are kinematically 
degenerate. 

The Born-level cross section da B can be (symbolically) written as a function of the 
jet-defining function Fj in the following way 

da B = d& m) \M m \ 2 Fj m) , (2.8) 

where <i$( m ) and Ai m respectively are the full phase space and the QCD matrix element 
to produce m final-state partons. The corresponding expression for the real cross section 
da R is: 

da R = d& m+ ^ \M m+1 \ 2 Fj m+1) . (2.9) 

^The presence of singularities in a QCD cross section computed in perturbation theory does not mean 
that the theory itself is inconsistent. It simply means that one is considering a cross section that cannot be 
reliably estimated using the perturbative expansion. At any energy scale, it is affected by non-perturbative 
phenomena that are as big as the perturbative ones. 
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a 



NLO 



rn+l 



+ 



da v + / do- 



The structure of Eq. (|2.9|) and the fundamental property (2/7) are essential for the 
feasibility of the subtraction procedure described in this Subsection. There are obviously 
many ways of approximating the matrix element Ai m +i i n the neighbourhood of its soft 
and collinear singularities. Correspondingly, one can approximate Fj m+1 ^ and obtain a local 
counter-term da A . The main point is that, due to the limiting behaviour in Eq. fl2.7|) , one 
can always find an approximation for Fj m+1 ^ such that the one-parton subspace leading to 
the soft and collinear divergences effectively decouples. Thus, one can perform the integral 
j l da A and the subtraction formula ( [2.6|) can, in principle, always be implemented. 



2.2 Dipole formulae and universal implementation of the sub- 
traction procedure for jet cross sections 

The key for the subtraction procedure to work in practice is obviously the actual form of 
da A . One needs to find an expression for da A that fulfils the following properties: i) for 
any given process, da A has to be obtained in a way that is independent of the particular 
jet observable considered; ii) it has to exactly match the singular behaviour of da R in 
d dimensions; Hi) its form has to be particularly convenient for Monte Carlo integration 
techniques; iv) it has to be exactly integrable analytically in d dimensions over the single- 
parton subspaces leading to soft and collinear divergences. 

In Ref. 0, a suitable expression for da A for the process e + e~ — > 3 jets was obtained by 
starting from the explicit expression (in d dimensions) of the corresponding da R and by 
performing extensive partial fractioning of the 3 + 1-parton matrix elements, so that each 
divergent piece could be extracted. This is an extremely laborious and ungeneralizable 
task, in the sense that having done it for e + e~— > 3 jets does not help us to do this for, say, 
e + e~—> 4 jets or for any other process. 

In Ref. [|j], the general properties of soft and collinear emission were first used (in the 
context of the subtraction method) to construct do~ A , for one- and two-jet production in 
hadron collisions, in a way that is independent of the detailed form of the corresponding 
da R . 

The central proposal of our version of the subtraction method is that one can give a 
recipe for constructing da A that is completely process independent (and not simply indepen- 
dent of the jet observable). Starting from our physical knowledge of how the m + 1-parton 
matrix elements behave in the soft and collinear limits that produce the divergences (see 
Sect. |), we derive improved factorization formulae, called dipole formulae (see Sect. ^), 
that allow us to write: 

da A = ]T da B <g> dV dipole . (2.10) 

dipoles 



The notation in Eq. ( |2.10 ) is symbolic. Here da B denotes an appropriate colour and spin 



projection of the Born-level exclusive cross section. The symbol (g> stands for properly 
defined phase space convolutions and sums over colour and spin indices. The dipole factors 
dVdipoie (which match the singular behaviour of da R ) are instead universal, i.e. completely 
independent of the details of the process and they can be computed once for all. In 
particular, the dependence on the jet observable is completely embodied by the factor da B 
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of Eq. (fnCD , in the form of Eq. Q. 



There are several dipole terms on the right-hand side of Eq. ( p. 10 ). Each of them 



corresponds to a different kinematic configuration of m + 1 partons. Each configuration 
can be thought as obtained by an effective two-step process: an m-parton configuration is 
first produced and then one of these partons decays into two partons. It is this two-step 
pseudo-process that leads to the factorized structure on the right-hand side of Eq. ( |2.10| ). 

The reason for having several dipoles is that each of them mimics one of the m+l-parton 
configurations in da R that are kinematically degenerate with a given m-parton state. Any 
time the m + 1-parton state in do R approaches a soft and/or collinear region, there is a 
corresponding dipole factor in da A that approaches the same region with exactly the same 
probability as in da R . In this manner do A acts as a local counter-term for da R . 

Our expression for da A in Eq. (|2.10| ) is completely defined over the full m + 1-parton 
phase space (in particular, da A does not depend on any additional phase space cut-off]): 
momentum conservation is exactly implemented in each term on the right-hand side of 
Eq. ( |2~T0l ) and there is a one-to-one correspondence between each partonic configuration 
in da R and (each of the several) in da A . Therefore, in our case, da R — da A is straight- 
forwardly integrable via Monte Carlo methods: one generates an m + 1-parton event with 
weight da R and, correspondingly, one can obtain an m+l-parton counter-event with weight 
da A . 

Furthermore, the product structure in Eq. fl2.10|) allows us a factorizable mapping from 
the m + 1-parton phase space to an m-parton subspace (that identified by the partonic 
variable in da B ) times a single-parton phase space (that identified by the dipole partonic 
variables in G?Vdi po ie)- This mapping makes G?Vdi po ie fully integrable analytically. We can 
write (again, symbolically): 

da A = J2 I da B ®j dV dipole = f \da B ® i] , (2.11) 

m+1 dipoles Jm Jl Jm 

where the universal factor I is defined by 

1= E / dipole , (2.12) 
dipoles 



and contains all the e poles that are necessary to cancel the (equal and with opposite sign) 
poles in da v . 

The structure of the final result is given as follows in terms of two contributions 
a NLO{m+\}^ (J NLO{m} m _j_ i_p ar ton and m-parton kinematics, respectively) which 

"This is quite a non-trivial feature of our approach. The most naive way of setting up subtraction 
procedures based on universal properties in the soft and collinear limits would lead to the introduction of 
energy and angular cut-offs, thus breaking Lorentz covariance at intermediate steps. Alternative and less 
naive methods for imposing soft and collinear cut-offs have their own disadvantages, too. For instance, 
the cut-off can be related to some kinematic invariant of the process (see Ref. Q), or one can introduce a 
Lorentz covariant cut-off on parton-parton invariant masses (see Ref. (8|) rather than energies and angles. 
In the first case one could spoil the universality of the subtraction procedure making it process dependent. 
In the second case, it is quite difficult to arrange the cut-off in such a way that the subtraction term is 
exactly integrable analytically to any accuracy in the cut-off itself. 
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are separately finite (and integrable) in four space-time dimensions: 

NLO 



a' 



a NLO{m+l} _j_ a NLO{m} 



(2-13) 



m+l 



da 1 



e=0 



- £ da 



B 



dVr 



dipole 



idipoles 



e=0 



+ 



da v + da 1 



t=0 



Equation (|2.13| ) represents our practical implementation of the general subtraction formula 



In this paper we provide explicit expressions for both the universal factors <iVdi po ie an d I- 
Having these factors at our disposal, the only other ingredients necessary for the full NLO 
calculation, according to Eq. (|2.13|) , are the following (reading Eq. (|2.13|) from the right to 
the left): 



• a set of independent colour projections^ of the matrix element squared at the Born 
level, summed over parton polarizations, in d dimensions; 

• the one-loop contribution da v in d dimensions; 

• an additional projection of the Born level matrix element over the helicity of each 
external gluon in four dimensions; 

• the real emission contribution da R in four dimensions. 



These few ingredients are sufficient for writing, in a straightforward way, a general-purpose 
NLO Monte Carlo algorithm. Note in particular that there is no need to extract a proper 
counter-term da A starting from a cumbersome expression for da R in d dimensions. The 
NLO matrix element contributing to da R can be evaluated directly in four space-time 
dimensions thus leading to an extreme simplification of the Lorentz algebra, particularly 
if one makes use of helicity amplitudes |22| to control the rapid increase in the number of 
Feynman diagrams as the number of parton grows. 



2.3 Factorization of collinear singularities and general algorithm 
for processes with identified hadrons 

The discussion in Sects. p.l| and [2.2| applies to all the processes with no initial-state hadrons 
(for instance, e + e~ annihilation). However, perturbative QCD can be used also for the cal- 
culation of jet cross sections in lepton-hadron and hadron-hadron collisions^]. The main 
difference is that the presence of initial-state hadrons (partons), carrying a well defined 
momentum, spoils the cancellation of the collinear singularities arising in the perturba- 
tive treatment. The left-over singularities can be factorized and reabsorbed into non- 
perturbative and universal (process-independent) distribution functions, the parton densi- 
ties of the incoming hadron. Provided this factorization procedure is consistently carried 

**Actually, if the total number of QCD partons involved in the LO matrix element is less than or equal 
to three, one simply needs its incoherent sum over the colours (see Appendix A). 

^Throughout this paper, whenever referring to initial-state hadrons, we implicitly also include hadron- 
like particles such as photons. 
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out, one can thus define parton-initiated jet cross sections (see Sect. |]) that are free from 
singularities and can be computed with a subtraction procedure similar to that described 



in Sect. 2.1 



Similar features appear when the jet observable depends on the actual value of the 
momentum of one or more hadrons observed in the final state (the inclusive one-particle 
distribution in e + e~ annihilation is the simplest example). Also in this case there are 
left-over collinear singularities that can be reabsorbed into non-perturbative and universal 
distribution functions, the fragmentation functions of the outgoing hadron. As a result, 
one can again define partonic cross sections that are free from singularities and computable 
in perturbation theory (see Sect. 

Because of these common features, in this paper the processes with initial-state hadrons 
and those involving fragmentation functions will be referred to as processes with identi- 
fied hadrons (partons). The hadronic cross section is obtained by convoluting partonic 
cross sections with non-perturbative distribution functions. As stated above, the NLO 
partonic cross sections can be evaluated using the subtraction method and, thus, trying 
to implement the subtraction formula in Eq. ( |2.6| ). There are nonetheless some additional 
complications with respect to the case with no identified particles. These complications 
regard the construction of the approximated cross section da A . 

Just as when there are no identified hadrons, the real cross section da R is singular 
whenever a pair of the m + 1 final-state partons become collinear. However in addition, it 
is also singular in the region where one of them becomes collinear to an identified parton. 
Moreover, the phase space integration has to be performed in the presence of additional 
kinematic constraints, related to the fact that the momenta of the identified partons have 
to be kept fixed (or, at most, rescaled by an overall momentum fraction). As for the 
approximated cross section da A , it follows that, on one side, it should act as a local counter- 
term also in the new singular regions and, on the other side, its integral j l da A should still 
be computable analytically even in the presence of the additional phase space constraints. 

The dipole formalism presented in this paper provides a simple and general solution to 
these problems. Indeed, we are able to write the cross section da A in the following form 
(see Sects. |-0) 

da A = ]T da B <g> (dV dipolc + d^ ip oie) . (2.14) 

dipolcs 



Equation fl2.14|) is completely analogous to Eq. (|2.10|) . The additional dipole terms dV dipol 



on the right-hand side match the singularities of da R coming from the region collinear to 
the momenta of the identified partons. Moreover, these dipole terms are still (i.e. even 
if the momenta of the identified partons are fixed) fully integrable analytically over the 
one-parton subspace leading to soft and collinear divergences. 

These are peculiar features of the dipole approach. As discussed below Eq. ( [2.10|) , each 
dipole contribution is effectively obtained by first producing an m-parton configuration 
and then letting one parton to decaying into two partons. The dipole formulae implement 
this two-step procedure by enforcing exact momentum conservation. Actually there are 
equivalent ways of doing that, corresponding to different ways of treating the momentum 
recoil in the m-parton configuration. This freedom allows us to define alternative versions 
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of the factorization formulae (Sect. |5|) and, hence, different dipole factors (like <iVdi po ie 
and dV,x pole in Eq. ( |2.14| )). These differences are then used to match (and overcome) the 
phase-space constraints that are encountered in the calculation of QCD cross sections with 
identified particles. 



Having introduced the counter-term da A in Eq. (|2.14j) , we can proceed to its integration 
as in Eq. ( [2.11| ). We thus obtain the singular factor J in Eq. (|2.12|) and additional singular 
terms, which are reabsorbed into the non-perturbative distribution functions. 

The final result of our subtraction procedure is given in terms of the NLO partonic 
cross section a NLO (p), where the dependence on the momentum p symbolically denotes 
the functional dependence on the momenta of the identified partons. This cross section is 
obtained by a formula that is similar to Eq. ( |2.13|) , namely (see Sects. jSl-[Tl~l) 



a 



.NLO 



'(?) = 



a 



NLO { m + l } (p) + a NLO {m} ,j + f 1 dx d NLO { m} , ) 

JO 



m+1 



da R (p) 



e=0 



+ 



X; da B (p)®(dV dipole + dV^ oley 

^dipoles j £ _q 

»1 



(2.15) 



da v (p) + da B (p) <g> I 



dx 

e =0 Jo -'in 



da (xp) ® (P + K + H) (x) 



e=0 



Here, the contributions a NL °^ m+1 ^(p) and a NLO ^ (p) (with m + 1-parton and m-parton 
kinematics, respectively) are completely analogous to those in Eq. Q2.13| ). 



The last term on the right-hand side of Eq. ( p. 15 ) is a finite (in four dimensions) re- 
mainder that is left after factorization of initial-state and final-state collinear singularities 
into the non-perturbative distribution functions (parton densities and fragmentation func- 
tions). This term involves a cross section a NL °^ m ^(x; xp) with m-parton kinematics and an 
additional one-dimensional integration with respect to the longitudinal momentum fraction 
x. This integration arises from the convolution of the Born- type cross section da B (xp) with 
x-dependent functions P, K, H that are similar (but finite for e — > 0) to the factor /. The 
functions P, K and H are universal, that is, they are independent of the detail of the scat- 
tering process and of the jet observable: they simply depend on the number of identified 
partons. Our algorithm provides the explicit expressions for these functions. Therefore, 
writing a general-purpose NLO Monte Carlo program for processes with identified particles 
does not require any further conceptual or analytic effort with respect to the case with no 
identified particles. 
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3 Notation 



3.1 Dimensional regular izat ion 

In general we use dimensional regularization in d = 4 — 2e space-time dimensions and con- 
sider d—2 helicity states for gluons and 2 helicity states for massless quarks (i.e. fermions are 
four-component spinors). This defines the usual dimensional-regularization scheme. Other 
dimensional-regularization prescriptions can be used. However, since the regularization de- 
pendence is unphysical (i.e. it cancels in physical cross sections), within our formalism it is 
more convenient to parametrize it in terms of simple coefficients that enter in the one-loop 
contribution (see Sect. |3~"3"1) . 

The dimensional-regularization scale, which appears in the calculation of the matrix 
elements, is denoted by //. Physical cross sections do not depend on //, although, when 
evaluated in fixed-order perturbation theory, they do depend on the renormalization scale 
Hr and on factorization scales /ip. In other words, the dependence on \i cancels after having 
combined the matrix elements in the calculation of physical cross sections. Therefore, in 
order to avoid a cumbersome notation, we set /i = fiR, while \x and //p will differ in general. 

The d- dimensional phase space, which involves the integration over the momenta {p±, 
p m } of m final-state partons, will be denoted as follows 

{2Tr) d S^(p 1 + ....+ Pm -Q) = d(f> m { Pl ,..., Pm ;Q) . (3.1) 

When there is no ambiguity on the number of final-state partons, we shall drop the subscript 
m in d(j) m . 



m d d pi 

5(2^T* +(P?) 



3.2 Matrix elements 

Let us first consider processes that involve only final-state QCD partons (e + e~-type pro- 
cesses). Non QCD partons (7*, Z°, W^, ■ ■ •), carrying a total incoming momentum Q^, are 
always understood. 

The (tree-level) matrix element with m QCD partons in the final state has the following 
general structure 

m;„ ,>,„) (3.2) 

where {ci, c m }, {si, ...,s m } and {pi,---,p m } are respectively colour indices (a = 1, 
iV| — 1 different colours for each gluon, a = 1, .., N c different colours for each quark or an- 
tiquark), spin indices (/1 = 1, d for gluons, s = 1, 2 for massless fermions) and momenta. 

It is useful to introduce a basis {\ci, c m > <8>|si, s m >} in colour + helicity space 
in such a way that 

M%'---' Cm ' Sl '---' Sm (p 1 ,...,p m ) = (< ci, ...,cj® < s h ...,s m \)\l,....,m> m . (3.3) 
Thus |1, m > m is a vector in colour + helicity space. 
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According to this notation, the matrix element squared (summed over final-state colours 
and spins) |A^ m | 2 can be written as 

\MJ 2 = m < l,....,m\l,....,m > m . (3.4) 



In the following we shall always consider matrix elements squared summed over final- 
state spins (the generalization to fixed-helicity amplitudes is feasible). 

As for the colour structure^], it is convenient to associate a colour charge T, with the 
emission of a gluon from each parton i. If the emitted gluon has colour index c, the 
colour-charge operator is: 

Ti = T t c \c > (3.5) 
and its action onto the colour space is defined by 

< Ob Ci, --Cm, c|Tj|6 1; bi, ..b m >= 5 Cl b 1 ....T^ ibi ...5 Cm b m , (3.6) 

where T® b = ifcab (colour-charge matrix in the adjoint representation) if the emitting par- 
ticle i is a gluon and T£ g = t a aQ (colour-charge matrix in the fundamental representation) if 
the emitting particle i is a quark (in the case of an emitting antiquark T^ B = i% g = —t% a ). 

The colour-charge algebra is: 

Ti ■ Tj = Tj -Ti if i ^ j- T\ = d, (3.7) 

where Cj is the Casimir operator, that is, Cj = Ca — N c if i is a gluon and Ci = Cf = 
(N% — 1)/2N C if i is a quark or antiquark. 

Note that by definition, each vector |l,....,m > m is a colour-singlet state. Therefore 
colour conservation is simply 

m 

Y^Ti |l,....,m> m =0. (3.8) 

4=1 



Using this notation, we also define the square of colour-correlated tree-amplitudes, 
|A4^ fc | 2 , as follows 



\M 



,fe|2 



= m < 1, ...,m| Ti ■ T k i 



1, m >, 



tu t U M';y ■-"> ■■"■■' (m ,>,„) . (3.9) 



In the case of hard processes with QCD partons in the initial state, in addition to m 
final-state partons, the relevant matrix element is: 

^:Z Cm ' Co '-" ;ai '-" ,W "'"-(Pl:---.P»»;Pa,..0 , (3-10) 



^Within our formalism, there is no need to consider the decomposition of the matrix elements into 
colour subamplitudes p3[, as in Ref. M. 
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and the corresponding vector in colour + helicity space will be denoted in the following 

way 



|1, — , 77t, G, ... ^771 a... — 



1 



n c (a) . . 



(|ci, .., c m , c a , .. > .., s m , s a , .. >) 



••'•"(Pi,...,p m ;p a , •••) • 



(3.11) 



Here the labels a, ... refer to the initial-state partons. The normalization of the state vector 



in Eq. (|3.11|) is fixed by including a factor of 1/ ' Jn c {a) for each initial-state parton carrying 
n c (a) colours. 

Note that the colour-charge operator of an initial-state parton a is defined by crossing 
symmetry, that is, {T a ) c hd = if bcd if a is a gluon and {T a ) c a/3 = t c aj3 = -t c pa if a is a quark (if 
a is an antiquark, (T a )^ = t c a p). The analogue of the colour-conservation condition ( |3.8| ) 
is: 



J2Ti + T a + ... |l,....,m;a, ... > m , Q 



. 



(3.12) 



Owing to the normalization of the state vector in Eq. ( 3.11 ), the square of colour- 
correlated tree-amplitudes is: 

\^ma..\ 2 = m,a...< l,...,m;a, ...\T T ■ T j \1, m; a, ... > m ,a... 



nJa 



where the indices J, J refer either to final-state or initial-state partons. 
We refer to Appendix A for more details of the colour algebra. 



rpc rpc 

1 b I a I 1 bja J 



(3.13) 



3.3 One-loop matrix elements and scheme (in)dependence 

We denote by \<M.^^ e }(jpi, ■■■iPmiPa, ■■■)\\i~ioop) ^ ne one -l°op correction to the square of the 
tree-level matrix element in Eq. ( |3.10| ). This term enters in the computation of the vir- 
tual contribution da v to the NLO QCD cross section. Actually, da v is proportional to the 
renormalized one-loop correction \M. m ,a...(pi-> ■■■,Pm)Pa, ■■•)\li-ioop) anc ^ ^he latter is obtained 
from the corresponding bare quantity by simply adding an ultraviolet counterterm. More 
precisely, if the (tree-level) matrix element squared \M. m ,a...(pi, ■■■■>Pm\Pa-> ■■■)\ 2 is propor- 
tional to the n-th power of the QCD coupling as, the renormalized one-loop correction is 
given by 

\M m , a ...(pi,-,Pm;Pa, •••) I = l-M$j£?(pi, Pm] Pa, •••) | fl-loop) ( 3 - 14 ) 

a s (47r) € f Pq ~ \ . . . . . l2 

~ U 27r r(l - e) 1~ + J \M m > a "-v )1 >''->P m >P a >'-)\ ' 

where /3 = (HC^ — 4T R Nf)/6 is the first coefficient of the QCD /^-function (Nf is the 
number of quark flavours and Tr(t a t b ) = 5 ab T R , i.e. Tr = 1/2) and /3q is an e-independent 
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and process-independent coefficient that defines the renormalization scheme. In Eq. (|3.14 ) 
and in the rest of the paper, a s stands for as(/i), the NLO QCD running coupling evaluated 
at the renormalization scale fi. The actual value of the QCD coupling as(fi) depends on 
the renormalization scheme. The customary MS scheme is obtained by setting (3q = in 
Eq. (fug) . 



of \M\f 



The detailed expression of \M^ e }(p u ...,p m ;p a , -)l(i-/ oop ) ( and > hence, ^ \^v\ (1 _ loop)j 
depends on the dimensional regularization procedure used for evaluating the loop inte- 
gral. Since we are using conventional dimensional regularization, we need the result for 
\^rncf > \\i~ioop) within this regularization scheme. If the one-loop correction is known in a 
different scheme, we have to introduce a correction factor proportional to the correspond- 
ing tree-level amplitude. More precisely, we have to perform the following replacement in 
Eq. (fUJl): 



j Pm i Pa > •■■ )\ 



(1—loop) 



\M 

m,a... (Pi,- > Pm j Pa > 
m 

Yl 7i + la + ■ 



as 
2tt 



i=l 



\(l-loop) 
\M m ,a...(Pl, 



(3.15) 



■ > Pm i Pa j 



where the universal coefficients ■ ■ ■ depend only on the flavour of the QCD partons. 

The actual values of these coefficients for several different regularization schemes can be 
found, for instance, in Ref. i24l . 
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4 Factorization in the soft and collinear limits 



4.1 Soft and collinear singularities in tree- level amplitudes 

Let us consider a generic tree-level matrix element M m +x,a... with m + 1 QCD partons 
(Fig. 1). The dependence of \Ai m +i, a ...\ 2 on the momentum pj of a final-state parton 
j is singular in two different phase-space regions: a) in the soft region, defined by the 
limit pj = \q, A — > (where q is an arbitrary four momentum), |-M m +i, a ...| 2 behaves 
as 1/A 2 ; in the collinear region, defined by the limit pj — > (1 — z)pi/z (where pi is the 
momentum of another QCD parton in A4 m +i,a...), \M. m +i,a...\ 2 behaves as l/(piPj). This 
singular behaviour of \M m +i A ...\ 2 leads to the soft and collinear divergences of the NLO 
contribution J m+1 da R in Eq. (|2.3|) if the phase-space integration over pj is performed in 
four dimensions. 

The starting point of our method for constructing the counter-term da A in Eq. ( |2.5| ) 
is the observation that the singular behaviour of |A^ m+ i ia ...| 2 is universal, that is, it is not 
dependent on the very detailed structure of M. m +x,a... itself. The origin of this universality is 
in the fact that, for its singular terms with respect to the momentum pj, the tree amplitude 
•M-m+x,a... can always be considered as being obtained by the insertion of the parton j over 
all the possible external legs of a tree- level amplitude Ai m ,a... with m QCD partons (Fig. 1). 
Thus, the singular behaviour of M. m +i, a ... is essentially factorizable with respect to M. m ,a... 
and the singular factor only depends on the momenta and quantum numbers of the QCD 
partons in M. m ,a...- Actually, according to the external- leg insertion sketched in Fig. 1, it 
is evident that the singular factor we are looking for is quasi-local, in the sense that it only 
depends on the momenta and quantum numbers of three partons: the parton j that is 
inserted onto M. m ,a... and the partons i and k in A4 m ,a...- 

This feature of the soft and collinear singularities will be used to obtain factorization 
formulae with the following symbolic structure (Fig. 2) 

\M 

m+l,a... | 2 \M m ,a...V ® V ijik • (4.1) 

Here Vij t k is the singular factor, which depends on the momenta and quantum numbers 
of the three partons i,j,k. As explained in detail in Sect. two of these partons (e.g. 
i and j) will play the role of 'emitter' and the third parton (e.g. k) that of 'spectator'. 




j m+1 

Figure 1: Diagrammatic representation of the external- leg insertion rule. The blobs denote 
the tree-level matrix elements and their complex conjugate. The dots on the right-hand 
side stand for non-singular terms both in the soft and collinear limits. 
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Figure 2: Pictorial representation of the dipole factorization procedure. When the partons 
i and j become soft and/or collinear, the singularities are factorized into the term Vy- fc (the 
dashed box on the right-hand side) which embodies correlations with a single additional 
parton k. 



Because of this structure, the factorization formulae described in Sect. |^ will be called 
dipole factorization formulae. In order to explicitly construct these formulae, we should 
first recall the known properties |7|j25[| of the tree-level QCD matrix elements in the soft 
and collinear limits. 



4.2 Soft limit 



Let us consider the matrix element M. m +\,a... an d the corresponding vector (in colour + 
helicity space) \l,....,m+ 1; a, ... > m +i, a .. defined in Eq. ( 13.111 ). Let us denote by pj the 
momentum of a final-state gluon in M. m +i,a...- I n the soft limit, which we parametrize in 
terms of an arbitrary four vector g M and a scale parameter A: 

■pfj = \q» , A - , (4.2) 

the matrix element squared behaves as follows |7]] 

m+l,a.. 

< l,....,m+ l;a,...\\l,....,m+ l;a, ... > m+ i, a ..— > 
-— 4nfi 2e a s m, a ..< 1, ...,m + 1; a, ..| [J^g)] 1 J^(q)\l, ....,m+ 1; a, ... > m , a .. -(4.3) 

Here we have neglected all the contributions less singular than 1/A 2 . The m-parton matrix 
element on the right-hand side of Eq. (O) is obtained from M. m +\,a... by simply removing 



the soft gluon pj. The term J^(q) is the eikonal current for the emission of the soft gluon q. 
Its explicit expression is given in terms of the momenta and colour charges of the partons 
in |l,....,m+ l;a, ... > m ,a..- 

J»(q)=Y, T i — + — + ■ ( 4 - 4 ) 

i Pi-q Pa-q 



The formula in Eq. (|4.3|) is well known. Here we limit ourselves to recalling a few 
properties of Eq. ( |4.3j ), which are relevant for understanding the structure of the improved 
factorization formulae that we shall introduce in Sect. |5|. 
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Although first derived in the four-dimensional case, Eq. (|4.3f) is actually valid in any 
number d = 4 — 2e of space-time dimensions. The only dependence on d is in the overall 
scale factor /z 2e on the right-hand side. 

The m-parton matrix element is not exactly factorized. Since the eikonal current in 
Eq. fl4.4|) depends on the colour charges of the hard partons, it leads to colour correlations 
on the right-hand of Eq. 



In the actual calculation of cross sections, Eq. ( f4.3| ) cannot be used as a true factorization 
formula not only because of these colour correlations. In fact, the tree-level matrix elements 
are unambiguously defined only when momentum conservation is fulfilled exactly. Since, in 
general, the m + 1- parton phase space does not factorize into an m-parton times a single- 
parton phase space, the right-hand side of Eq. (|4.3| ) is unequivocally defined only in the 
strict soft limit A = 0. Away from the point A = 0, care has to be taken in implementing 
momentum conservation. 



The form of the eikonal current in Eq. ( f4.4|) is actually valid both for massless and 
massive partons. Squaring the current as in Eq. fl4.3|) , and taking the massless limit, one 
obtains: 



+...) . (4.5) 



(p a q)(piq) 



Each term PkPi/(Pkq)(Piq) on the right-hand side of Eq. (|4.5|) leads to collinear singularities 
when the soft momentum q is parallel either to pi or to pt or p a . These collinear singularities 
can be disentangled by using the following identity 

PkPi PkP l PkPi (AR . 

+ 7 ; \ — , ( 4 -6) 



(Pkq) (piq) Pk ■ q (pi + Pk) ■ q v%-q {v% + Pk) ■ q 

and likewise for the terms PaPi/(p a q)(.Piq) , ■ ■ ■■ Thus, Eq. ( f4.3|) can be rewritten as follows 

m +i,a..< l,....,m + l;a, ...| |1, m + 1; a, ... > OT +i,a..-> 8nfi 2e a s (4.7) 

1 , -i . , ,T k -Tip k pi 
• > — > m a ..< 1, ...,m + 1; a, .. — ■ - — 1, ....,m + 1; a, ... > m a .. + . . . , 

i Ptq t+i [Pi+Pk)q 

where the dots stand for similar contributions that involve the initial-state parton a, . . .. 
The dipole structure mentioned at the end of Sect. [44] starts to emerge from Eq. (|4.7|) . 
Each term on the right-hand side of Eq. (|4.7|) depends on the radiated soft momentum q, on 
the 'emitter' momentum pi (whose direction signals the presence of a collinear singularity) 
and on the 'spectator' parton k (which accounts for colour correlations). 



4.3 Collinear limit 

Let us consider the momenta Pi and pj of two final-state partons in M. m +i, a ...- The limit 
where Pi and pj become collinear is precisely defined as follows 

pi = ZJ f + k»- , = (1 - z)p» -e ± - 



z 2pn 3 1 — z 2pn 
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In Eq. (|4.8| ) the light-like (p 2 = 0) vector p M denotes the collinear direction, while n fl is 
an auxiliary light-like vector which is necessary to specify the transverse component k± 
{k\ < 0) (k±p = kj_n = 0) or, equivalently, how the collinear direction is approached. In 
the small-/c_i_-limit (i.e. neglecting terms that are less singular than l/k±), the m + 1-parton 



matrix element behaves as follows p5 | 



m+l,a.. 

< 1, ....,m + l;a, ...||1, ....,m+ l;a, ... > m+ i, a ..— > 

— 47r/i 2e a s m ,a..< 1, ...,m + 1; a, ..| Puj),i{z, k±; e) |1, ....,m + 1; a, ... > m>a .. .(4.9) 
PiPj 

The m-parton matrix element on the right-hand side of Eq. ( [4.9|) is obtained by replacing 
the partons i and j in -M m+ i iCl ... with a single parton denoted by ij. This parton carries the 
quantum numbers of the pair i + j in the collinear limit. In other words, its momentum is 
p M and its other quantum numbers (flavour, colour) are obtained according to the following 
rule: anything + gluon gives anything and quark + antiquark gives gluon. 



The kernel Pi 



in Eq. 



is the d-dimensional Altarelli-Parisi splitting function. It 



depends not only on the momentum fraction z involved in the collinear splitting ij — > i + j, 
but also on the transverse momentum k±_ and on the helicity of the parton ij in the m- 
parton matrix element. More precisely, Puj) t i is a matrix acting on the spin indices of the 
parton ij in m>0 ..< 1, m + 1; a, ..\ and |1, ....,m+ 1; a, ... > m ,a..- Because of these spin 
correlations, the square of the m-parton matrix element cannot be simply factorized on the 
right-hand side of Eq. 



(4.10) 



The explicit expressions of P a b{z, k±; e) for the splitting processes 

a(p) -> b(zp +k ± + 0(kD) + c((l - z)p -k ± + 0(k 2 ± )) 



are as follows 



< s\P qg (z, k±; e)\s' >= 5 SS , C F 



1 + z 2 



6(1 - Z) 



(4.11) 



< s\P qg (z, k \ ; e)|s' >= 8 SS > Cf 



ez 



(4.12) 



< n\P gq {z, k±;e)\u >= T } 



R 



< fi\P gg (z, kj_; e)\v >= 2C A 



-9 



4,(1 - zf-j^ 



1 - z 



+ 



l-z 



2(l-e)z(l -z 



rv J_ >v | 



(4.13) 



(4.14) 



where the spin indices of the parent parton a have been denoted by s, s' if a is a fermion 
and /i, v if a is a gluon. 
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Equations ( 4.11|) - (|4.14 ) lead to the more familiar form of the (i-dimensional splitting 
functions only after average over the polarizations of the parton a. The <i-dimensional 
average is obtained by means of the factors 



1 



for a fermion, and (the gauge terms are proportional either to p^ or to p" 



(4.15) 



d-2 



d^u(p) 



2(1 -e) 



-g^u + gauge terms) 



with 



g txu d flv (p) = d-2, p»d flu (p) = 



(4.16) 



(4.17) 



for a gluon with on-shell momentum p. Denoting by < P a b > the average of P a b over the 
polarizations of the parton a, we have: 



< P qq (z;e) >=C f 



1 + z 2 
1 - z 



e(l-z) 



(4.18) 



< P qg (z;e) >= C F 



i + a-z) 



— ez 



(4.19) 



< Pgq{z\e) >'- 



<P gg (z;e) >= 2C A 



2z(l - z) 



1-6 

1 - Z 



1 - z 



+ 



z 



+ z(l-z) 



(4.20) 
(4.21) 



So far, we have considered the case in which two final-state partons in M. m +i,a... become 
collinear. In general, one has to deal also with the case in which a final-state parton % 
becomes collinear to an initial-state parton a. Here the collinear limit is defined as follows 



=(l-z)p£ + A£ 
k 2 



n 1 



1 — x 2p a n 







x 



(4.22) 



and the corresponding splitting process a — > ai +i involves the transition from the initial- 
state parton a to the initial-state parton ai with the associated emission of the final- 
state parton i. The quantum numbers of the parton ai are assigned according to their 
conservation at the QCD tree-level vertices: if a and i are partons of the same species, 
then ai is a gluon; if a is a fermion (gluon) and i is a gluon (fermion) then ai is a fermion 
(antifermion). 



The analogue of Eq. (19) in the collinear limit ( 4. 22]) is the following 



m+l 



)0 ..< 1, ....,m + l;a, ...\\1, ....,m+ l;a, ... > 



m+l,a.. 



(4.23) 



1 1 

X PiPa 



Atc/j, e a s m ,ai..< l,...,m + l;ai,..\ P a r a{) (x,k ± ;e) |1, m + 1; ai, ... > 
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Now, the m-parton matrix element on the right-hand side is obtained from M. m +\,a... by 
removing the final-state parton i and replacing the initial-state parton a with the parton 
ai. Note two main differences with respect to Eq. ( |4.9|) : on the right-hand side of Eq. (|4.23| ) 
there is an additional factor of 1/x and the initial-state parton ai carries the momentum 

xi 

As in the case of Eq. ( |4.3|) , Eqs. (|4.9|) and (|4.23| ) have to be regarded as limiting formulae 
rather than factorization formulae. Their implementation in the calculation of QCD cross 
sections indeed requires a careful treatment of momentum conservation away from the 
collinear limit. 

Note that the splitting functions in Eqs. ( [4.11| - P~bi| ) are divergent for z — *> 0, 1. These 
divergences are the soft singularities already discussed in Sect. [O. When using Eqs. (|4.7| ) 



and ( |4.9|) to approximate the singular behaviour of M. m +i,a... care has to be taken in order 
to avoid double counting the soft and collinear divergences in their overlapping region. 



Note also that Eqs. (|4.9| ) and ( |4.23| ) do not depend solely on the collinear momenta 



Pi,Pj and Pi, p a - In fact, the Altarelli-Parisi splitting functions produce spin correlations 
with respect to the directions of the other momenta in the matrix element. The dipole 
structure of these limiting formulae is thus hidden in the azimuthal dependence of these 
correlations. 
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5 Dipole factorization formulae 



In this Section we present in detail our basic formalism to construct the local counter-term 
da A for the NLO cross section in Eq. (|2.4J) . We introduce improved factorization formulae 
for the QCD matrix elements. These formulae are based on a dipole structure with respect 
to colour and spin indices and have the following main features. Our dipole factorization 
formulae coincide with Eqs. ( |4.7|) and ( |4.9| ) (or ( |4.23| ) ) respectively in the soft and collinear 
limit. These limits are approached smoothly, thus avoiding double counting of overlapping 
soft and collinear singularities. This smooth transition is possible because the dipole formu- 
lae fulfil exact momentum conservation. Actually, we present several alternative versions of 
the factorization formulae that differ from one another in the implementation of momentum 
conservation away from the soft and collinear limits. These differences are then used to 
match the phase-space constraints that are encountered in the calculation of different kinds 
of QCD cross sections. In this manner, for any QCD process, we achieve the analytical 
integrability of the counter-term da A over the single-parton subspace leading to soft and 
collinear divergences. 

We start by considering the case of final-state singularities for matrix elements without 
(Sect. |5.1|) or with (Sect. |5.2| ) initial-state partons. Sections |5]3| and [53J deal with initial- 
state singularities in the case of one or two initial-state partons, respectively. Cross sections 
with identified hadrons in the final state require the introduction of fragmentation functions. 



Factorization formulae suitable for these fragmentation processes are presented in Sects. |5T4 



and 5.6 



The main properties of the dipole factorization formulae are considered in detail in 
Sect. |5.1j. In the following Subsections we limit ourselves to writing down the formalism 



and emphasizing the relevant differences with respect to the case discussed in Sect. 5.1 



5.1 Final-state singularities with no initial-state partons 

The dipole factorization formula in the limit Pi ■ pj — > for the matrix elements with no 
partons in the initial state is the following 

m+ i< 1, ....,m + 1||1, ....,m + 1 > m+ i= ^2 X>y,fc(pi, ...,p m +i) + . . ■ (5.1) 

where . . . stands for terms that are not singular in the limit pt • pj — > and the dipole 
contribution is given by 

Vij,k (Pl> -,Pm+l) = -7T~ ( 5 ' 2 ) 

2pi • Pj 

~ ~ Tfc • T ■ ■ ~ ~ 

■ m< 1) ••) ij ') ••) k, .., ?7i + 1| — 2 V ij t k |1, - , ••, k, .., m + 1 > m . 

ij 

The m-parton matrix element on the right-hand side of Eq. (|5.2j) is obtained from the 
original m + 1-parton matrix element by replacing a) the partons % and j with a single 
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Figure 3: Effective diagrams for the different dipole formulae introduced in Sect. [5| The 
blobs denote the m-parton matrix element. Incoming and outgoing lines respectively stand 
for initial-state and final-state partons. 

parton ij (the emitter) and b) the parton k with the parton k (the spectator)^. 

All the quantum numbers except momenta are assigned as follows. The spectator parton 
k has the same quantum numbers as k. The quantum numbers of the emitter parton ij 
are obtained according to their conservation in the collinear splitting process ij — > i + j 
(cfr. Sect. |4.3| ). This rule applies to Eq. ( p.2| ) as well as to all the dipole formulae we shall 
introduce in the following Subsections. 

The momenta of the emitter and the spectator are defined in different ways in different 
dipole formulae. In Eq. ( |5.2| ) we have 

Pk = T—~ P k ' Pij =Pi+Pj ~ I V ^ Pk , ( 5 - 3 ) 

J- Vij,k J- Vij,k 

*In general, we use the following notation in the dipole formulae (Fig. 3). A pair of indices like ij 
denotes the emitter parton and a single index like k denotes the spectator parton. If these indices appear 
as subscripts or superscripts, they respectively indicate final-state or initial-state partons. 



23 



where the dimensionless variable y^^ is given by 

PiPj 



Vij,k 



PiPj + PjPk + PkPi 



(5.4) 



Note that both the emitter and the spectator are on-shell (pfj = fr k = 0) and that, 
performing the replacement {i,j,k} —>■ {ij,k}, momentum conservation is implemented 
exactly: 

Pi + Pj +Pk= Pij + Pk ■ ( 5 - 5 ) 
These are common features of all the dipole formulae in the paper. 

In the bra-ket on the right-hand side of Eq. ( |5.2|) , T k and Tjj are the colour charges of 
the emitter and the spectator and Vij t k are matrices in the helicity space of the emitter. 
These splitting matrices, which depend on and on the kinematic variables Zi,Zf 



PiPk 



PiPk 



, Zj 



PjPk 



PjPk 



1 — Zi 



(5.6) 



PiPk + PiPk PijPk PjPk + PiPk PijPk 

are related to the d- dimensional Altarelli-Parisi splitting functions in Eqs. (|4.11|) - (|4.14 ). 



For fermion + gluon splitting we have (s and s 1 are the spin indices of the fermion ij 
in < and > respectively) 



< s\V qigj>k (zi;y ijik )\s' > 



^qigj,k &ss' • 



1 - Zi(l - y ij>k ) 



- (1 + ~z t ) - e(l - Zi) 



(5.7) 



For quark + antiquark and gluon + gluon splitting we have (/x and v are the spin indices 
of the gluon ij in < and .. > respectively) 



< fi | V qiq k (zi)\u >= 8vr/x e a s Tj 



PiPj 



[ziPi ~ ZjVfj (ziPi ~ ~ZjP V j) 



qiij,k ' 



(5.8) 



< ^\Vg igj ,k(zi;yij,k)W >= 167r/x e a s C A 
1 _\ . 1 , 



+ 



1 ~Zj(l -Vij,k) 



-2 + 1-e 



PiPj 



\1 - Zi(l - y ij;k ) 

ZiPi - ZjP%) (ZiPi - Zjrf) 



yin, 

9i9j,k 



(5.9) 



Soft and collinear limits 

Note that the matrices Vij, k do not lead to two-particle singularities in any of the limits 
Pi ' PjiPi ' PkiPj ■ Pk — > 0. This is because the only non-polynomial dependence on Zj is in 
the factors 

1 = PiPj + PjPk + PkPi , 5 1Q x 

1 - Zi(l - y ijik ) {Pi + Pk)Pj 

Therefore the dipole term on the right-hand side of Eq. ( |5.1| ) contains only collinear and 
soft divergences for pi ■ pj — >• 0. 
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In the collinear limit of Eq. (|4.8| ), the dipole variables become: 

kj 

Vthk ~* 2z(l - z) PPk ' Zl ~ h ^ Z ' 

% - Pk, P%^P» ■ (5.H) 

Thus, the m-parton matrix element in Eq. ( |5.2|) coincides with that in Eq. (|4.9| ). Moreover, 
the vector — ZjPj in Eqs. ( |5.8| , p.9D becomes [ziZ — Zj(l — z)]p^ + k^_ and, since its 
p^-component gives a vanishing contribution to Eq. ( |5.2j) (p M < fi\\..,p,.. > m = because 
of gauge invariance), the matrices Vijk become proportional to Altarelli-Parisi splitting 
functions: 

V ij;k -> 87r/i 2e a s P(ij),i( z > ki; e ) • ( 5 - 12 ) 
In particular, the only dependence on k that survives in Eq. (|5.2|) is that on T k . Therefore, 
one can perform the sum over the colour charges on the right-hand side of Eq. (|5.1|) and, 
using charge conservation (J2k^i,j Tk = —Tij), one can check that this equation reproduces 
the collinear behaviour in Eq. ( f4.9|) . 



In the soft limit of Eq. ( f4.2|) , we have 

Vij,k > , Z% ► 1 , Zj > , 

Pk ^ Pk , Pij Pi ■ ( 5 - 13 ) 

Thus the m-parton matrix element in Eq. ( |5.2| ) coincides with that in Eq. (|4.7| ). Moreover, 
the only singular factor in is due to the term in Eq. (|5.10|) : 



' ' Vm (5.14) 



1 - Zi(l -yijjk) A {pi + Pk)q 
which gives 

XV igj , k 167r/i 2e a s T 2 . (5.15) 

Inserting this expression into Eqs. ( |5.2|) , we see that the dipole term in Eq. ( |5.1| ) reproduces 
the soft limit in Eq. (pL"7D . 

This discussion proves that the dipole formula ( p.l| ) provides a point- wise approximation 
of the m + 1-parton matrix element in the singular region Pi - pj — > 0. Note that to achieve 
this, the helicity dependence of the splitting kernels is essential. The azimuthal 

correlations due to this dependence cancel after integration over Pi (see Eqs. ( |5.26| , |5.27D ) 
and hence, provided that the counting of helicity states is consistently performed in d 
dimensions, they are not relevant for reproducing the correct poles in 1/e in the contribution 
J m+ ida A to the NLO cross section in Eq. (|2.5|) . Nonetheless, these correlations have to 
be properly taken into account in constructing the local counter-term do~ A , which makes 
the contribution [da R — da A ] integrable in four dimensions. Indeed, the parton azimuthal 
correlations due to this dependence are not only essential in the most general case when 
the jet cross section explicitly depends on them, but even when it does notQ 

^In this case the evaluation of J m+1 da R in four dimensions usually involves double angular integrals of 

the type dcosd f Q ™ dip cosip/(l — cos0), where <p is the azimuthal angle. These integrals are mathe- 
matically ill-defined. If their numerical integration is attempted, one can obtain any answer whatsoever, 
depending on the detail of the integration procedure. Performing the integral analytically before going to 
4 dimensions, one obtains dcosQ dip cos (p/(l — cos6>) sin~ 2£ 6>sin~ 2e ip = 0. 
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Phase space factorization 

The definition ( |5.3| ) of the dipole momenta is particularly useful because it allows us to 
exactly factorize the phase space of the partons k into the dipole phase space times a 
single-parton contribution. Indeed, let us consider the following 3-particle contribution to 
the final-state phase space: 

d<j){p uPhPk ;Q) = ( 2 ^-i 5 +(Pj) 5 +iPk) (2n) d 5 {d) (Q-p t - Pj -p k ) . 

(5.16) 

In terms of the momenta pij, pk and Pi, this phase-space contribution takes the factorized 
form: 



t,Pj,Pk',Q) = d(j)(pij,Pk;Q) [dpi(pij,Pk)] , (5.17) 
where 

d d Pi 

[dpi(Pij,Pk)} = (o^y-l S +(Pi) >7(Pi>Pi3>Ph) > ( 5 - 18 ) 
and the Jacobian factor is 

\d-3 

1 - z. 

In terms of the kinematic variables defined earlier, we have 

[dpi(Pij,Pk)} = ^iq^I — (2^)1-2* dli dyij ' k ( i ^ 1 _ ^ (^. fc ( 1 ~ 

• {z i {l-~z i )r (l-^) 1_2£ % e fe , (5-20) 
where dV^ d ~^ is an element of solid angle perpendicular to p^ and pk and thus 



J(Pi;Pij,Pk) = ©(l - Zi) ©(l - Vij,k) (1 : ^ 3 • (5.19) 



f d Q(^) = J* . (5.21) 
J TKT(l-e) 



Integration of the splitting functions V 



ij,k 



To evaluate the integral / m+1 do A in Eq. (j2.5|), we can compute the integral of the 
splitting functions V^*. over [dpi(pij , pk)) once and for all. The only non-trivial point 
involved in this integration regards azimuthal correlations. 



Note, however, that the spin correlation tensor 

1 



PiPj 



ZiPi ~ ZjPj) 



(5.22) 



in Eqs. ( |5.8| , |5.9| ) is orthogonal to both p^ and p\y Using this property and Lorentz covari- 
ance (the integral may depend only on p^ and pk), it follows that the azimuthal integration 
of the spin correlation tensor gives a contribution of the type 



K, 



A 



Pij Pk + VkVij 
Pij ■ Pk 



(5.23) 
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On the other hand, the (on-shell) matrix element |1, .., ij, .., k, .., m + 1 > m is conserved 
(gauge invariance), that is, 



Pii < A*||l, k , ..,m+ 1 > m = 



(5.24) 



Thus, only the term —A in contributes to the dipole formula after integration over 
[dpi(pij ,Pk)]- This term can be singled out by taking the <i-dimensional average over the 
polarizations of the emitter: 



d-2 



1 



d-2 



(5.25) 



In conclusion, after integration of the dipole V^ipi, ■■■,Pm+l) over [dpi(pij,pk)], only colour 
correlations survive, in the form: 



[dpiipij^k)] V ijik (pi, ■■■ yPm+1, 



ij,k m 



< l,..,ij,..,k,..,m + l 



\l,..,ij,..,k,..,m + 1 >, 



where 



V^ = j[dp t {p^p k )] ^ <V lhk >= g^^f^) Me) 



\2pijPkJ 



(5.26) 



(5.27) 



and < Vij t k > denotes the average of Vij,k over the polarizations of the emitter parton ij. 

The function VV,(e) depends only on the flavour indices i and j. Using Eq. Q5.20| ), from 
the definition of V^e) in Eq. (|5.27|) we obtain 



Vij(e) = / dzi (zi(l - Zi))~ 



dy 

o y 



l-2e - e < Vij,k( z ^y) > 



57TQ;g/i 



2c 



where the spin-averaged splitting functions are: 



2c 



1-2(1-3/) 



[l + z)-e(l-z) 



{5.21 



(5.29) 



< V q q >k (z;y) > 



2c 



T, 



R 



1 - 



< V g g !k (z;y) > 



2c 



2C> 



+ 



25(1 - z) 
1 - e 

1 



l-z(l-y) l-(l-2)(l-y) 



2 + z(l-z) 



Performing the integration in Eq. ( |5.28|) , we find 



V 95 (e) 



r 3 (i-e) 

r(l-3e) 
T 3 (l - e 



1 1 3 + e 



e 2(1 - 3e) 



2e 



R 



1 



2(1 -el 



V 59 (6) 



r 3 (i-6 
r(i -3t 



2C, 



e (l-3e)(3-2e) 
1 1 11 - 7e 
e2 + 7 2(l-3e)(3-2e) 



~-2C A 



16 



3e 9 



+ 0(e) 



(5.30) 
(5.31) 

(5.32) 
(5.33) 



1 11 50 7T 2 ^. , 



(5.34) 

Note that in each case, the first result is exact in any number of dimensions, d = 4 — 2e, 
while the latter is valid for e — > 0. 
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5.2 Final-state singularities with initial-state partons 



In the presence of initial-state partons a, . . ., the m + 1-parton matrix element has both 
final-state (pi ■ pj — > 0) and initial-state (p a ■ pj — > 0) singularities. Let us first consider 
the dipole factorization formula for final-state singularities. Neglecting terms that are not 
singular when pi ■ pj — > 0, we have: 

m+i, a ..< 1, m + 1; a, ...||1, m + 1; a, ... > m +i,a.. = ^ T^ij,k(Pi, •••> 



+ 



+ 



(5.35) 



The first term on the right-hand side is the same contribution as in Eqs. ( |5.1| , |5.2|) , while the 
terms in the square bracket represent additional dipole contributions in which the emitter 
is the final-state parton ij and the spectators are the initial-state partons a, . . . . These 
dipole terms are explicitly given by 



V^ j (p 1 ,...,p m+1 ;p a ,..) 



1 1 



.< 1, -Jj, ..,m+l 



^Pi ' Pj %ij,a 

T T 

■a I — 13 V a . II 

ij 



(5.36) 



,ij,..,m+ l;a, ... > m . 



In Eq. ( |5.36|) the momenta of the spectator a and the emitter ij are defined as follows 

Pa = Xij,aPa , Pij = Pi + Pj ~ (1 ~ »ij,a) Pa , (5-37) 

_ PiPa + PjPa ~ PiPj 
[Pi+Pj)Pa 

and the corresponding splitting functions are: 

2 



< s\V a q . g .{zi;Xi j>a )\s' > = 8ir/j, 2e a s C F 



1 Zi ~\~ ( 1 •Eij,a) 



1 + Zi) - e(l - 2j 



(5.39) 



+ 



~0 



1 £j ~H (1 %ij,a) 



i - 5, + (i - x y , a ) - 2 ) + (1 - e) ^~ ~ ^ ~ ^ 



< yU |VV(z i )|z/>=87r^a s T J 



P;Pj 



^ipf - Z$) (Zip- - ZjPj) 



where 



PiPa 



PiPa 



PiPa + PjPa PijPa 



, Zj 



PjPa 



= PjPa 

PjPa + PiPa PijPa 



1 — Zi 



(5.40) 
(5.41) 

(5.42) 
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In the definition of the splitting functions V"-, the kinematic variable plays a role 
similar to that of y^ ^ in Vijk'- it provides a smooth interpolation between the soft and 
collinear limits. Following the same argument as in Sect. |5J], it is straightforward to check 
that Eqs. ( |5.35| , [5.36] ) reproduce the correct soft and collinear behaviour of the m + 1-parton 
matrix element in the limit Pi ■ Pj — > 0. 

Note that, comparing Eqs. ( |5.39| - |5T4T|) with the equivalent expressions when the specta- 
tor is a final-state parton, Eqs. (|5.7H5.9|) , they are not symmetric under crossing symmetry, 
Pk —>■ —p a - While the non-soft terms do become identical in the collinear limit, the soft 
terms do not in the soft limit. This is because the crossing- symmetric eikonal term is split 
between the two dipole terms that, separately, are not crossing symmetric. In fact, after 
adding Eqs. ( |5.39| - |5T4~l| ) and the equivalent expressions for initial-state singularities given 
below in Eqs. (|5 . 65H5768|) , the soft terms for pj — ► do obey crossing: 

PkPi PkPi PkPi 

Final-state parton p k : H = — , (5.43) 

PiPAPiPj + PkPj) PkPjiPiPj + PkPj) [PiPjKPkPj) 

T ... , , , , PaPi PaPi PaPi / c A ,\ 

initial-state parton p a : H = — . (5.44) 

PiPjiPiPj + PaPj) PaPjiPiPj + PaPj) {PiPj){PaPj) 

If, instead, naive crossing, p k — > —p a , had been used on the individual terms on the left- 
hand side of Eqs. ( |5.43| , |5.44"D , spurious singularities would develop at piPj = p a Pj- Similar 
comments concerning crossing apply to the formulae in all later Subsections. 

Considering the three-parton phase space 

d<j>( Pi , Pj ;Q +Pa) = J0-^4P1) J^bl^) (2 7 r) d ^(g + p a -p,-p J ), (5.45) 

the analogue of the phase-space factorization in Eq. ( |5.17| ) is the following phase-space 
convolution 

d(f)(pi, pj; Q +Pa) = [ dx d(f)(pij; Q + xp a ) [dpi(pij;p a , x)] , (5.46) 
J o 

where 

d d pi 1 

[dpi(Pij]p a ,x)} = ' hiPi) Q(x)Q(l -x) 5{x-Xi j>a )- , (5.47) 

\Lix) a 1 1 — Zi 

or, more explicitly, using the kinematic variables in Eqs. (|5.38| , |5~42"1) : 

[dpi(Pij]Pa,x)} = ^ J g°] — (S) 1 - 26 ^ dXij ' a ~ - x)) 

■ (zi(l - ~ Zl )Y e 5(x - x ij>a ) (1 - x)~ e , (5.48) 

where dVL^ d ~^ is an element of solid angle perpendicular to pij and p a . Note that we are 
interested in using Eqs. (|5.46|J5.47|) for a NLO calculation whose leading-order kinematic is 
Q + p a = P with p 2 = 0,po > 0. Thus the physical constraint Q 2 < is always understood 
in Eqs. (|5.46|J5.47|) . 

The integration of the splitting function V°- over the phase space in Eq. (|5.47 ) defines 
the functions Vy(x;e) (the treatment of the azimuthal correlations is similar to that in 
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Sect. [5.1|) : 
a s 1 



27rr(l-e) 
that is, 



Vij{x;e) = / [dpi(jpij]p a ,x)] 



> 



V i:i (x; e) = e(x)6(l-x) (j 



l+e „i 



< V%(z^x) > 



(5.49) 



(5.50) 



As usual, < Vfj > denotes the average of V?- over the polarizations of the parton ij. 

The e-dependence of Vy(a;;e) has to be interpreted with care. For e = 0, Vij(x;e) is 
well-defined everywhere except the point x = 1 where it is singular as MizfEi, T n order 
to avoid this singularity, the limit e — > has to be taken uniformly in x. In this manner, 
Vij(x;e) defines an x-distribution whose coefficients contain poles in 1/e. More precisely, 
we have: 

V^x; e) = [ Vi,(x; e) ] + 5(1 - x) [' dz V tJ (z; e) , (5.51) 



[V ij (x;e)] + =[V ij (x;e = 0)] + + O(e) , (5.52) 
where the '-(-'-distribution is defined, as usual, by its action on a generic test function g(x): 



dxg(x) [ V(x) ]_ 



dx[g{x)-g{l)]V{x) . 



(5.53) 



Note that in Eq. (|5.52|) , the expansion around e = is well defined, i.e. the 0(e) term on 
the right-hand side is integrable in x. 



The explicit form of the spin-averaged splitting functions < V?- > is the following: 

2 



<V a gg (z;x)> 



<V a gg (z;x)> 



2c 



c f 



l-z+(l-x) 



(l + 2)-e(l-5) 



2c 



2C/ 



-R 



+ - 



25(1 - z) 
1 -e 

1 



1 — z + (1 — x) z + (1 — X 
Inserting Eqs. ( |5.54| - |5~5lj| ) into Eqs. ( |5.50| - |5~5*2] ) we obtain: 

2 h, 1 



2 + 5(1-5) 



(5.54) 

(5.55) 
(5.56) 



V qg (x;e) = C F 

+ 5(1 -x 
2 



1 — x 1 — x 



V 95 (e) 



2\l-x 
+ 0(e) , 



+ i 

+ 1 — x 



ln(2 - x) 



V qq (x-e) = -T RKi _ x 



V gg (x-e) = 2C A 



In 



+ 5(l-x) 
1 



V qq (e) + -T R 



1 — x 1 — x 



11 



1 



+ 5(1 -x) 



6 VI — x 
+ 0(e) , 



+ 



+ 0(e) 
2 



(5.57) 
(5.58) 



1 — x 



ln(2 - x) 



(5.59) 



where Vy(e) are the functions in Eqs. ( |5.32| - |5T34"|) . Note that the terms 5(1 — x)Vij(e) 
account for all the e-poles of Vij(x; e), i.e. Vij(x; e) — 5(1 — x) Vy(e) are finite for e — ► 0. 
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5.3 Initial-state singularities: one initial-state parton 



In the limit p a -Pi — > 0, the dipole factorization formula for the m + l-parton matrix element 
with a single initial-state parton a is the following 

m +i,a< 1, ....,m+ l;a||l, ....,m+ l;a > m +i, a = X^fc l (Pi» ...,p m +i;p a ) + . . . , (5.60) 

k^i 

where the dots stand for non-singular contributions and the first term is the sum of the 
dipole contributions in which the emitter is the initial-state parton ai and the spectator is 
the final-state parton k. These dipoles are given by 



^k(Pl,--,Pm+l,Pa) = 



,„< l,..,k,.,m+ l;ai 



1 



2j?a ' Pi ^ik,a 



—^3 — y fe |1, ..,/c, ..,m + l;az > m , a 



(5.61) 



Note that, unlike the case in which the emitter is a final state ij, now the dipole is not 
symmetric under a <-> i. In fact, the momentum of the emitter parton ai is parallel to p a : 



where 



%ik,a 



Pai ~ x ik,aPa j 



PkPa + PiPa ~ PiPk 



(Pk + Pi)Pa 

On the other hand, the momentum of the spectator k is not parallel to p k . 

#=p£+rf-(l-**,a)tf, 



(5.62) 
(5.63) 

(5.64) 



The splitting functions V£ l are matrices in the helicity space of the parton ai. Their 
explicit expression is the following 



< s\V q k a9i (x ik , a - Ui )\s' > = 87r/i 2e a s C F 5 S 

— v k °ss' 1 



1 x ik,a ~l" l^i 



'1 + x ik „) -e(l -x ik n) 



(5.65) 



< sIV^^iM)!*' >= 8 ^ «s T fl [1 - e - 2x ifc , a (l - x ikA )\ 8 SS , = V 9aqi 5 SS , , (5.66) 



1 - X ik , a 2Ui(l- Ui) ( pf 



+ 



A 1 
Pk 



9^ ^ik,a 



Pi 



%ik,a 



PiPk 



Ui 



Pk 



1 — Ui I \Ui 1 — Ui 



(5.67) 



< n\ V 9 k a9i {x iKa - Ui )\v >= 16vr/i 2e a s C A 



-9 



_ 1 x ik,a "i" U, 



1 -|- Xjfc ia (l Xifc^a) 



+ (1-6) 



1 - ^ifc,a «i(l - ^) / p£ 



•Eik,a 



PiPk 



Ui 



1 - 



PI 

■u,- 1 — u, 



(5.68) 
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where 



PiPa (r rn \ 

Ui = . (5.69) 

PiPa + PkPa 



Note that the factor l/xik, a on the right-hand side of Eq. ( |5.61|) is relevant for repro- 
ducing the collinear limit in Eq. (|4.23|) . 



The analogue of the phase-space convolution in Eq. ( |5.46 ) is the following identity 
(similarly to Eqs. ( |5.46| , |5.47| ), we have Q 2 < in the physical region of the phase space we 
are interested in) 



d(p(pi,Pk]Q + Pa) = / dx d(j)(p k ; Q + xpa) [dpi(p k ; p a , x)) 
Jo 



(5.70) 



where 



„ d d pi 1 

[dpi(p k ;p a ,x)} = " hipl) 8(x)8(l - x) 5(x - x ik>a )- . (5.71) 

[Zn J 1 — Ui 



Using the kinematic variables in Eqs. (|5.63|J5.69| ), the phase space in Eq. (|5.71| ) can be 
written as follows 



{dpi(p k ;Pa,x)} 



16^2 ( 27r )l-2e dUi dXik > a ~ U i)) - X)) 

(ui(l-Ui))' e 5{x-x ik , a ) (l-x)~ e , (5.72) 



where d^ d 3 - ) is an element of solid angle perpendicular to p k and p a . 

It is also useful to introduce the following integral of the splitting function V™: 

as 1 / 4:iru 2 



2irT(l - e) \2p k p a/ 



V a ' ai (x\e) 



[dpi(p k ;pa,x)} _}_ Us ^ < V a k \x ikt a]Ui) > , 
2p a Pi n s {a) 



(5.73) 

where n s (a) (n s (ai)) is the number of polarizations of the parton a (ai) (n s = 2 for fermions 
and n s = d — 2 = 2(1 — e) for gluons) and < V™ > denotes the average of V™ over the 
polarizations of the emitter parton ai. 



From the definition in Eq. (|5.73 ), we have 

1 



V a ' at (x; e) = e(x)e(l-x) 



1 — XJ Jo Ui 



,,IU Ml-^r^^^^f. (5.74) 



n s (a) 8irasu 



The integration over Ui in Eq. ( |5.74|) leads to e-poles of the type f(x)/e, where f(x) 
can be either an x-integrable function for x — 1 or (1 — x)" 1 . In order to get a series in e 
whose coefficients are well-defined x-distributions, we write: 



V a ' ai (x: e) 



i n 



e ix 



exV a ' ai (x:e) 



+ e5(l-x) / dz zV a > ai {z-e) 



(5.75) 



The two terms in the curly bracket can separately be expanded in e and lead to contributions 
with the following structure 



P"-*"(.r) const. (- + 0(1)) Ml - .<) +C)(c) . 



(5.76) 
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where P a,at (x) is the four- dimensional regularized (real + virtual) Altarelli-Parisi proba- 
bility for the splitting process a(p a ) — > ai(xp a ) + z((l — x)p a ). 

In order to obtain this result, we first compute the spin- averaged splitting functions 
<Vf(x;u)>: 



n s (q)<V q k 9 (x;u)> 



n s (q) 87ra s /x 



2c 



c f 



1 — x + u 



(1 + x) - e(l - x) 



(5.77) 



n, 50 < W x) > 



n s (g) 87ra s /i 2e 

n s (g) < V 9 k 9 (x;u) > 
n s (g) 8na s n 2e 



T, 



R 



2C, 



1 — X + U X 



2x{l-x) 
1 -e 

H l + x(l-x) 



n s (g) < Vf(x) > 



C f 



(1 -e)x + 2 



1 — x 



X 



n s (q) 87ra s /i 2e 

Then, inserting Eqs. (|5.7?H5T80|) into Eqs. ( |5.74| , |5.75"D , we find: 

V q ' 9 (x; e) = - X - P q9 (x) + P q9 (x) ln(l - x) + C F x + 0(e) , 

V 9 ' q {x- e) = — P OT (x) + P 9q {x) ln(l - x) + T R 2x(l - x) + 0(e) 



(5.78) 
(5.79) 

(5.80) 

(5.81) 
(5.82) 



V q > q (x;e) 



- P qq (x) + 5(1 - x) 
e 



+ Ci 



4 m. 1 



1 — x 1 — Xj 4- 1 — X 



V«(e) + ( 3^ - 5j C F 
ln(2 - x) 



+ (1 - x) - (1 + x)ln(l -x) 



+ 0(e) , 



(5.83) 



V 9 ' 9 (x;e) = — P"(x) + 5(1 - x) 



+ ™N f T R 



4 In- 1 



1 — x 1 — x/ + 1 — X 
1 — X 



9 

ln(2 - x) 



+ 2(-l + x(l-x) + -)ln(l-x) 



x J 



0(e) 



(5.84) 



where the functions V^-(e) are given in Eqs. ( |5.32| - 5~3^ ) and the (customary) regularized 
Altarelli-Parisi probabilities are 



P19( X ) = C F 



1 + (1 -x) 



X 



(5.85) 



psi( x ) =t r x 2 + (1 - x 



(5.86) 
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pii( x ) = c F 



l + x 2 
1 — X 



(5.87) 



P"(x) = 2C, 



1 — x 



1 — X 



X 



— 1 + x(l — x) 



+6(l-x)(^C A -^N f T R ) . (5.88) 



To simplify our notation in the following Subsections, it is also convenient to introduce 
the regular (i.e. not singular for x — > 1) parts P^ g (x) of the Altarelli-Parisi probabilities as 
follows 



P T %(x) = P ab {x) - 8' 



ab 



2 <rb) 



(5.89) 



3 r U r 2 t at 



(5.90) 



that is, 



P*{x) = P a \x) 



yab I 



if a 7^ b 



P t %(x) = -C F (l + x) , P r %(x)=2C A 



1 — X 



1 + x(l — x) 



x 



(5.91) 



5.4 Dipole formulae with one identified parton in the final state 



The perturbative QCD treatment of cross sections that, in addition to jets, have identified 
particles in the final state requires the introduction of fragmentation functions. In order to 
deal with the collinear singularities related to the fragmentation process, it is convenient 
to consider dipole factorization formulae that treat the final-state partons in an unsym- 
metric manner. Thus, the factorization formulae we are going to present below rescale the 
momenta of identified and non- identified partons in a different way. 

Let us consider the matrix element \a, 1, m + 1 > m+ i +a , where a denotes a parton 
whose momentum p a is identified in the final state. In this case we can get singularities 
from the regions where Pi ■ p 3 • — > or ■ p a — > 0. In the first region we write the following 
dipole formula 

m+l+a< a, 1, ....,771+ l||a, 1, ....,771+ 1 > m +l+a 

= ^ T> ijtk (p a ,p 1 , ...,p m+1 ) +V ijta (p a ,p 1 , ...,p m+1 ) + . . . (5.92) 

where . . . stands for terms that are not singular in the limit pi ■ Pj — > 0. 

The first term on the right-hand side of Eq. ( |5.92| ) is analogous to that in Eq. ( |5.1| ), i.e. 
the dipole V ijik (p a ,Pi, -,p m +i) is: 

Vij,k(Pa,Pl,";Pm+l) = ~7> ~ 

2pi ■ pj 

~ ~ T T ■ ■ ~ ~ 

m+a < a, 1, .., ij, .., k, .., m + 1| — 2 v V ij>h \a, 1, , .., k, .., m + 1 > m+a (5.93) 

ij 
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The second term is similar to the first one, that is, 



ay Ply Pm+lj 



2pi ■ pj 



T a -Ti- 

hg < a, 1, .., ij, .., m + 1[ a rra 3 V ijja \a, 1, .., ij, .., m + 1 > m+a , (5.94) 



rjyl 



but now the spectator is the identified particle a. It is thus convenient to rewrite Eqs. ( |5.3j ) 

(Pi + pj)p a 



( f5.9p in terms of the following variable 



PiPj + PjPa + PaPi 



(5.95) 



The dipole momenta p a and are: 

1 



1 — Z- ■ 

z Pa ) 



■Pa> Pij=Pi+P'j 



(5.96) 



and the splitting functions are given in terms of and the variables Sj, £,-, 



PiPa 



_ PiPa 

(Pj+Pi)Pa PijPa 

by the following expressions 

< S \ V q t 9,A^ 1 _ Z ij,a)W >= 87T/i 2e a S C F 



PjPa _ PjPa 



(Pj + Pi)Pa PijPa 



1 - Zi 



(5.97) 



< [i\V gi q a (zi)\v >= 8tt/2 e a s T i 



R 



1 Z{ Zij >a 

2 



[1 + Zi) - e(l - Zi 



PiPj 



(ZiPi ~ ZjP^j) (ZiPi ~ ZjPj) 



(5.98) 
(5.99) 



< V\Vg igj ,a(Zi'i 1 _ Z ij,a)W >= 167T/1 e « S CU 



1 Zij^ 



+ 



- 2 + (1 - e)— (Zipf - zjpfj) {zip\ - Zjp u j) 



1 Zj Zij^ a J PiPj 



(5.100) 



The only difference in the treatment of the two dipole contributions in Eq. ( |5.92|) is that 
a is identified in the final state and, hence, in the physical cross section we do not have 
to integrate over its momentum p a . Thus we are led to consider the following convolution 
formula for the phase space (the factor l/z 2 ~ 2e is introduced for later convenience: see 
Eq. Q5T7D and Sect, g) 



d<P{Pi,Pj;Q -Pa) 



1 dz 

- 1 - T d(f)(pij] Q - p a /z) [dpiipi^pa] z)\ 

0-2 



where 



d d p 



w 2-2e 



[dpi(pij,p a ; z)\ = ^y_ x S+iPi) 9(2)6(1 - z) 5(z - zy,«)- 



(5.101) 



(5.102) 
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In terms of the kinematic variables in Eqs. ( 5.95 , 5.97 ), Eq. Q5.102 ) is 

[dpi(pij,Pa,z)\ 



16tt 2 (2tt) ] 
(zi(l-5i))~ e 5(z-z iija ) (z(l-2:))~ £ , 

where dil^ d ~ 3 ^ is an element of solid angle perpendicular to pij and p a . 



(5.103) 



The integration of the splitting functions V ' over the phase space in Eq. ( |5.103| ) 



defines the functions Vij(z; e): 



2vrr(l 
that is, 



i /4vy 

-e) \2pijPa) 



Vij(z; e) = / [dpiipij^pa, z)\ - < Vy, (2j; 1 - z ij>a ) > , (5.104) 

^PiPj 



(1 



< VijJzi; I- z)> 



V ij (z;e)=e(z)Q(l-z) M - / ^ (5,(1 - 5 t )P ~ ' ~ ' • ( 5 - 105 ) 



57ras/i 



As in the case of the function Vij(x;e) in Eq. ( |5.50| ), Vij(z;e) is a z-distribution with 
e-dependent coefficients: 



Vij(z;e) 



V ij (z;e = 0)\ + + 5(l-z) j dz' Vy(/;e) + 0(e) . 



(5.106) 



The explicit form of the spin-averaged splitting functions < Vj Jia > is the same as in 
Eqs. ( |5.29| - |53T|) , apart from the replacements k — > a , 1 — y — > z: 



< V gg:a (z; l-z)> 

— = Uf 



2c 



1 — ZZ 



(l + 5)-e(l-5) 



(5.107) 



< V qq>a {z) > _ 



< V ggta (z; l-z)> 



2c 



2C, 



2c 



R 



\ 25(1-5) 
1 -e 



— + 



-2 + 5(1-5) 



(5.108) 



(5.109) 



1 — zz 1 — (1 — z)z 

Inserting Eqs. (|5.107|) -( |5.109| ) into Eq. ( |5.105|) , performing the 5-integration and expanding 
in e as in Eq. ( |5.106| ), we find 



V qg (z; e) = C F 



2 In- 1 



3 / z 



1-z 1 - zj+ 2 VI - z 



+ 5(l-z) V qg (e)+0(e) , (5.110) 



Z/ + 

11 / z 



V,(.-; e ) = 2C, yi _ z - l _ zJ+ 6 V1 _, 
where Vy(e) are the functions in Eqs. (|5 . 32j — |5734j) . 



5(1-0) V w -(e)+0(e) , (5.111) 
+ 6(1- z) V gg (e) + 0(e) , (5.112) 
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Let us now consider the singularities from the region where Pi ■ p a — > 0. In this case we 
write: 

m+i+a< a, 1, ....,m + l\\a, 1, ....,m + l > m+ i +a = ^ ^ai^iP^Pi, —,Pm+i) + ■■■ , (5.113) 



a; Pi) •••5 Pm+1, 



2>Pi ■ Pa 



'm+a 



< ai,l,...k,..,m + l\ — — 3 — V aijfc |ai, 1, ...fc, .., m + 1 > m+a . (5.114) 



The parton momenta in the dipole ( |5.114j ) are: 



Pa t 



z ik,a 



'Pa, Pk=Pi+Pk 



Pa 



where 



z ik,a 



{Pi + Pk)Pa 



PiPk + PkPa + PaPi 



(5.115) 



(5.116) 



Note that, unlike the case of the final-state dipole {ij, k) with no identified partons, here 
the momentum of ai (k) is (is not) parallel to a (k). 



The splitting functions V a i,k are: 

< s \Vqa9i,k( Z ik,a] Ui)W >= 87T/i 2e a S C F 5 S 



< s\V ga q i>k (z ik , a )\s' >= 8iru e a s C F 5, 



1 - Zik,a( l ~ U i) 

"1 + (1 - z lKa f 



- (1 + Zik,a) ~ e(l - Zik,a)) 

(5.117) 



z ik,a 



£ Zik,a 



(5.118) 



< A* I V qa qi,k{ z iha]Ui)\v >= 8iru 2e a s T R 

Ui jl-Uj) fpf_ 
PiPk \ Ui 



(5.119) 



^ifc,a)" 



A 1 
Pk 



1 - Ui 



Pk 



l-Ui 



1 



.1 - z ikA (l - Ui 



2 + 



z ik,a , 



+ (1 — e )(l ~~ z ik,a) z ik,a 



Ui(l - Ui) (rf f k \ (p\ p u k 



PiPk 



Ui 1 — Uj I V Ui 1 — Ui 



where the variable m is defined as follows 



PiPa 



(Pi + Pk)Pa 



(5.120) 



(5.121) 
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The phase-space for the dipole {ai, k} can be written in terms of the rescaled momenta 
PkiPai by using the following convolution formula 



d(f>(pi,p k ;Q -p a ) 



i dz 



z 



2-2e 



dcj){p k ;Q - p a i 'z) [dpi(p k ,p a ;z)] 



where 



[dpi(p k ,p a ;z)] = p n y-i - z ) K z ~ z ik,a)^ 



,2-2e 



U: 



(5.122) 



(5.123) 



Using the kinematic variables in Eqs. ( 5.1 16 J5T121|) , Eq. ( p. 123] ) can be written as 



[dpi(p k ,p a ; z)] 



( Ui (l - Ui )Y e 5{z - z iKa ) (z(l - z))~ e , (5.124) 



where dVt^ d 3 - ) is an element of solid angle perpendicular to p k and p a . 



The integration of the splitting functions V a i, k over the phase space in Eq. ( p. 12 1| ) 
\ 

as 1 / 47r/x 2 



defines the functions V a i, a (z; e ) : 



27rr(l - e) \2p k p aj 
that is, 



I ' [dpi(p k ,Pa,z)]—^— < V a i tk (z ikta ;ui) > , (5.125) 



V ai , a (z; e) = 6(^)0(1 - z) [z(l - z)r t — - «*))"' - V t k ^ > ■ ( 5 - 126 ) 

Jo Ui onasfi 



The spin-averaged splitting functions < V a i,k( z ; u ) > are: 

2 



< Vqa 9 i,k{z;u) > _ 

o — ^-f 



1 - z(l-u) 



(l + z )-e(l-z) 



(5.127) 



< V 9a mA z ) > 



2c 



c f 



l + (l-z) 



— ez 



(5.128) 



< V«„ft,k(*) > 



87ro;s/i 
< V ga9ijfc (z;M) > 



2c 



2z(l - z) 
1 -e 



87rag/i 



2c 



2C, 



1 7\ 7-2 + - + ^(l-^) 

1 — z[l — u) z 



(5.129) 



(5.130) 



Inserting Eqs. ( 5.127 - 5.130 ) into Eq. ( |5.126| ) and expanding in e, we obtain the following 
^-distributions 



V q , g (z; e) = -i P qg (z) + P qg (z) \n[z(l -z)}+C F z + 0(e) 



(5.131) 
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V g , q (z; e) = - X - P gq (z) + P gq (z) \n[z(l - z)\ + T R 2z(l - z) + 0(e) , (5.132) 



V M (*;e) = -~P qq {z) + 5{l-z) 



+ Cj 



4 In- 1 



1-z 1-z 



-5KV 

\nz 



1-z 



+ (l- z )-(l + z) \n[z{l-z)} 



(5.133) 



P 99 (^) + 6(1 -z) 



+ ™N f T R 



+ C A 



V gg (e) + N f V qq (e) 



-IT 



1 \ 2 , 

m ) + in z 



1-z 1-z 



+ 2 (-2 + z(l -z) + -) ln[*(l-z)] 



1 - z 



+ 0(e) 



(5.134) 



where the functions Vy(e) are given in Eqs. ( |5.32| - ^h3~4T ) and P a b(z) = P ab (z) are the regu- 
larized Altarelli-Parisi probabilities in Eqs. ( |5.85 - fT5£| ). 



5.5 Initial-state singularities: two initial-state partons 



Let us come back to the treatment of initial-state singularities in the case with two partons 
a and b in the initial state. In the limit p a ■ pi — > 0, the dipole factorization formula is 

m +i,ab< l,....,m+ l;a,b\\l,....,m+ l;a,b > m +i,ab 

= E^(Pl, -,Pm+l,Pa,Pb) + V ai '\ Px , ...,p m+ X]Va,Pb) + ■■■ ■ (5.135) 

This equation generalizes Eq. ( |5.60| ). 

The second term on the right-hand side represents a new dipole contribution in which 
the emitter is the initial-state parton ai and the spectator is the other initial-state parton b: 



V ahb (p h ...,p m+1 ;p a ,pb) 



1 



1 



^Pa ' Pi ^i,ab 



ab < 1, l;ai,b\ 



T b -T a 



|1, ...m + 1; ai, b > 



m,ab 



(5.136) 



This dipole contribution differs from that in which the spectator is a final-state parton 



because, when computing the cross section (see Sect. JT0|), it is convenient to leave the 
momentum p b unchanged. Thus the matrix element |1, ..,m + 1; ai, b > m ,ab involves an 
initial-state parton ai with momentum parallel to p a : 



Pai 



•Ei.ab 



•Ei,abPa > 

PaPb ~ PiPa ~ PiPb 
PaPb 



(5.137) 
(5.138) 
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and all other final-state momenta kj (and not only the momenta pj of the QCD partons!) 
rescaled as follows 



(K + K) 2 

where the momenta and are defined by 



K 2 



(5.139) 



Pa+Pb-P? , 
Pai+Pb ■ 



(5.140) 



Note that the momentum conservation constraint in the m + 1-parton matrix element 



is 



(5.141) 



Therefore we have 2 Y,j k, • K = 2K 2 and 2 % ■ (K + K) = 2K 2 + 2K ■ K = (K + K) 2 
and it is straightforward to check that momentum conservation is exactly implemented in 
the m-parton matrix element on the right-hand side of Eq. ( |5.136|) , that is, 



(5.142) 



Note also that Eq. ( EU39p can be rewritten in the following way: 



(K + K) 2 



K 2 



(5.143) 
(5.144) 



and thus the matrix K tl v {K,K') generates a Lorentz transformation (actually, a proper 
Lorentz transformation) on all the final-state momenta. 



The splitting functions V ai ' b in Eq. (gl35|) follows 

2 



< s\V^ 9i ' b (x i>ab )\s' >= 8n^a s C F 6 SS , 



1 — x 



i,ab 



- (1 + x i>ab ) - e(l - Xi, a6 )) 



(5.145) 



< s|V^^' b (x iiafe )| S ' >= 87r/i 2e a s T R [l-e- 2x i>a6 (l - x i>a6 )] 5, 



(5.146) 



Xi,ab 



1 - X i)afe 2p a • p 6 / Pip a \ f „ Pip a 

+ ! Pi Pb )\Pi Pb 

Xi,ab Pi -PaPi-Pb V PbPa J \ PbPa , 



(5.147) 



< fi\ V 9a9i '\x i>ab )\u >= 167r/i 2e a s C A 

1 - Xi,ab Pa ■ Pb 



-9^ f -, ^ \- Xi !ab {l - Xi jab ) 



1 - Xi 



ab 



, | , ' • i. 'in i n i'n I ii PiPa ti\ I u PiPa v 

+ I 1 - e) ! p£ ft P b 

^i.ab Pi 'PaPi-Pb V / V y 



(5.148) 
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The Jacobian factor associated with the Lorentz transformation Q5.143 ) acting on the 
final-state momenta is equal to unity. Therefore the phase space for the dipole {ai, b} has 
a trivial convolution structure: 



d<p(jpi,k 1 ,...;p a +p b )= / dx d(f)(k 1 ,...;xp a +p b ) [dpi(p a , p b , x)} 



where 



[dpi(p a ,p b ,x)) 



d d pi 



(27T) 



jS + (pi) 0(x)9(l - x) 5{x -X^ab) 



(5.149) 



(5.150) 



The phase space in Eq. ( |5.15Qj ) can be written as follows 

[dpi(p a ,Pb,x)} 



( 27i y-2e dVi dXi < ab "(^l 1 ~ X >) ®{ V i) e { 1 ~Y 



16tt 2 



(l - x) 2e 5(x - x iAb ) 



l — x 



1 — X 



Xy 

(5.151) 



where x ijab is defined in Eq. Q5.138 
perpendicular to p a and p b . 



PaPi/PaPb and dQ( d 3 -* is an element of solid angle 



The following integral of the splitting function V at ' defines the x-distribution y a ' a * ; 
2tt T(l - e) \2p a p b ) J 2p a pi n s {a) 

(5.152) 



Using the phase space in Eq. (|5.151|) , we thus have: 

ir 2 (i- e ) 



V a ' ai (x:e) = - 



eT(l-2e) 



e(x)e(i -x) (i-x) 



_ 2e n s (ai) < V ai ' b (x) > 



n s (a) 87rag/z 



2, 



(5.153) 



where the spin averages < V m,b (x) > are exactly proportional to the spin average of the 
ci-dimensional Altarelli-Parisi splitting functions in Eqs. fl4.18| - f4~21~D : 



nJai) < V ai ' b (x) > 



n s (a) 87ras/i 



2c 



=< P ~(x: e) > . 

a, era v ' ' 



(5.154) 



Performing the e-expansion in Eqs.( |5.l53"| , |5.154| ) according to the procedure in Eqs. ( |5.51| ) 
and ( |5.52| ), we find 



V a ' b (x;e) = V a ' b (x;e) + 8 ab T. 



a,b/ 



^ab rji2 

a 



2 In 1 



1 — x 1 — X 



+ ln(2 - x) 

1 — x 



+ K ab (x) + 0(e) 



(5.155) 

where V a ' b (x; e) are the functions defined in Eqs. (|5.8lH5T84l) and K ab (x) are given in terms 
of the regular part (see Eq. ( |5.89| )) of the Altarelli-Parisi probabilities as follows 



K ab {x) = P T ab g (x) ln(l-x) 



+ 5 ab T 2 



(5.156) 
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5.6 Factorization formulae with many identified partons 



In the most general case one deals with QCD cross sections involving initial-state partons 
and many identified partons in the final state. Here the NLO matrix element has collinear 
singularities when pi ■ pj — > (i and j being unidentified final-state partons) and when pi ■ 
Pa — > (a being either an initial-state parton or an identified final-state parton). The former 
singularities can be factorized in terms of the dipoles V^^, T>%, T^ij, a in Eqs. ( |5.2| . |5.86| . |5.94j ). 
As for the latter singularities, one should consider two different possibilities. If the spectator 
is an unidentified final-state parton k, one can factorize in terms of the dipoles T>'^,'D ai ^ 
in Eqs. (|5.61|j5TT4]) . If the spectator b is an initial-state parton or an identified final-state 



parton, it is convenient to introduce new dipoles DS.Dif.Pf'^pl")^ in which the 
momentum of the spectator is left unchanged. Actually, with respect to the momentum 
dependence, these objects are 'pseudo-dipoles' rather than dipoles, in the sense that they 
depend on the momentum p ai of the emitter, on the momentum p b of the spectator and on 
an additional momentum n. This momentum n is: 

n»=pl- £ Pa , (5-157) 

agfinal state 

where p^ n is the total incoming momentum and the second term on the right-hand side 
is the sum of all the momenta of the identified partons in the final statefj] Note that, 
by momentum conservation, n M is equal to the sum of the momenta of the final-state 
unidentified particles (QCD partons or not). Therefore the momentum n M is time-like (and 
with positive definite energy). Furthermore, since we only consider non-trivial quantities in 
which the lowest order process has at least one unidentified parton, it cannot be light-like. 

As we shall see, the dipole T> ai,b considered in the previous Subsection corresponds 
to X5(™) a *> 6 for the particular case with no identified partons in the final-state (i.e. when 

n = Pa+Pb)- 

Let us first consider the singularities pi ■ p a — > when a is an initial-state parton. In 
this case the factorization formula is: 

m+i...,a...< + 1, a..||l, ....,m + l,...;a.. > m +i...,a... 

= £Df ...,Pm+l, -;Pa, ••) +££f )( >l, :,Pm+l, -APa, ••) + ••• ■ (5-158) 

ky^i b^a 

The first term on the right-hand side of Eq. ( [5.158 ) is the same as that in Eq. ( 5.135 ). 



The new dipole contributions are given by the second term on the right-hand side of 
Eq. ( |5.158|) . Note that the sum over the spectators b refers to all initial-state partons as 



well as to all identified partons in the final state and no distinction is made between these 
two cases. The explicit expression for these dipole terms is: 

1 1 



u b U 7 ! ) ■■■iP 



'm+l> •••) Pa) 



2Pa ■ Pi X, 



arn 



(0 ...< l,..,m + l,...;ai...\ Th ^ m V { b n)ai |1, ..,m+ 1, ..; ai.. > m .., a .. ■ (5.159) 



■''Note that, in general, one can change p-^ by adding some momentum transfer that does not involve 
QCD partons. For instance, in deep-inelastic lepton-hadron collisions l(k) + h(p) — > l'(k') + . . . one can 
replace = + p^ with — + p^ where Qn = k„ — kl 
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The assignment of momenta in the matrix element |1, .., m + 1, ..; ai.. > m ..,a.. is the 
following. The momentum of the emitter is parallel to p a : 

K, = z«inK, (5-160) 
Pa • n 

and all other initial-state and identified final-state momenta are left unchanged. All other 
final-state momenta kj (QCD partons or not) are transformed according to the following 
Lorentz transformation 

kf = A^K, K) k] , (5.162) 

„ M . ,.-M±jgg& + *g. , ( , 16 a) 



(5.164) 



where the momenta and are defined by 

= n M _ (i _ Xa . n ) p M . 
Note that Eq. ( |5.162| ) actually defines a proper Lorentz transformation. 
The splitting functions V b ai in Eq. ( |5.158|) £1X6 cLS follows 
< s\V { b n)qa9l (x ain ; v iiab )\s' >= 8vr/i 2e a s C F 5 SS > [2v i>ab - (1 + x ain ) - e(l - x ain )\ , (5.165) 

< s\V { b n)9aqi {x ain )\s' >= 87r/i 2e a s T R [1 - e - 2x ain (l - x ain )\ 6 SS , , (5.166) 



< n\ Vt )qaq \x am )\u >= 87T/i 2e « s C F 



av 1 ^ n 

~9 -^ciin 



^Pi ■ Pa ( np a \ ( np a 

1 Pi ~ ^ Pi ~ n 



2(p a ■ n)(pi ■ n) - n 2 pi ■ p a \Pip a 1 I \PiP a 



, (5.167) 



< 



A*l vP^^ain^i^W >= 167r/i 2e o;s C A -g^ (^v i<ab - 1 + x ain (l - x ain )j 

,« v 1 - 2^ • p a / np a W np a \ 
f- (1 - e) — pt* ~ n ^ ~ n 

Xain 2(p a • n)(pi -n)-n 2 pi- p a \Pip a ) \PiPa J 



(5.168) 

where we have defined 

Vi,ab = i - : r • (5.169) 

Since all the momenta kj are obtained by means of the Lorentz transformation (|5.162p , 
the 'pseudo-dipole' phase space has a trivial convolution structure (note that in the physical 
region of interest, kj) 2 — Q 2 > 0), namely 

d(f>(pi, h, ...;p a + Q) = [ dxd<i>(ki,...]xpa + Q) [dpi{n = Q + p a ,p a ,x)] , (5.170) 

Jo 
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where 



d d pi 



[dpi(n,p a ,x)} = T^wzi 5+(Pi) 6(ac)6(l -s) 5(s -x ain ) 



(5.171) 



The phase space in Eq. (|5.171| ) can explicitly be written in terms of the kinematic variables 
x ain in Eq. (|5.161| ) and v { = p a Pi/p a n : 



[dpi(n,p a ,x)] = ^ P "™} 2 ^u_ 2e dl, i dx ™n 0(x(l - x)) 6(^) (l - — 



n 2 Vi 



16tt 2 (2vr) ] 



(1 - X) 26 8(X - X ain ) 



n 2 Vj 



1-x \ 2(1 — x)p a n / 



x)p a n J 
(5.172) 



where dVL^ d 3 - ) is an element of solid angle perpendicular to the light-like momenta 
n n _ n 2 p^ (2p a n) and p^. 



Performing the integration of the splitting kernels V^ a% over the phase space in 
Eq. ( |5.172|) , we introduce the functions V a,ai (x;e;p a ,p b ,n) : 



as 1 / 47ryU 2 



2vrr(l-e) \2p aPb/ 



V a ' ai (x;e;p a ,p b ,n) = / [dpi(n,p a ,x)] 



1 n s (ai) 
2>PaPi n s (a) 



v (n)ai ( . s 
^ r b y^aim u i,abj " 



(5.173) 



The spin averages of the kernels V b ai are related in a simple way to the corresponding 
averages of the (ci-dimensional) Altarelli-Parisi splitting functions, that is, 



n s (ai) < V { h n)a \x]v^ ab ) > 



n,[a 



57ras/i 



2r 



i 



VLab 



1 — X 



(5.174) 



Therefore, the only non-trivial integration in Eq. (|5.173|) is that which involves the term 
Vi,ab (see Appendix B) and leads to the following result 



V a > b (x;e;p a , Pb ,n) = V a ' b (x; e) + £ a >\x;p a ,p b ,n) + 0(e) . 



ja.b 



a,b I 



(5.175) 



The first term on the right-hand side is given in Eq. (|5.155|) and the second term is defined 
by 



C a \x;p a , Pb ,n) = 5 ab 5{l-x)2T 2 a 



ab 



Li 2 1 



[1+v) (Pa+Pb) -n 



+ Li 2 1 - 



[1 - v) (Pa+Pb) ■ n 
2 Pa ■ n 



Pa-n J 

_ pab (r)]n n 2 (Pa-Pb) 

regiJ 2{p a -nf 



(5.176) 



_ n 2 (p a +p b ) 2 
\ [(Pa+Pb)-n}< 



(5.177) 



P?eg( x ) being the regular part of the Altarelli-Parisi probabilities in Eq. ( |5.89| ) and Li 2 ( 
is the dilogarithm function: 



x 



Li 2 (x) 



c dz 
n z 
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ln(l - z) . 



(5.178) 



Note that unlike the V-functions considered in the previous Subsections, Eq. ( |5.175j) 
depends not only on the momentum fraction x, but also on the momenta p a ,pb,n. This 
momentum dependence is entirely accounted for by the C a,b function on the right-hand side 
of Eq. ( |5.175| ). In the case with no final-state identified partons, p b is the momentum of an 
incoming parton, n = p a +Pb, so C a,b vanishes, thus recovering the results already discussed 
in Sect. |5~5|. 



In order to deal with the singularities for pi ■ p a —>■ 0, when a is an identified parton in 
the final state, we introduce the following factorization formula 

m+la...... 

< 1, ....,m + l, a...; ...||1, ....,m+ I, a..; .. > m+la ...... 

= Y, V aiAPU -,Pm+l,Pa, ••; •••) + -,Pm+l,Pa, ■■] •••) + • • • , (5.179) 

where V ai ^ is the dipole in Eq. (|5.114|) and the new dipole contribution is given by 



2p a ■ Pi 



■ ma..,...< 1, ..,m + l,ai..; ...\ b 2 V^lT, ..,m + I, at..; .. > mo ...... . (5.180) 

ai 

As in Eq. ( |5.158j ), the sum over the spectators b in Eq. ( |5.179| ) refers to all initial-state 
partons as well as to all identified partons in the final state with no distinction between 
these two cases. 

The assignment of momenta in the matrix element |1, .., m + 1, ai..; .. > ma ...... is the 

following. The momentum of the emitter is parallel to p a : 

Pi = —P"a, (5-181) 



"am 



Pa ■ n 



(Pa +Pi)-n 



(5.182) 



and all other initial-state and identified final-state momenta are left unchanged. All other 
final-state momenta kj (QCD partons or not) are transformed according to the Lorentz 
transformation in Eqs. ( |5.162| - |5".164j) . 

The splitting functions in Eq. ( |5.179| ) are given in terms of the variable v^ a b in 

Eq. ( |5TT69D as follows 



< S \ V qa gi ,b( Z ain] V itab )\s' >= 87Vj2 2t a S C F 5 S 



2^-(l + z .)- e (l 



, (5.183) 



< s\V 



(n) ( 
g a qi,b\ 



\s' >= 8nfi 2e as Cf 



l + (l-z a 



5ss> , 



(5.184) 



< V\ V tli,b( Z ain)W >= 87T/i 2e a S T R 



_ /^ _ A z . (I - z 



Pi 'Pa 



2(p a ■ n)(pi -n)-n 2 p t ■ p a \Pip a 



n Pa p a 

Pi — n 



nPa v 
Pi - n 

,PiPa 



(5.185) 
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Pi 'Pa 



/J.U 



VLab 



- 1 + 



1 — Z„ 



+ 2(1 - e)z ain {l - z a 



np a 



2 (Pa • n)(pi -n) -n 2 pi- p a \pip a 



Pi - n 



A i nPa ^ 



K PiPa 



Pi -n 



(5.186) 



The phase space for the dipole (|5.180 ) has the following convolution form 

- 1 dz 



fci, Q-Pa) 



where 



o z 



d d p 



2-2e 



d(j)(ki,...;Q -p a /z) [dpi(n = Q - p a ;p a ,z)] , (5.187) 



[dpi(n;p a ,z)\ = / 27r w_! h(Pi) &(z)Q(l - z) 5(z - z ain ) z 



2-2e 



(5.188) 



Introducing the kinematic variables z a i n in Eq. ( |5.182| ) and Vi = PaPi/p a n } Eq. ( |5.188 
be written as follows 



can 



n 2 VjZ 



[dpi(n;p a , z)) = ° a ^ 2 f ^y_ 2e dvi dz ain Q(z(l - z)) Q(vi) 0l 1 - ^TfZ^ 



16tt 2 (2tt) ] 



n 



1 — z) 2e 5(z — z a 



VjZ 

1 - z 



1 - 



n 2 VjZ 



2(1 - z)p a n 



(5.189) 



where dVL^ d 3 - ) is an element of solid angle perpendicular to the light-like momenta n M 
n ii _ n 2 p^ (2p a n) and p%. 



Thus, the integration of the splitting kernels V a ™\ defines the functions V a i ja (z; e; p a , p b} n): 



27rr(l-e) \ 2p a p h/ 



Vai, a (z;e;p a ,p b ,n) = / [dpi(n;p a , z)) 



1 



ZPaPi 



Vai!bi Z ain'i v i,ab) > 



(5.190) 

As in the case of initial-state singularities, the spin averages < V„™ \ > are simply related 
to the (d- dimensional) Altarelli-Parisi splitting functions: 



<V { £ b (z;v hab )> 



2c 



< P~ (z: e) > +2 5 ai ' a T 2 



VLab 



Z 1 — Z 



(5.191) 



and, explicitly performing the integration of the Vi >ab contributions in Eq. ( |5. 191 ) (see 
Appendix B), we find that the functions V a ,b are equal to the analogous functions V a,b in 
Eq. Q5TT751) : _ 

Va :b (z;e;p a ,p b ,n) = V a ' b (z;e;p a ,p b ,n) . (5.192) 
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6 QCD cross sections at NLO 



In the following Sections we describe in detail our subtraction method for evaluating QCD 
cross sections. To this end, it is useful to recall the general and precise definitions of the 
NLO cross sections. 

In the case of processes with no initial-state hadrons, for instance in e + e~ annihilation, 
the partonic cross section is| 

a = a LO + a NLO , (6.1) 
a LO = f da B , a NLO = f da R + f da v , (6.2) 

Jm J m+1 Jm 

where da B , da R , da v are the cross sections in the Born approximation and to one-loop order 
(R: real emission; V: virtual correction). If a is a jet cross section (no final-state hadrons 
observed), hadron- level and parton- level cross sections coincide. 

In the case of processes with one initial-state hadron carrying momentum p M (for in- 
stance, deep inelastic lepton-hadron scattering), the calculation of the QCD cross section 
requires the introduction of parton distributions. If we denote by faijIi^V) the density of 
partons of type a in the incoming hadron, the hadronic cross section is given by 

°{p) = E f ^ fa(v, /4) [ai°(vp) + of^fap; /4)1 , (6-3) 
and the corresponding parton-level cross sections are: 

* L a°(p) = I da B (p) , (6.4) 

Jm 

°a L °(p;ti) = I da R (p)+ f daV(p)+ f da^^). (6.5) 

Jm+l Jm Jm 



The notation B, R, V is as in Eq. fl6.2p . The contribution da^ represents the collinear- 
subtraction counterterm and is explicitly given by the following expression 



e \ PF / 



da^(zp) . (6.6) 



The partonic contributions on the right-hand side of Eq. (|6.5| ) are separately divergent 
for e — > 0. Their sum a^ LO is finite for e — >• but depends on the factorization scale and 
on the factorization scheme of collinear singularities. Both dependences are contained in 
the definition of da^: fip is the factorization scale and the actual form of the kernel K } 



ab 



specifies the factorization scheme. Setting K£ b & (z) = defines the MS subtraction scheme. 
The functions P ab (z) in Eq. ( |6.6|) are the four dimensional Altarelli Parisi probabilities in 
Eqs. ( |5.85[ - |5TH^) . The parton densities f a (r],fip) are also scale/scheme dependent, so that 
this dependence cancels in the hadronic cross section of Eq. (|6T 



Note that 

W dxx P a \x) = 1 , (6.7) 



§We are using the same notation as in Sect. Thus, m is the number of unobserved final-state partons 
for the leading-order process. 
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and that, in the MS scheme, momentum conservation reads as follows 



V f dxX f a {x,fip) = 1 . 

In other factorization schemes the generalization of Eq. (|6.8|) is: 



(6.S 



V / dx x f b (x,fi F ) 
77 Jo 



cba «S 



Z7T JO 



(6.9) 



The extension of Eq. ( |6.3|) to processes with two initial-state hadrons is straightforward. 
Denoting by f a and f b the parton densities of the two incoming hadrons, we have 



CT (P> P) = / d V fa(V, /4) / dV fb(V, t4) \ a ab°(VP, VP) + °ab L °(VP, VP', /4 



u,6 



(6.10) 



4?(P,P) 



d °ab(P,P) 



^ab L °(p,P;^ 2 F ) = / da a R b (p,p) + da^ b (p,p)+ da^ b (p, p; /4) 
where the collinear counterterm is: 

■ U d 6(l - z) 
+ 8 ac S(l-z) 



(6.11) 
(6.12) 



cd 

-U¥)' p " w+asw 



(6.13) 



It is completely trivial to generalize the resulting formulae to the case in which one intro- 
duces different factorization scales for the two hadrons, as one might in photoproduction 
for example. The replacement f b (fj,fi F ) — > fb{fj,pp) m Eq. (|6.10| ) is simply accompanied 
by fip — ► ftp in the second term in the curly bracket on the right-hand side of Eq. ( |6.13| ) . 

Let us now consider fragmentation processes. The one-hadron inclusive cross section 
C(ind) (p) in the case with no initial-state hadrons is: 



O-(md) {P)=J2 — d a (V, /4) \rfincl)a(p/ V) + rfincl) a(p/m /4) 



(6.14) 



where d a (r], up) is the fragmentation function of the parton a into the observed hadron and 
the partonic cross sections are: 

°(iZi)a(p) = / daf ind)a (p) , (6.15) 
^ind)a(p) = J dcrf ind)a (p) + J m da( incl)a {p) + J daf- nd)a (p; fi 2 F ) , (6.16) 
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C I 2 \ a s 1 i-^ f 1 dz 



1 /47T/X 



2 N 



/4 



P 6a (z) + tf£f(z) 



d°(incl)b(p/ z ) 



(6-17) 

where the Altarelli-Parisi probability P a b(^) describes the time-like splitting process a(p) — > 
b(zp). Note that P ab {z) = P ab (z). The kernel H^(z) in Eq. . 1 7[ ) defines the factorization 
scheme (H^(z) = in the MS subtraction scheme). 



Note that the one-particle inclusive cross section (J(i nc i) (p) in Eq. ( |6.14[ ) is defined with- 



out integrating over any component of the momentum p^ of the observed hadron (the 
unusual convolution measures drj/rj 2 and dz/z 2 ~ 2e in Eqs. ( |6.14| ) and ( |6.17| ) follow from 
that). Thus the following integral 

d A P r /_2 



(2tt)* 



5 + (p ) a {ind) (p) (6.18) 



is equal to the associated multiplicity times the total cross section. The corresponding 
associated multiplicities at partonic level (i.e. in d dimensions) are (/ = B, R, V,C,): 

Jr^y-l S+(P 2 ) d(j\ ind)a (p) (6.19) 

In the most general case, one should consider multi-particle correlations, that is, one 
deals with processes of the type 

p + p-^ qi + ...q n + X , (6.20) 

where p, p are the momenta of two incoming hadrons, qi,...q n are the momenta of n 
hadrons detected in the final state and X stands for unobserved final-state particles or jets. 
Note that, by definition, the momenta qi,...,q n are supposed not to be parallel to each 
other or to the incoming momenta p and p. The hadronic cross section is: 

V(incl)(p,P;qi,~,qn) = J2 I dr lfa(V,^F)[ dfj f b (fj,fJ. 2 F ) 

a,b J ° J ° 

a^aJO Vl JO T)l 

I ® ab, (incl)ai ,..,a n {vp,vp; qi/vi,-,qn/vn) 

+ ^abf(?nci) ai ,.., an ivp,vP;qi/vi,---, qjvn] i4>, t4, ■■, vl)} , (6-21) 

where we have introduced a single factorization scale n 2 F for the parton densities and n 
different factorization scales fi 2 , for the fragmentation functions. 

The corresponding cross sections at parton level are the following 



^ab,(incl)au..,a n (P'P'^^-'1n) = dff^f^^ (p, p]q U ...,q n ) , (6.22) 

Jm 

a ab L (incl)a 1 ,..,aS'P->P'i ?1> ■"' 9n', l^Fi A*!' "> Mn) = / ^Siinc/W ,..,a„ (P> P\ <?1 j • • • j Qn) (6.23) 

Jm+1 

+ / da^ b ^ incl ' )ai ^ n (p,p;qi,...,q n )+ da ab ( incl \ ai>an (p,p-, q±, q n ; fx F , .., /4) , 

Jm Jm 
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where the collinear counterterm da is given by 



da c 



/ - .22 2\ _ 1 

ab,(incl) ai,..,a n 

{p,p; qx, ...,q n ; /ip,/^, --^J = 

1 /47T/2 2 



+ E / dz 



e V up 



P aa (z) + K£ (z) 



e \ H'F J 



da a'b,{incl) ai,..,a n { Z Pi P\ QU ••) 9n) 
d^ab' ,{incl)a 1 ,..,a n {P-' Z P'-> ?1> — > 



i=l J0 ^ 



2-2e 



d(T ab,{incl) a 1 ,..,a' i ,..,a n (P^ P'l Ql) "J Qi/ ' Z i "i In) 



(6.24) 



Note that, computing the hadronic cross sections in Eqs. ( |6 . 3| , |6 . 1 0| , |6 . 1 4| , |6T2~1"| ) , the factorization- 
scale evolution of the parton distribution functions has to be consistently carried out at 
NLO. For the parton densities f a (rj, fip) and the fragmentation functions d a (rj, n 2 F ) we have 



din fi 2 F 



v- i l dz f i i 2n as 04) 



Tj Z 



2tt 



p 6 a (z) + Os04) pi l )ba , z) 



(6.25) 



dd a (rj,fj, 2 F ) 
din ftp 



E 



1 dz as(np) 



n Z 27T 



Pab(z) + 



2tt 



d b (r)/z,n F ) 



(6.26) 



where P^ ab (x) and P^{x) respectively are the space-like and time-like NLO Altarelli- 
Parisi probabilities. They depend on the factorization scheme and this dependence is given 
in terms of the flavour kernels Kf^ and Hlf as follows 



1 dz 



P {1)a \x) = P£! ab (x)+i: / -\P ac (z)K*(x/z)-K™(x/z)P cb (zj\ -2irP K^(x 



(1) MS 
ab 1 



X 



x Z 



1 dz 



(6.27) 



E / - [Hlt(xfz) P cb (z) - P ac {z) Hl£{x/z)\ - 2kPoB£{x) 



x Z 



_ (6.28) 
where P^ a0 (x) and P^ )MS (x) are the corresponding probabilities evaluated in the MS 



subtraction scheme 26 
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7 Jet cross sections with no initial-state hadrons 



In processes with no initial-state hadrons, the QCD cross section for jet observables is given 
by Eqs. (|6.1|J6.2| ). In terms of the QCD matrix elements, the Born-level cross section in d 
dimensions is the following 

da B =J\f in ^2 d(j) m {p u ...,p m ;Q) — !— \M m (pi, ■■.,p m )\ 2 Fj m) (pi, ...,p m ) , (7.1) 
{m} b i m } 

where J\f in includes all the factors that are QCD independent, J2{ m } denotes the sum over 
all the configurations with m partons, d(f> m is the partonic phase space in Eq. (|3.1|), S{ m y is 
the Bose symmetry factor for identical partons in the final state and Ai m is the tree-level 
matrix element. 

The phase space function Fj(pi, ...,p m ) defines the jet observable in terms of the mo- 
menta of the m final-state partons. In general Fj may contain ^-functions (thus, Eq. ( |7.1| ) 
defines precisely a cross section), ^-functions (Eq. ( |7.1| ) defines a differential cross section), 
numerical and kinematic factors (Eq. ( |7.1| ) refers to an inclusive observable), or any combi- 
nation of these. The essential property of Fj is that the jet observable we are interested 
in has to be infrared and collinear safe. From a formal viewpoint this implies that Fj fulfils 
the following properties 

F^ n+1 \ Pl ,.., Pj = \q,..,p n+ i) -+ Fj n) ( Pl ,...,p n+1 ) if , (7.2) 

Fj n+1) (pi,~,Pi,:,Pj,..,p n +i) ^ FP(p h ..,p, ..,p n +i) if Pi^zp, pj->(l-z)p (7.3) 
for all n > m, and 

Fj m) ( Pl ,...,p m )^0 if Pi-pj^O . (7.4) 



Equations (|7.2|) and ( |7.3| ) respectively guarantee that the jet observable is infrared and 
collinear safe for any number n of final-state partons, i.e. to any order in QCD perturbation 
theory. The n-parton jet function Fj 1 ^ on the right-hand side of Eq. ( [7.2|) is obtained from 
the original Fj ra+1 - ) by removing the soft parton pj, and that on the right-hand side of 
Eq. ( |7.3|) by replacing the collinear partons {pi,pj} by p = Pi + pj. 

Equation (|7.4j ) defines the leading-order cross section, that is, it ensures that the Born- 
level cross section da B in Eq. ( |7. 1| ) is well-defined (i.e. finite after integration) in d = 4 
dimensions. 

The cross section da R has the same expression as da B , apart from the replacement 
m — > m + 1. 



7.1 Subtraction term 

In order to compute a NLO we write the following identity 



a NLO 



I (da R - da A ) + \[ da A + f da v 



(7.5) 
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where da A is a local counterterm for da R , i.e. da A has the same (unintegrated) singular be- 
haviour as da R . An explicit and general form of do A is provided by the dipole factorization 
formula introduced in Sect. |B"7[| Thus we can define: 

do A = N in -,Pm+l,Q) ^ (7-6) 

{m+1} ^{m+1} 

' S V idAV\i -,Pm+l) Fj (Pi, ~Pij,Pk, -,Pm+l) , 

pairs fe^jj 

U 

where V^ipi, ■■■,p m +i) is the dipole contribution in Eq. (|5.2D , namely 

1 

?>ij,k {Pu -,Pm+l) = —z (7-7) 

2pi ■ pj 

m< 1) ■•; ijt ••) k, .., m + 1| — 2 - Vij : k |1, .., ij, .., fc, .., m + 1 > m , 

ij 

and Fj m \pi, ..pij,p k , ..,p m+ i) is the jet function for the corresponding m-parton state 

{Pi, -Pij,Pk, --jPrn+l}- 

We can check that the definition (|7.6|) makes the difference (da R — da A ) integrable in 
d — 4 dimensions. Indeed, its explicit expression is 

da R -da A = N in ^ d(/) m+1 (pi, ...,p m+1 ;Q) — — 

{■m+l} ^{m+l} 

■ {\M m+1 ( Pu -,Pm + i)\ 2 F!r +1 \pu- 

~ H H V ij:k (p u ...,p m+1 ) F!) m \p 1} ..p ij} p k} ..,p m+1 )} . (7.8) 

pairs k^l,j 

Each term in the curly bracket is separately singular in the soft and collinear regions. 
However, as discussed in Sect. |5J], in each of these regions both the matrix element M. m +\ 
and the phase space for the m + 1-parton configuration behave as the corresponding dipole 
contribution and dipole phase space: 

\M m+1 {pi, ...,p m+1 )\ 2 V ij;k (px, ...,p m +i) , (7.9) 
{p 1 ,..pi,..p j ,..p k ,..,p m+1 } {p u ..p ij} p k} ..,p m+1 } . (7.10) 

Thus, because of Eqs. ( |7T2|) and ( |7.3|) , the singularities of the first term in the curly bracket 
are cancelled by similar singularities due to the second term. On the other hand, each 
dipole T>ij k in Eq. ( [7.7|) has no other singularities but those due to the m-parton matrix 
element |1, .., k, ..,m + 1 > m . However, because of Eq. ( [7.4|) , these singularities are 
screened (regularized) by the jet function Fj(pi, ..pij,p kj ..,p m+ i) in the curly bracket of 
Eq. 



Note that this cancellation mechanism is completely independent of the actual form of 
the jet defining function but it is essential that da R and da A are proportional to Fj m+1 ^ 
and F^ respectively. Nonetheless, both da R and da A live on the same m + 1-parton phase 
space d<p m+1 . Thus the numerical integration (in d = 4 dimensions) of Eq. ( |7.8| ) via Monte 
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Carlo techniques is straightforward. One simply generates an m + 1-parton configuration 
and gives it a positive (+ |.M m+ i| 2 ) or negative (— J2k^i,j ^y,fc) weight. The role of the 
two different jet functions _pj m+1 ) and Fj 71 " 1 is that of binning these weighted events into 
different bins of the jet observable. Any time that the generated configuration approaches 
a singular region, these two bins coincide and the cancellation of the large positive and 
negative weights takes place. 



7.2 Integral of the subtraction term 



Having discussed the four-dimensional integrability of (da R — da A ), the only other step we 
have to consider is the <i-dimensional analytical integrability of da A over the one-parton 
subspace leading to soft and collinear divergences. 



We start by noting that the dipole contribution in Eq. (|7.7|) can be written as follows 



T^ij,kiPli ■■■iPm+l) 



0n tJ '* 7^2 \M% k (pi, ..Pij,Pk, -,Pm+l] 
z Pi " Pj -L ij 



(7.11) 



J h 



where M. m ,k is a colour-correlated m-parton amplitude (see Eq. (|3.9|)) depending only on 
Pi, -Pij,Pk, -,Pm+i while V ijtk /pi ■ pj depends only on Pi,Pj,p k or, equivalently, Pi,Pij,pk 
(the subscript h in the square bracket of Eq. (|7.11|) means that Vy^ and M. m ^ are still 
coupled in helicity space). Using the phase space factorization property in Eq. ( p,17| ), we 
can thus completely factorize the Pi dependence in the following form 



da J 



m+l 



f Mn ^2 d</>m(PU~Pij,Pk,;Pmrhl]Q) -q— 

Jm {m+l} P airs k^i,j ^{m+1} 



Fj m \pi } ..Pij,Pk,-;Pm+l) 
7fW \M% k (pi, ..pij,p k , ...,Pm+l) 



(7.12) 



V 



ij,k 



1 2 Pi-Pj\ h 



[dPi(pij,Pk)} 



and we can perform the integration over p^. In particular, according to the discussion in 
Sect. |5J], the azimuthal correlation between Vy^ and |.M^ ,A; | 2 vanishes after integration 
over [dpi(pij , pk)] and we get (see Eq. (|5.27|) ) 

/ da A = - f N in XI Y f d(l) m (p 1 ,..p ij ,p k ,.,p m+ i;Q)-— — 

Jm+1 Jm {m+l} Pairs k^i,j ^{m+1} 



• Fj^iPl, -Pij,Pk, -,Pm+l) |A4^(pi, ..Pij,p k , ...,p m+1 )\' 



as 1 ( 4ttm 2 y 1 
27rr(l-6) ^2p«pJ T 2 ^' lej - 



(7.13) 



In order to rewrite Eq. ([7.13 ) in terms of an m-parton contribution times a factor, we 
have to perform the counting of the symmetry factors for going from m partons to m + l 
partons. 
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Consider a Born-level m-parton configuration with rrif quarks of flavour /, mj anti- 
quarks of flavour /, and m g gluons. From this parton configuration we can obtain an 
(m + l)-parton configuration by changing 



a m n — > m. 



+ 1 or b) rrif — > rrif + 1, rrif — > rrif + 1, m g — > m 9 — 1 . (7.14) 



The corresponding ratios of the symmetry factors for identical particles are 

1 



q(o) 



S 



{m+l} 
S{m+1} 

Thus we write 

1 



[nig + iy. rrig + 1 
.... nif \ mj! rrigl 



(7.15) 



m fi 



(m/ + 1)! (m f + 1)! (m 5 - 1)! (m f + l){m f + 1) 



{m+l} p. 



airs 
j 



E 

m+] 

+ E 



(a) 1 



{m+ l} ^m} m 9 + 1 
(b) 1 



pairs 



and then 



E ••■+ E •••+ E 

pairs pairs 

i=qf.g U=9f, 

{^1} 5 W ( m / + + X ) , ^ 

#(«, j)m+l 



pairs 



\ 
) 



E 

pairs 



#(y)r 



E 



(7.16) 



(7.17) 



where j) m +i denotes the number of pairs i, j in the configuration with m + l partons 
and #(ij) m denotes the number of partons with flavour ij in the corresponding m-parton 
configuration. Since we have 



#(g/»fl)m+i r m f (m g + 1) #(g/,ff) m +i = m/(m g + 1) 
#(#,#)m+i m g (m 9 + l)/2 #(g/, g~/) m +i _ (m/ + l)(m/ + 1) 



(7.18) 



we end up with 



m n 



1 



(a) 1 



E o^— £••• = £ w — I £•••+£ 

{ m +l} °{m+l} pa 



airs 



{m} 



2 i 1 - 



K ij=q f i]=q f 



V=9 



{m} °W 

Inserting Eq. ( |7.19| ) into Eq. ( |7.13| ), we obtain: 

/ da A = -f N in ^d<j) m {pi,..,p m ,Q) — — Fj m) (p!,...,p„ 

Jm+l Jm r i u/m.) 



(7.19) 



{m} 



S{ m } 
47T// 2 



EE^W-,^, 2 , r(1 _ e)V2pi . Pfc ; T 



7^V 4 (e) , (7.20) 
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where we have defined 



V*(c) 
VAe) 



if i = q,q 



l^ 9 M+N f V m {e), if i = g . 



(7.21) 
(7.22) 



Equation ( |7.2Q ) explicitly shows that the subtraction contribution da A can be integrated 
in closed analytic form over the subspace leading to soft and collinear divergences. These 
divergences are indeed collected in terms of e poles into the factors V;(e). 



The final result for J m+1 da A in Eq. ( fT^DQ can be written as follows 



da' 



m+l 



da B ■ /(e) 



da B ■ 1(e) 



(7.23) 
on the right- 



Comparing Eq. ( [7.201 ) with Eq. (|7.1|), we see that the notation 

hand side means that one has to write down the expression for da B and then replace the 
corresponding matrix element squared at the Born level 



\M 7 



m < 1, — , J7i|l, — ,fn > m , 



by 



m < 1, ....,m| /(e) |1, ....,m > m , 



(7.24) 



(7.25) 



where the insertion operator /(e) depends on the colour charges and momenta of the m 
final-state partons in da B : 



I(px, ...,p m ;e) 



27rr(l-e)^ T\ A > \2 Pi .p K 



(7.26) 



The singular factors Vj(e), defined in Eqs. ( [7.2 H , |7.22|) , are given by (see Eqs. (|5.32| - |5~34]) ) 



where ji is defined in Eq. ( ^.90|) and we have introduced the following constants 

(1 ^ 



i=q,q 



I C TP , Kj — n 

2 6 9 



18 6 / 9 1 



(7.27) 



(7.28) 



Note that all the terms in Vj(e) have a simple interpretation. The coefficient of the 
double pole 1/e 2 is the square of the colour charge of the parton i, that of the single pole 
1/e is related to the integral of the non-soft part of its four-dimensional Altarelli-Parisi 
splitting function, and the 7r 2 -term is a customary phase space factor. The constant K g 
typically appears in the resummation program of higher-order logarithmic corrections of 
Sudakov type [2~7| if one uses dimensional regularization and the MS renormalization scheme 
(we do not know of an analogous role for the constant K q ). Actually, in the context of 
these calculations the constant K g can be absorbed by a redefinition of the renormalized 
coupling, thus introducing a 'more physical' renormalization scheme, accidentally called 
the 'Monte Carlo' scheme in Ref. |28f . 
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7.3 One-loop corrections 



The NLO QCD cross section in Eq. ( |7.5|) is finite by definition, i.e. because of the infrared 
and collinear safety of the jet observable. Since the first term on the right-hand side is 
finite, the second term 



m+l 



do A + [ da v = f {da v + \da B ■ 1(e)] } 

.I'm .I'm L -J 



(7.29) 



has to be finite as well. Thus all the e poles in da v must be cancelled by those in da B ■ 1(e). 
The virtual contribution da v has the following expression in terms of the (renormalized) 
one-loop matrix element 



da V = Mn E #m(Pl, -,Pm, Q) c 

{m} b i m } 



\M m (Pl, -,Pm)\(i —loop) J 

(pi,...,p m ) ■ (7.30) 



Comparing Eq. (|7.30|) with Eqs. ( |7.20|]7T2"3"D , we find that the singular terms of the one-loop 
matrix element have the following universal structure 



\M m (pi, -,Pm)\li-ioop) = ~ m< 1, ....,m| 1(e) |1, ....,m > m + 

' Att/j 2 



2^r(l-e)? T? Vl{e) S \2 Pl - Pk 



\M l m h (Pi,-,Pm)\ 2 + ... , (7.31) 



where |.M„ fc | 2 = m < l,....,m|Tj • T k \l,....,m > m is the square of the colour-correlated 
tree-amplitude in Eq. ( |3.9|) and the dots stand for contributions that do not contain any 
e poles. Thus, using the finiteness property of the NLO cross section, we have obtained 
as by-product of our algorithm the general expression (|7.31| ) for the singular terms of the 
one-loop QCD amplitudes. 



Alternatively, we can use the results of Ref. E| (see also Ref. [|]) on the singular 
behaviour of the one-loop amplitudes to prove that our algorithm correctly produces the 
cancellation of all the e-poles in Eq. ( 7.29Q , thus leading to a finite NLO cross section. As 
a matter of fact, using Eqs. (|7.27|) and keeping only the e poles, we can rewrite Eq. Q7.31]) 
as follows 



\M m (pi, --,Pi 



m) | (i— loop) 



2nT(l 



EE 



Att/i 2 
2pi ■ Pk 



8 < 1, ....,m| Ti -T k \l, ....,m > m + 



27rr(l-e) 
1 



E 



i 



47T/i 2 



i,k 



E7i 



v 2pi -p k 
< 1, ....,m|l, 



,< 1 ml T, T k II, 



1 7i_ 
e T? 



m >, 



m > m + 



(7.32) 



where in the last expression we have used colour-charge conservation (J2k^iTi-T k = —T 2 ). 
Equation (|7.32|) is completely equivalent to Eqs. (2.3) and (2.9) in Ref. p9|Pl 

^Note a misprint in Eq. (2.9) of Ref. j29): the overall sign on the right-hand side should be plus instead 
of minus. 
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7.4 Final formulae 



The final results of the application of our algorithm to the calculation of jet cross sections 
with no hadrons in the initial state are summarized below. 

The full QCD cross section in Eq. ( |6.1|) contains a LO and a NLO component. Assuming 
that the LO calculation involves m final-state partons, the LO cross section is given by 



a 



LO 



[ da B = [dQ™ \M m (pi,..., Pm )\ 2 F^\ Pl ,...,p m ) , (7.33) 



where Ai m is the tree-level QCD matrix element to produce m partons in the final state and 
the function F^ defines the particular jet observable we are interested in (see Eqs. (|7.2| - 
7~% for the general properties that Fj 71 ^ has to fulfil). The factor d& m ^ collects all the 
relevant phase space factors, i.e. all the remaining terms on the right-hand side of Eq. ( |7. 1| ) . 
The whole calculation (phase space integration and evaluation of the matrix element) has 
to be carried out in four space-time dimensions. 



According to the notation in Eq. fl2.13|) , the NLO cross section is split into two terms 
with m + 1-parton and m-parton kinematics, respectively: 



a 



NLO 



a NLO{m+l} _|_ a NLO{m} 



(7.34) 



The contribution with m + 1-parton kinematics is the following 



a 



NLO{m+l} 




dipole 



KLo- 

j rf$( m+1 ){ \M m+1 ( Pl ,..., Pm+1 )\ 2 F^ +1 \ Pl ,..., Pm+1 ) 

-Y, v ljik ( Pl ,..., Pm+1 ) F^ m) ( Pl ,.. Pij , Pk ,.., Pm+1 ) | 

pairs ky^i,j 



(7.35) 



where the term in the curly bracket is exactly the same as that in Eq. ( |7.8| ): M m+ i is the 
tree-level matrix element, T>ijh is the dipole factor in Eq. (|5.2|) and F^ is the jet defining 



function for the corresponding m-parton state (note, again, the difference between the two 
jet functions Fj m+1 ^ and Fj m ^ in the curly bracket). Despite their original (^-dimensional 
definition, at this stage the full calculation is carried out in four dimensions. 

The NLO contribution with m-parton kinematics is given by 



a 



NLO {m} 



da + da 



B 



e=0 



d& m) {\M m { P X,..., Pm )\\l-loop) +m< 1, 



, ml /( 



m > 



(7.36) 

;Pm) ■ 



The first term in the curly bracket is the one-loop^ renormalized matrix element square 
to produce m final-state partons. The second term is obtained by inserting the colour- 
charge operator of Eq. ( 7.26 ) into the tree-level matrix element to produce m partons as in 



"Remember that, according to our calculation of the insertion operator I, the one-loop matrix element in 
Eq. ( 7.36|) h as to be evaluated by using conventional dimensional regularization. We refer to the discussion 
in Sect. 3.3 for the use of different regularization schemes. 
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Eq. (|7.25|) (see also Appendix C). These two terms have to be first evaluated in d = 4 — 2e 
dimensions. Then one has to carry out their expansion in e-poles (the expansion for the 
singular factors Vj(e) is recalled in Appendix C), cancel analytically (by trivial addition) 
the poles and perform the limit e — > 0. At this point the phase-space integration is carried 
out in four space-time dimensions. 
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8 Jet cross sections with one initial-state hadron 



Sections |S|-|TT| are devoted to the generalization of the results of the previous Section to 
processes with identified hadrons (cfr. Sect. |2.3|) . In each of these Sections, we first describe 
the implementation of our subtraction procedure by following closely (although with less 
detail) the steps of Sect. [7|. We start by defining the jet cross sections for each class 
of processes, then we introduce the explicit expression for our subtraction term da A and, 
finally, we perform its integration, calculate the appropriate combinatorial factors, and show 
how the ensuing contribution can be combined with the virtual term da and the collinear 
counter-term da c to provide a finite NLO partonic cross section. The final formulae of our 
algorithm are summarized at the end of each Section. 

Let us start by considering hard-scattering processes with a single incoming hadron 
(cfr. Eqs. (|6.3|- |6.6| )) like, for instance, deep- inelastic lepton-hadron scattering. In the case 
of unpolarized scattering, the Born-level partonic cross section with one parton of flavour 
a and momentum p a in the initial state has the following expression in terms of the QCD 
matrix elements 

dcr a(Pa) = Mn 7 \ ) w T #m(Pl, ■■;Vm\Va + Q) 



n s (a)n c (a)<$>(p a ) ^ } "' ' S {m} 

■ \M m ,a(pi,...,p m ]p a )\ 2 F { J n \p l ,... 

1 Pirn Pa ) > V / 

Here the factor l/(n s (a)n c (a)) accounts for the average over the number of initial-state 
polarizations and colours and &(p a ) is the flux factor. Since p 2 a = 0, the flux factor fulfils 
the following scaling property 

HVPa) =VHPa) ■ (8.2) 

The jet function Fj m \pi, ...,p m ;p a ) is infrared and collinear safe (see Eqs. (|7.2|-[7.4|)) 
and, moreover, it fulfils the property of factorizability of initial-state collinear singularities. 
From a formal viewpoint this implies that 

Fj (Pi, -,Pi, -,Pn+i;Pa) -> Fj^iPl, -,Pn+l,XPa) , ^ Pi ~> (1 ~ %)Pa (8.3) 

for any number n of final-state partons (the n-parton jet function on the right-hand side is 
obtained from the n + 1-parton function on the left-hand side by removing the parton i) 
and 

F < f } (p u ...,p m -p a ) ^0 , if Pi-Pa^O (8.4) 
for the leading-order process (i.e. n = m). 

All the other factors in Eq. (|8.1|) are analogous to those in Eq. (]7[ 



8.1 Implementation of the subtraction procedure 



In order to compute the NLO cross section in Eq. flB.5| ), we write the following identity 

Jm+1 K J 

\ d°a{Pa)+ da^(p a )+ da°{p a ]H 2 F ) 

Jm+1 Jm Jm 



+ 



•5) 
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where, according to the dipole formulae in Sects. 5.1-5.3, a local counterterm da^(p a ) is 
provided by: 

da a(Pa) = Mn „ \ Y d( Pm+l(Pl, -,Pm+l,Pa + Q) 



n s (a)$(p a ) S {m+ i } 

H V H,k(jPX, -,Pm+i;Pa) ^J^bl, -Pij,Pk, ■;Pm+i;Pa) 

pairs k^i,j 
ij 

• \Pij > •• j Pm+1 j Pa) 



pairs 



+ I] Y, V k(Pl,-iPm+liPa) Fj (Pl,~Pk,~,Pm+l,Pa 

i kjti 



(8.6) 



Here V ij>k (p 1 , ...,p m+1 ), X>y(pi, ...,p m+ i;p a ) and £>jf (pi, ...,p m+ i;p a ) are respectively the 
dipoles in Eqs. Q, (|5.36|) and (|5.61|) and F} m \p u ..pij,p k , ..,p m +i;p a ), F} m \p u .Jp ih .., 
p m +i;p a ) and Fj m \pi, ..p k , ..,p m+ i;p a j) are the jet defining functions for the corresponding 
m-parton states {p u ..pij,Pk, ..,p m +i;p }, {pi, -Pii, -,Pm+i;Pa} and {pi, ..p k , ..,p m+1 ;p a i}. 

Since the dipole contributions exactly match the soft and collinear divergences of the 
square of the matrix element |.M m+ i ia | 2 , the subtracted expression (da^(p a ) — da^{p a )) in 
Eq. (|S.5|) is integrable in d — 4 dimensions. 



In order to compute the contribution in the square bracket of Eq. ( |8.5| ), we write da^{p a ) 
as follows 

dffi(p a ) = dafiPa) + daf'iPa) + da?"( Pa ) , (8.7) 

where the three terms on the right-hand side of Eq. ( |8.7| ) are in one-to-one correspondence 
with those in the curly bracket on the right-hand side of Eq. ( |3.6D . 

The integration of da^'{p a ) can be performed analogously to that of da A in the previous 
Section, thus leading to the following result 

f a f 1 1 

/ da a '(p a ) = - Mn , w . Y d(j>m(pi, ~>Pm\P* + Q) TT~ 
Jm+1 Jm nAa)<2(Vn) 7^ Ofr, 



{m} 



j Pm j Pa )YY\M i m k a(PU- i Pm i Pa ) 



a s 1 



{m} 



Vi(c) 



i 27rr(l-e) \2 Pl . Pk J T\ 

where the functions V*(e) are defined in Eqs. (|7.21| , |7.22"|) . 

Let us now consider the integration of da^"{p a ). We first rewrite Eq. (|5.36|) as follows 



> Pm+li Pa) 



V" 1 .. _ 

7^2" l-^m,a(Pl ; "Pij, -,Pm+i;Pa)| 2 



^Xij^a Pi • Pj T i 



19) 



■J /i 



where the notation is similar to that in Eq. ( |7. 1 1| ) , and then, using the phase space convo- 
lution properties in Eqs. Q5.46J5.47 ), we obtain: 

/ daf'ipa) = - [ M in Y Y [ dx d(f) m (pi, ..pij, ..,p m+1 ;xp a + Q) — 

J m+1 Jm t m+ n pairs JO 0{m+l} 
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1 

x 



n s (a)$(p a ) 
\M%? a (pi, -Pij, ..,p m +i;xp a )\ 2 J 



2pi ■ pj 



(8.10) 

[dpi{pij\Pa,x)\ . 



As in the case of the splitting functions Vij,k, the azimuthal correlations due to V?- vanish 



after integration over p$ (keeping p^ and x fixed) and we find 

/ daf( Pa ) = - [ M in E t dx ( w 



•d<f> m (pi, .-Pij, ..,p m+1 ;xp a + Q) 
■ \M^ a (p!, ..Pij, ..,p m+1 ;xp a )\ 



Fj\pi,..Pij,..,p m+ i;xp a ) 



S{m+1} 
2 "S 1 



477/i 2 



27rr(l - e) 



I rp2 Vij(, X ' € ) ) 
Pa J -L ij 



(8.11) 



where the functions Vij(x;e) are given in Eqs. (|5.57| - |5^59|) and we have used Eq. (|8.2| ) to 
replace x<&(p a ) with $(;rp a ). 

Equation (|8.11 ) is similar to Eq. (|7.13| ) in Sect. |7|, apart from the replacement k — > a 
as spectator parton. In order to rewrite Eq. ( |3.11|) in terms of an m-parton contribution 
times a factor, we have to perform the counting of the symmetry factors for going from m 
partons to m + 1 partons in the final state. Since, as shown in Sect. [7], this counting is 
independent of the spectator parton, we immediately get 



/ daf(p a ) = - ( fV m I 

Jm+l Jm JO 



dx 



o n s (a)$(xp a ) 



d<t> m (pi, •• ) Pm ) %Pa + Q) 



j Pmi Xp a i Pmi Xp Q 



S{m} 

a s 1 



/ 47T/X 2 ' 

27rr(l-e) \2pi-p at 



Vi(x; e) 



where we have defined 



V l {x\e) = V qg (x;e) , 



if i — q,q 



V l {x; e) = - V gg {x; e) + N f V q q(x; e) , if % = g . 



(8.12) 

(8.13) 
(8.14) 



Recall that all the e-poles of V{(x; e) are accounted for by 5(1— x)Vi(e) terms (see Eqs. (|5.57| - 
5.59| )), where K(e) are the functions defined in Eqs. ( |7.21| , |7.22|) . 

Let us now consider the integration of da^"'(p a ). We first rewrite the corresponding 
dipole contribution (see Eq. ( |5.61| )) as follows 



£T(Pl> -iProfliPa) = - 



1 



2^jfc,a Pi ' Pa T 



m+1 ) Fat i 



(8.15) 



Thus, using the phase space convolution in Eqs. ( |5. 70 J5TTID and performing the integration 
over pi (keeping pk and x fixed), the azimuthal correlation due to V™ vanishes and we 
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obtain 



/ dafiPa) = - / Mn C dx £ £ £ ( ~L \ 
Jm+i Jm Jo {m+1} ■ k ^n s {ai)<$>{xp a ) 



^j m) (Pi,-,Pfc,--,P m ;a;p a ) 



^{m+l} 

|M"'~(pi, -Pfc, ..,p m ;a;p a )| 2 °" 



27TT 1 - e 



Ann 2 

Ipk ■ Pa 



-^V a > ai {x;e) , (8.16) 



where the functions V a,at (x; e) are given in Eqs. (|5.81 - |^53 ). 



The right-hand side of Eq. ( 8.16p can easily be rewritten in terms of a sum over m-parton 
configurations. To this end, we have to perform the corresponding counting of symmetry 
factors, which, nonetheless, is trivial in this case. If the Born-level m-parton configuration 
has rrii partons of type i, the corresponding m + 1-parton configuration has m^ + 1 partons 
of the same type, so that: 



E v _ 

{m+l} °{m+l} i 



£■■■ = E E 



... mil 



1 



{m+l} i 



...(rrn + iy. S. 



{m} 



(8.17) 



However, there are mj + 1 possible ways of choosing the parton i in the m + 1-parton 
configuration and, hence, we obtain 



E y^E 

{m+l} °{m+l} i 



E 



... mil 



rrii + 1)! 
1 



[m. 



+ i)E 



{m} S{m} 



EE, 

i {m} °{ m i 



lit 



where the J2i i n the last line of Eq. (|8.18|) simply denotes the sum over the flavours i in the 
m-parton configuration. We can thus rewrite this sum as a sum over the flavours b = ai in 
the initial state and Eq. (|8.16f) becomes: 

f daf"(p a ) = -[ J\f in J2 I dx d<j) m (pi,..,p m ;xp a + Q) - — 

J m+l Jm {m} &{m} 



Fj m \pi,..,p m ;xp a ) ££ 



\M b A(p l: ...,p m ; xp a )\' 



k b 



1 / 47T/X 2 ' 

27rr(l-e) \2p k -p ai 



n s (b)®(xp a ) 



V a ' b (x:e) . 



(8.19) 



Note that, in addition to 5 ab 5{l - x)V b {e) terms like those in Eqs. (|8.13| , |8.14|) , V a ' b (x; e) 
contains e-poles in terms of the form P ab (x)/e (see Eqs. ( |5.8lH5^84l) ). 



Collecting Eqs. (M , CT , g7i9|) and adding Eq. (gjf), we find 

f daa(Pa)+[ dof(p a ) = -£ I dx [ AfinY,d4>m(pi,~,p m ;xp a + Q) 

J m+l Jm u JO Jm rl 



{m} 



Fj^iPl, ..,Pm',Xp a ) 



S{m} 



n s (b)$(xp a ) 27rr(l-e) 
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+ 5 ab Y: \M% b ( Pl ,..., Pm ;x Pa ) 

i 



47T/1 



2 \ e 



1 



\2p t - P aJ T\ 
47T/i 2 V 1 



T 2 



+ 



n c (b) 



\M mth {pi, ... 



,2pi -p a 

1 Avr/x 2 



^2 Vi(x;e) 

i 

V a ' b (x; e) 



, , P a6 (x) + Kf s (a;) 



(8.20) 



We see that J m+ ida^(p a ) + J m da^(p a ) is obtained from the leading-order expression 
J m da^(xp a ) by replacing the Born-level matrix element squared 



n. a 



i,a^ ••■) Pa | ••■) Pa ^m,a 



by 



£ 



6 < ...;xp | i" a ' 6 (x;e) \...;xp a > mfi 



n s (b) 



(8.21) 



(8.22) 



and performing the x-integration. Here the insertion operator I(x; e) depends on the colour 
charges, momenta and flavours of the QCD partons 



r' b (p 1 ,..., Pm ;p a ,x;e]ii 2 F ) 



27rr(l - 



^(l-z)££lVT fe (^-Y iv^ + ^ET^Ti ' l7r/r> 



2 Pi - Pk T\ 



2 Pi -p a l T\ 



V,(x; e) 



^8.23) 



This form of the insertion operator can be simplified in the limit e — > 0. We start by 
rewriting Eq. ( |8.23|) as follows 



where 



i" a ' (p 1 ,... }Pm ;p a ,x;e;iJ, F )=5 a 5(1 - x) I(pi, ...,p m ,p a ; e 



+ S ab I {1) ( Pl ,...,p m ;p a ,x;e) + -,P m ;Pa, x\ e; /4) 



a,b i 



I{p 1 ,...,p m ,p a ;e) = 

/ 47T/i 2 



+ T l -T a 



2^r(l-e){^T 2 Vi(e) 



1 



2Pi-Pa, 



(8.24) 



(8.25) 



I m {px,...,p m ;p a ,x;e) 



a s 1 



£2VT a 



1 



27rr(l-e)Y ' " Wi-PaJ Tt 
Vi (x;e)-5(l-x) Vi(e)l , 



(8.26) 
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■^(2)(Pl' ■■■iVm'-iVai £\ l^p) 



1 



V a >\x-e)-5 ab 5(l-x) V a (e) 



1 /47T/i 2 V 

e 



P afe (*)+iCa(* 



47T/i 2 
06 , 



(Tl2 

1 b 



(8.27) 



The operator I(pi, ...,p m ,p a ; e) in Eq. ( |8.25| ) is fully symmetric with respect to the 
dependence on colour charges and momenta of the m + 1 partons {pi, ...,p m ,p a }. In other 
words, this operator is identical (apart from depending on the additional initial-state parton 
a) to that in Eq. ( |7.26| ). Since crossing the momentum of partons from the final to the 
initial state does not change the singular terms in the one-loop QCD amplitudes, it follows 



that the insertion operator ( |3.25|) cancels all the singularities in the virtual contribution 
f m d(7a (j>a). As a consequence, the two other operators J(i) and J( 2 ) should contribute as 
finite counterterms. Actually, 1^ and 1(2) are separately finite in the limit e — > 0. 



In order to show that, we can use Eqs. (|53^ - [539| ) , ( |8~T3| , |8TT4] ) and ( f?T21]f?T22D and thus 
obtain 



V l {x-t)-5{l-x)V i {e) 



2 m 1 



1 — X 1 — X 



(t^) + + ^-*> 



1 — X 
+ 0(e 



ln(2 - x) 



where 7^ are given in Eq. ( [5.901 ). Then we can write: 



I(i) (pi, ...,p m ;p a ,x; e) 



+ — ln(2-ar)-^ 
1 - x V ; T 2 



— In 

x I — x 

+-0(6) . 



(8.28) 



(8.29) 



Coming to 7(2) , let us rewrite Eq. ( |8.27| ) in the following form 

1 < ' ^ 



T a,bf 2 \ _ a S 



"(2) 



27rr(l-e) 



47T/T 



V a '"(x; e) + - P af, (x) - <T5(1 - x) V a (e) 



'47T/! 2 



/4 



Ip^(x) + ^f a (x) 



5.30) 



Using Eqs. (|5.81 - 5.84 ), we see that the first square bracket contribution on the right-hand 
side of Eq. ( p.30| ) is finite for e — > and given by 



V a ' b (x; e) + - P a \x) - 5 ab 5(l - x) V a (e) = K a \x) + P a \x) In a 



2 1 \ 2 

In + ln(2 - x) 

1 — x 1 — xj + 1 — X 



+ 0(e) 



(8.31) 
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-u6 



where the K (x) functions are defined so as to simplify the final formulae, 



K g9 (x) 



K 9q (x) 



K qq (x) 



K"(x) 



P q9 (x)\n- — - + C F x , 

x 

P 9q (x) In — - + T R 2x(1 - x) 

x 



C f 



— ^— In- — -(l + a;)ln- — - + (1 - x) 

1 — X X J + X 



5(1 - x) (5 - 7T 2 ) C 
2C A 



1 1 - x\ fl — x „ _ A, 1 - x 
In ) +( l + x(l-x) In 



1 — x x 



a; 



5(1 -x) 



9 / 9 * " 



(8.32) 
(8.33) 

(8.34) 
(8.35) 



As for the second square bracket contribution on the right-hand side of Eq. (|8.30|) , using 
colour conservation (X^Tj = —T b ) and expanding in e, we get an (9(e) term: 



b m2 



47TyU 2 

2pj -p a . 



1 - /4 



4 



^E^.-T fc -,ln-p- + 0( e 2 ) . (8.36) 



Inserting Eqs. ( |8.31| . |8.36| ) into Eq. ( |8.3(J| ), adding the Jm contribution in Eq. ( |8.29| ) and 
again using colour charge conservation, we end up with 



$ a I (i)(Pi, -,p m ;Pa,x; e) + / J) (pi, ...,p m ;Pa,x] e; fi F ) 
= K a ' b (x)+P a > b (p u ...,p m ,xp a ,x;^ F ) + 0(e) , 

where we have defined 



li 



K a ' b (x) = J± K a \x) - K*{x) + 5 ab £ Ti ■ T b7 ^ 



1 — x 



+ 5(1 -x) 



(8.37) 



(8.38) 



noft, N 1 „ , ftp 



P a > b ( Pl , ...,p m , xp a , x; &) = ^ P ab (x) — 2 £ T 4 • T b In 



2xp a • 



5.39) 



The final result for I a,b (x; e) is the following 

I a ' b {pi, ...,p m ;p a ,x; e; /4) = 5 ab 8(1 - x) I( Pl , ...,p m ,p a ; e) 
+K a > b (x) + P a ' V, ...,p m ; xp a , x; n F ) + 0(e) , 



(8.40) 



where the insertion operators I,K a ' b and P a,b are given in Eqs. (|8.25 ), ( 8.38 ) and ( |8.39|) . 
Therefore, using a notation similar to that in Eq. (|7.23|) , we can write 



d^( P ) + I daC(p;v F )= f [da B a (p) ■ 1(e) 

-1 Jm Jm L 



(8.41) 



+ E / dx / \ Ka ' b ( x )- d(7 bM]+J2 dx / \P a,b (xp,x;fi 2 F ) -da^xp) 



Jm 



Jm 
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Note that all the insertion operators 1(e), K a ' b (x), P a,b (xp, x; fip) depend on the colour 
charges and flavours of the QCD partons. However, while this dependence is fully symmetric 
in 1(e), the operators K a ' b (x) and P a ' b (xp, x; ji 2 F ) do depend asymmetrically on the flavour 
and colour charge of the incoming parton p. In addition, 1(e) depends on the parton 
momenta, K a ' b (x) depends on the momentum fraction x (but not on the parton momenta) 
and P a ' b (xp,x; ji 2 F ) depends on x, parton momenta and factorization scale. 

As in the case of processes with no initial-state hadrons, the term da^(p) ■ 1(e) in 
Eq. flS.41| ) cancels all the e-poles in the virtual contribution da^(p), thus making the NLO 
cross section in Eq. (|S.5|) finite. The other two terms on the right-hand side of Eq. Q8.41) are 



finite remainders that are left after factorization of the initial-state collinear singularities 
into the parton densities. 

The operator P a,b (xp, directly related to the Altarelli-Parisi probabilities. In 

particular, from Eq. (j8.39|) we can check that it fulfils the relation: 

gg^l^p£iii) = | p- (x) i^r, . Tb= _g p- (x) , (8 .42) 

where we have used colour-charge conservation (X^Tj = —T b ). It follows that it cancels 
the similar (and with opposite sign) factorization-scale dependence of the parton densities 
f a (x,ii F ) thus making the hadronic cross section ( |6. 3[ ) /z^-independent to NLO accuracy. 

The operator K a,b (x) contains ( ) + and 5 distributions with coefficients 7^ of the same 



type as those appearing in the singular operator /(e) (see Eqs. ( |7.27y8.25|) ). These terms, 
due to colour correlations between the incoming parton and the final-state partons, are the 
heritage of the initial-state collinear divergences originally present in the real contribution 
do~a(p)- Moreover, K a,b (x) depends on the flavour (and colour diagonal) kernels K^ & (x) 
and K (x). The kernel K^\(x) is related to the definition of the factorization scheme 
(cfr. Eq. (|6.6|) ) while K (x) has a close relationship with the parton splitting functions. 
As a matter of fact, if we define Pj b (x) by expanding in e the <i-dimensional Altarelli-Parisi 
splitting functions in Eqs. (f4.18H4.2l"|) : 



< P ab (x; e)> = < P ab (x; e = 0) > - e Pj b (x) + 0(e 2 ) , (8.43) 
we can rewrite Eqs. fl3.32| - j3T3"5| ) as follows 

K ab (x)=P n ' h (x)+P r t(x) In- 

+ 5° 



x 

r / o 1 _ ^\ 

\-o6 



T "(T^'^) -*(l-*)(7. + *.-^H 



(8.44) 



where P^ b g (x) is the non-singular part of the (four-dimensional) Altarelli-Parisi probabilities 
(see Eq. (|5.89f) ) and the coefficients K a are defined in Eq. (|7.28|) . The contribution Pj b (x) 
on the right-hand side of Eq. ( p. 44 ), is directly related to the dimensional regularization 



of the initial-state collinear singularities. Indeed, it comes from the interference of the 1/e 
collinear pole and the (9(e)-contribution to the (i-dimensional splitting functions. Also the 



other terms on the right-hand side of Eq. ( |8.44 ) have a simple interpretation. After having 
factorized the initial-state collinear divergences, the finite remainder K (x) is proportional 
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to the phase space available for the emission from the incoming parton. The factor ln(l — 
x)/x multiplying P^ g (x) has this kinematic origin. The same factor controls radiation that 
is both collinear and soft and thus it enters into the ( ) + -distribution in the square bracket. 



8.2 Final formulae 

Summarizing the results derived in the previous Subsection, we obtain the following final 
formulae for jet cross sections involving one initial-state hadron. 



The partonic cross section on the right-hand side of Eq. ( |6.3| ) consists of a LO and a 
NLO component. In the case of an incoming parton of flavour a and momentum p a , the 
explicit expression for the LO component is: 

°a°{Pa)= I da*( Pa )= !d& m \p a ) \M mA {p l) ..., Pm]Pa )\ 2 F ( r\p 1) ... )Pm] p a ) , 

Jm J n c \0,) 

(8.45) 

where M. m ,a is the tree-level matrix element to produce m final-state partons, Fj^ is the jet 
defining function (it fulfils the properties in Eqs. (|8.3|J8.4| ) in addition to those in Eqs. (|7.2| - 
7.4|)), the factor l/n c (a) comes from the average over the colours of the initial-state parton 



and d&( m \p a ) collects all the remaining factors (phase space, flux, spin average) on the 
right-hand side of Eq. ( |8.1| ). The evaluation of the LO cross section ( |8.45| ) is carried out in 
four space-time dimensions. 



The NLO partonic cross section is split into three terms, as in Eq. (|2.15|) : 



^ LO (p a ,ti)=^ LO{m+1} (Pa) + a^ LO{m} (Pa)+ / dxa N a L0 ^\ X -xp a ^l) . (8.46) 

Jo 

The term with m + 1-parton kinematics is given by 

a NLO {m+ i }{pa) = 1^ \(da*( Pa )) e=o - ( £ da*{p a ) ® ( dV dipole + dV d ' ipc 

ydipoles 

• , Pm+1 ;Pa)\ 2 Fj m+1) fa , . . . , p m+ i ; p a ) 




e=0 



- E (V-F^)( Pl ,...,p m+l]Pa )\ , (8.47) 

dipoles 

where Ai m +i,a is the tree-level matrix element with m+1 partons in the final state and 
Edipoles (V-F^)( Pu ... , Pm+i', Pa) is the sum of the dipole factors contained into the curly 
bracket on the right-hand side of Eq. (|S.6|). Note that the m + 1-parton matrix element 



is multiplied by F^ m+1 \ the jet function for m + 1 final-state partons, while the dipole 
contributions involve the m-parton jet function Fj m \ All the terms in Eq. ( |8.47| ) are 
evaluated and integrated in four dimensions. 

The NLO contribution with m-parton kinematics is exactly like that in Eq. ( |7.36| ) for 
e + e~-type processes, apart from the additional dependence on the colour and momentum 
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of the incoming parton. Indeed, combining the virtual cross section with the first term on 
the right-hand side of Eq. ( |8.41| ), we obtain 



d°l (Pa) + d° a iPa) 



e=0 
2 



7 T \M m: a(pi, ■■■ j Prri) Pa ) \ (\—loop) 

+ m , a < 1, ....,m; a\ 1(e) |1, ....,m; a > m>a [ F ( f l \pi 1 ... 

i Pm j Pa J 

J <E = 



(8.48) 



where \M. m ,a\fi-i oop ) is the one-loop matrix element square and the colour charge operator 
/(e) is given in Eq. (|8.25|) (see also Appendix C). The two terms in the curly bracket 
have to be separately computed in d = 4 — 2e dimensions and the limit e —>■ performed 
after having cancelled analytically the e poles. At this point the phase space integration is 
carried out in four dimensions. 

The NLO component involving the one-dimensional convolution with respect to the 
longitudinal-momentum fraction x is given by the last two terms on the right-hand side of 
Eq. flOU) : 



Jo 



dx ^NLO{m} ( 



%Pa , ^F) 



E 



dx 

J m 



daf(xp a ) 



K 



P) a ' fe (x) 



e=0 



= W 1 dx [ d& m \xp a ) Fi m \ Pl ,...,p m ;xp a ) 

• mfi < 1, ....,m;xp a \ {K a ' b (x) + P a ' b (xp a , x; /i|)) |1, m; xp a > TOi6 . 

The colour-charge operators K and P are respectively defined in Eqs. ( fcj.38| ) and 
Their explicit expressions, as well as those of the related flavour kernels P ab (x), K (x) and 
K® b & (x), are recalled in Appendix C. The calculation of Eq. ( ^.49| ) is directly performed in 
four space-time dimensions. 



.49) 

!)• 



The partonic cross sections in Eqs. ( |8.45| - p749D have to be convoluted with the parton 
densities as in Eq. ( |6.3|) , in order to compute the corresponding hadronic cross sections. 



Note that, because of the x-dependence of the operators K and P in Eqs. ( |8.38| , |8.39| ), 
the cross section component in Eq. ( |3.49|) is boost-invariant with respect to the direction 
of the incoming momentum p a . Therefore, in the evaluation of the hadronic cross section, 
this contribution enters in the form of multiple convolution of a Born-type partonic cross 
section with the kernel K or P and with the parton densities. 



68 



9 Jet cross sections with one final-state identified hadron 



Let us now consider fragmentation processes, starting with the simplest case, which does 
not involve initial-state hadrons. According to the definition in Sect. |^, the inclusive cross 
section to produce a parton of flavour a and momentum p a has the following expression at 
the Born level 

dcr fincl)a (Pa) = Mn Z (Ph -, Pm] Q ~ Pa) 

{m} 

■ \M m +a{Pa,Pl,--,Pm)\ 2 F^j m) ([p a ],p 1 ,...,p m ) . (9.1) 

The notation in Eq. ( |9.1|) is similar to that in Eqs. ( |7.1|) and (|8.1|) . The only relevant 
difference regards the jet defining function Fj m \[p a ],pi, ...,p m ). Besides being infrared 
and collinear safe (see Eqs. (|7.2|- [7.4| )), it should guarantee the factorizability of final-state 
collinear singularities. This implies the following general (i.e. for any number n of partons) 
property 

f j (\Pa],Pi, -,Pi, -,Pn+i) ->F} n \\p a /z],pi,...,p n+ i) , if pi -»• {l/z - l)p a , (9.2) 
and the leading-order constraint 

F!j m \ip a },p u ...,p m )^0 , if Pi-pa^O ■ (9.3) 



9.1 Implementation of the subtraction procedure 



As usual, in order to evaluate the NLO partonic cross section, we rewrite Eq. ( |6.16| ) in the 
following form 



a (ind)a(Pa) ~ \ (d^Uncl) a(Pa) — d^find) a(Pa) 
-'m+l x 



(9.4) 



+ 



/ daf ind)a (p a ) + daY incl)a (p a )+ daf md) a (p a ; fi 2 F ) 



where, using the dipole formulae in Eqs. ( |5.92| ) and ( [5.113D , we define the local counterterm: 

[pi, ...,p m+1 ; Q — p a ) — 

{m+l} 



dcrfincl)a(Pa) = Mn Z #m+l(Pl, Pm+i; Q - Pa) - Q 

'-'{m+l} 



Z Z V ijAPa,Pl, -,Pm+l) F} m \\p a ],pi, ..pij, ..p k , ..,Pm+l) 

pairs k=£i,j 
i.j 

+ Z V ijAPa,Pl, -,Pm+l) ^/([Pa],Pl, -Pij, -,Pm+l) 



pairs 



Z J2 V ai,k(P<nPl' ->Pm+l) F j(\Pai\,Pl, -Pk, -,Pm+l] 



(9.5) 



i k^i 



The four-dimensional integrability of (da^ ncl ^ a (p a ) — d<jf incl ) a (p a )) follows in exactly the 



same way as for Eqs. ( |7.5|) and flS75] 



(incl) ( 
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To compute the term in the square bracket of Eq. ( |9.4|) , we write da A ncl y{p a ) as follows 

drfincl)a(Pa) = d °(incl)a(Pa) + d(T A " d)a (Pa) + d(J A '" d)a {j> a ) , (9.6) 

where the three terms on the right-hand side are in one-to-one correspondence with those 
in the curly bracket on the right-hand side of Eq. (|9.5| ) . 



The integration of da^ ncl ^ a {p a ) can be carried out analogously to that of da A in Sect. [7| 
and, thus, we obtain 

/ daf- ncl)a (p a ) = -f MnZ)d0m(pi,.-,Pm;Q-Po) 
Jm+1 Jm r | >J{m} 

■ F j(\Pa],Pl, -,Pm)J2J2 \ M m+a(Pa,Pl, -;Pm)\ 2 



as 1 I 47T/1 2 



^2 Vi(e) , (9.7) 



2vrr(l -e) \2 Pi -p k J T 

where the functions Vj(e) are defined in Eqs. ( |7.21|j7T2"2|) . 

The integration of da A '^ d ^ a {j> a ) is similar to that of da A "(p a ) in Sect. ||. The main 
difference is that the phase space convolution in Eq. ( p.lUl| ) replaces that in Eq. 
Thus, we find: 



/ dau'n d )a{Pa) = ~ / Mn / "4^7 E d <Pm(pi , •• , Pm\ Q ~ Pa/ Z ) 

J m+1 Jm JO Z , _ 



{m} 

Fj m) (\Pa/z],Pl,~,Pm)Yl \ M m+a(Pa/z,Pi,.. } p„^ 2 



J m) 



S{m} 

a s 1 ( Arrfi 2 V 1 ^ 



where we have defined 



Vfae) = V gg (z;e) , if % = q, q , (9.9) 

Vi(z;e) = i V M (z;e)+W,V OT -(z;e) , if i = g . (9.10) 



The integration of da A '" d ^ a (p a ) is again analogous to that of da A "'(p a ) in Sect. ||, apart 
from the different phase space convolution in Eq. ( |5.122| ) for the dipole {ai, k}. We find: 



da Al " 



m+1 



(incl)a 



(Pa) 



[ MnE / 

Jm r -, JO 

{m} 



1 dz 



Z 



2-2e 



d(j)m(Pl, -,P m ; Q -Pa/Z) 



S{m] 



F J (\Pa/z\,Pl,..,Pm) \M^ +h {pJ Z,p X , ...,p T , 

k b 

' 1 



tt's 



1 / 4tT/^ 

27rr(l-e) \2p k -p a 



-a Vb,a{z\ e) 
J- k 



(9.11) 



where the functions V a , a i(z] e) are given in Eqs. ( p,131| - [5\134j ). 
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Combining the results in Eqs. (9.7,3.8, 9.11 ) and adding the collinear counterterm in 
Eq. ( |6.17| ), we have 



/ daf incl)a (p a )+ [ daf incl)a (p a ) = -J2 I [ ^ 

Jrn+l y ' Jm y ' JO Z Jm 

(Pi, - Q-Pa/z)- Fj m) ( \p a /z],pi, ..,Pm) 77^^71 \ 

{m} S {m} 27rr(l-ej 



5 ab 5(1 ~Z)J2J2 \ M m+b(Pa/z,PU-,Pn 

( 4tc[i 2 



+ 5 ab J2 \M l m+b (p a /z,p u ...,p r , 

i 

+ \^m+b{Pa/z,Pl,...,Pm)\ 2 



+ \M m +b{Pa/z,Pl,--,Pr, 



\2pi-p a ) T\ 



2 (s h v ^ 

1 _ 

=2 v i(^; e ) 



\2pi "Pa 
1 /47T/1 2 



(9.12) 



Hence, Xm+i df7 (Ld)a(Pa) + J m da find)a(Pa) is obtained from the leading order expression 



(incl) 



fmd°~finci)a(.Pa) by replacing the Born-level matrix element squared 



m+a 



<Pa,-\Pa, - > 



m+a i 



by 



m+b< Pa/Z, ....| Ib,a(z; e) \Pa/z, ... > 



(9.13) 



(9.14) 



where the insertion operator I(z; e) depends on the colour charges, momenta and flavours 
of the QCD partons: 



* b.a 

(P a; Pit ■■■■> Pmi #i 



2vrr(l - e) 



- 2 ne) + S ab Y,T,-T b 



2 \ 6 



2pi-Pa T 



' Vi(z:e) 



2 K IV*) 



(9.15) 



This insertion operator is similar to that in Eq. (|8.2 



for the cross section with a 
single incoming parton, apart from the replacements Vj(x;e) — > Vj(z;e), V a,fc 



V, 



b.a 



z:e 



Pab(x) -> AaU), # 



a 6 



Hi 



6a 



(note, in particular, the transposition of 



the flavour indices that is involved in these replacements). Therefore, using the analogues 
of Eqs. ( gggggl[ ), namely 



V i (z;e)-8{l-z)V i {e) = li + T\{— In ' 



1-z 1 - z 



- li 



+ 0(e) 



(9.16) 
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V M 0; e) + ~ P ha (z) - M(l - *) V„(e) = i^ a (z) + 2P ba (^) ln^ 
- ^(^In^) +0 (e), 



(9.17) 



and performing the same algebraic manipulations as in Sect. |8|, we end up with the final 
result: 



m+l 



d<?tind)a(p) + daf ind)a (p; fi 2 F ) = Ida f incl)a (p) ■ 1(e) 

Jm Jm L 



(9.18) 



+ L ^ I daf incl)b {p/z)-H bA (z)\+Y, -J / daf incl)b {p/z)- P bA (p/ Zl z]n 2 F ) 



Z Jm 



^ Jm 



where the insertion operator /(e) is exactly the same as in Eq. (|8.25| ) and the insertion 
operators H b)CL (z) and P^ a (p/ z, z; fi F ) are defined as follows 



Z7T L 



7< 



1 



1 - z 



6(1 -z) - 1 



(9.19) 



Pb,a(Pl, ■■■,P m ;Pa/z,Z;/J 2 F ) 



z/sp 



^ P b ( z ) J_y T , r T b In 
2tt bal Tfe Y 2p a - Pi 



(9.20) 



Equation (|9.18 ) is the time-like (a single identified parton in the final state) analogue 
of Eq. ( |8.41|) for the space-like (a single parton in the initial state) case. The contribution 
a (p) ■ 1(e) cancels all the e-poles in da^ incl ^ a (p) thus making the NLO cross section 



d^Lci) 



(in t 



finite in the four- dimensional limit. 



The operators H bA (z) and P bjCl (p/z, z\ fip) are instead finite for e — > (for this reason, 



in Eq. (|9.18 ) we have replaced the phase space factor dz/z 2 2e of Eq. (|9.12j ) with dz/z 2 ) and 



are similar to the operators K a,b (x) and P a,b (xp,x; ^ 2 F ) in Eq. ( 8.41| ). Actually, apart from 
the momentum rescaling xp — > p/z, the only other difference between the two P operators 
is in the transposition of the flavour indices a and b. Therefore, in spite of the identity 
Pab( z ) — P ab (z) of the Altarelli-Parisi probabilities for time-like and space-like splittings, 
in the case of the operators P we have P b ^ a (p, z; fi F ) ^ P b,a (p, z; [i 2 F ). This difference is due 
to the colour correlations Tj ■ T b or, more precisely, to the fact that momentum fraction 
and colour flow in opposite direction in the time- like and space-like cases (Fig. 4). 



Comparing Eq. fl3.38|) and ( |9.19|) , we see that this transposition of the flavour indices 
also affects the difference between K and H. Moreover, H differs from K by an extra 
term 2>P ba (z) \nz — 6 ab Y,i Ti ■ T b 'Ji/T 2 , which can be attributed to the kinematic crossing 
from initial- to final-state partons. The appearance of the two different kernels K® b a and 
H b f is trivially related to the choice of the factorization scheme of collinear singularities 
(cfr. Eqs.(p|) and (gl?D ). 
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Space - like 



Figure 4: Altarelli-Parisi insertion operators for the time-like (Eq. ( |9.20| )) and the space-like 
(Eq. (ET3SD ) cases. 



9.2 Final formulae 



The calculations carried out in Subsection |9.1| for jet cross sections with a single iden- 
tified hadron in the final state lead to results that are very similar to those described in 



Subsection 3.2 for the kinematically-crossed process with a single hadron in the initial state. 



In summary, the LO parton-level cross section flB.15 ), to be convoluted with the non- 
perturbative fragmentation function as in Eq. fl3.14j ), can be written as follows 



a 



LO 



u„d)a(Pa) = / d °(incl)a(Pa) = / d& m) (p a ) \M m+a (p aj p U . .. , p m )\ 2 ([ Pa ], Pu . .., Pm ) . 

Jm J 

(9.21) 

Here a and p a respectively denote the flavour and momentum of the identified parton, 
M. m +a is the tree- level matrix element to produce m unidentified partons in the final 
state, Fj m ^ is the jet defining function (its general properties are listed in Eqs. (|9.2|J9~ 



and Eqs. (J_^-J_A)) and d& m \p a ) stands for all the remaining phase-space factors on the 
right-hand side of Eq. (|9.1|). 



The NLO partonic cross section is decomposed into three terms. Following the symbolic 
notation in Eq. ( 2.15|) , we write: 



NLO i ...2\ NLO {m+1}, s NLO{m} ( v 

a (incl)a\Va, ^F) ~ a (ind)a \Pa) + ° Unci) a \Pa) + 



(incl) t 



1 dz ^ NLO {m} 



a 



Z" 



MrteP*M) ■ (9-22) 



The NLO contribution with m+l-parton kinematics has the following explicit expression 



a 



NLO {m+1} 

(incl) a 



(Pa) 



m+1 



d& find) aiP a, 



e=0 
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E da (vncl)a(Va) ® ( oVdipole + dV& 



dipole 



i dipoles 



e=0 



J d^ m+l \p a ) { \M m+1+a { Pa , Pl ,...,p m+ i)\ 2 Fi m+1 \\p a },p l ,...,p m+1 ) 



e (v-F^)( Pa , Pl ,...,p m+1 )} , 

dipoles 



(9.23) 



where j\4 m+ i +a is the tree-level matrix element with m + 1 unidentified partons in the final 
state and £ dipolcs (v ■ F^) (p Pi, ...,p m+ x) is the sum of the dipole factors contained into 
the curly bracket on the right-hand side of Eq. ( |9.5|) . 



The NLO term with m-parton kinematics is obtained by adding the virtual cross section 
and the first contribution on the right-hand side of Eq. fl9.18|) . This term is completely 
analogous to that in Eq. (|8.48|) and is given by: 



(mcZ)i \Pa) — J d(J( ind ) a (p a ) + da" inc ^ a {p a ) 

J d& m \ Va ) | \M m+a (p a ,p 1 , ...,p m )\ 2 {1 _ loop) 



e=0 



+ 



< a, 1 ml 1(e) la, 1 m 



(e) |a,l,....,m> m+a j F { J m \[p a },p 1 ,...,p m ) . (9.24) 



where the colour-charge operator /(e) is defined in Eq. (|8.25| ) (see also Appendix C). 

The third contribution on the right-hand side of Eq. (|9.22|) involves the integration of 
an m-parton cross section with respect to the fraction z of the longitudinal momentum 
carried by the identified parton. This contribution is given by the last two terms on the 
right-hand side of Eq. ( |9.18| ) : 



£ % C^W^I) = E £ % I [d*fa)Mz) ®(H+ P) ba (Z) 
JO Z , JO Z Jm L 



- 1 dz 



e=0 



E 



1 dz 



o z 1 - 



d& m \p a /z) F^\\p a /z], Pl ,...,p m ) 



(9.25) 



+b < p a /z,l,... } m\ (H bta (z) + P b , a (p a /z,z;iJ, F )) \p a /z,l,...,m> m+b 



where the colour-charge operators H and P are respectively defined in Eqs. ( |9.19| ) and 
( P-2C| ) (see also Appendix C). 

The actual evaluation of Eqs. (|9.21|) , ( p.23| ) and ( |9.25| ) is directly performed in four 
space-time dimensions. As for Eq. ( |9.24|) , one should first cancel analytically the e poles 
of the one-loop matrix element with those of the insertion operator I, perform the limit 
e — > and then carry out the phase-space integration in four space-time dimensions. 
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10 Jet cross sections with two initial-state hadrons 



In the case of unpolarized scattering, the Born-level cross section with two incoming partons 
of flavours a and b and momenta p a and pb is the following 



1 



n (a)n (b)n (a)n ( bWp . Pf) T, d MPn -,Pm,Pa + Pb + Q) 

1 b s\ Uj } 1 L s\ u ) 1 b c\ Uj ) 1 L c\ u j^xl/a fb) r m j. 



\M m: ab(Pl, ■■■ j Pmi Pai Pb)? Fi m \p u ...,p m ;p a ,p b ) 



S 



(10.1) 



Here the factor 1 / (n s (a)n s (b)n c (a)n c (b)) accounts for the average over the number of initial- 
state polarizations and colours and $(p a • pb) is the flux factor. The flux factor fulfils the 
following scaling property 

HvPa-Pb) = vHPa-Pb) ■ (10.2) 

The function Fj (pi, ...,p m ;p a ,p b ) defines the jet observable and has the same properties 
as the function Fj m \px, ...,p m ;p a ) in Sect. || (more precisely, the factorizability of initial- 
state collinear singularities has to be valid with respect to both p a and pb). All the other 
factors in Eq. ( |10.1[) are analogous to those in Eq. ( |3.1|) . 



10.1 Implementation of the subtraction procedure 

In order to compute the NLO cross section, we can write it as follows 
Vab L °(Pa,Pb;V 2 F ) = / (da^(p aj p b ) -da£ b (p a ,pS 



+ 



/ da£ b (p a ,p b ) + da^ipaiPb) + do c ab {p a , p b ; fi 2 F ) 

Jm+l Jm Jm 



(10.3) 



where the local counterterm da^ b (p a ,pb) is given by 

dcTab{Pa,Pb) = Nin 



/wW V Z d(f) m+1 (p 1 ,...,p m+1 ;p a +p b + Q)- 

n s {a)n s {b)${p a p b ) S {m+1} 



J2 J2 V ij,k(Pl'--iPm+l>Pa,Pb) *J (Pi, -Pij,Pk, ..,Pm+l,Pa,Pb) 

pairs k^£i,j 
ij 

Z [ V ij(PU -,Pm+i;Pa,Pb) F^iPl, -Pij, ■;Pm+i;PaiPb) + {a <-> 6) 

pairs 

Z! Z [ V k(Pl, -,Pm+i;Pa,Pb) (Pi, -Pk, ~,Pm+i;Pai,Pb) + {a «-> 6) 



i k=£i 



+ [v ai > b ( Pl ,...,p m+1 ;p a ,Pb) Fj m) (p 1 ,...,p m+1 ;p ai ,p 6 ) + (a <- 6)] | . (10.4) 

While the first three terms in the curly bracket exactly correspond to those in the curly 
bracket of Eq. (|8.6|) , the last term is the new dipole contribution introduced in Eq. (|5.135|) . 
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In order to compute the integral of da^ b (p a ,pb), we write it as follows 

d(Ta b {Pa ,Pb) = da^b (Pa , Pb) + (Pa , Pb) + d(T%" (p a , p b ) + da^b" (Pa,Pb) , (10.5) 

where the four terms on the right-hand side are in one-to-one correspondence with those 
in the curly bracket of Eq. (|10.4j) . 

The integration over p, in the first three terms on the right-hand side of Eq. ( |10.5| ) is 
completely analogous to that carried out in the Sect. || Thus we obtain: 

Mn — TTT^ C E d(pm(pi, --jPmlPa+Pb + Q) 



d(Tab(Pa,Pb) 



m+l 



n s (a)n s (b)$(p a p b ) ^ } 
-J—F ( j m \p l , ...,p m ,p a ,Pb)YY\M^(p 1 ,...,p m - 1 p a ,p h )f 



a s 1 



47T/i 2 



2vrr(l -e) \2pi-pkJ T\ 



Viie 



2 1 \ J ' 



10.6) 



/ da%'(p a ,p b ) = - [ N in f dx 

Jm+1 Jm JO 



1 



E d4> m {pu ■■,p m ;xp a + Pb + Q) 



S{m} 

a s 1 



o n s (a)n s (b)$(xp a p b ) ^ } 

Fj(Pu -,Pm\ XPa,Pb) E \ M m,ab(PU ~,Pm\ %Pa, Pb)\' 

i 

2 \< 1 



Vi(x;e) + (a^b) 



m+l 



27rr(l-e) \2pi-p a J Tj 
d&ab"(Pa,Pb) = - I MlnJ2 [ dx d(f> m (pi,..,p m ;xp a +p b + Q) 

Jm r i JO 

{m} 

■ ~ET— Fj m \pi,--,Pm]XPa,Pb) 
0{m} 



(10.7) 



EE 

k c 



1 



\M%!' db (jpi,...,p m -,xp a ,p b )\' 



n s (c)n s (b)$(xp a p b ) 



a s 1 



47T/i 2 



2.r(l- £ )V2 W .W 5? V-(«; £ ) + , (10.8) 

where the functions Vj(e), Vi(x; e), and V a,b (x; e) are respectively defined in Eqs. ( |7.21| , |7.22| ), 
Eqs. (PngjOD and Eqs. ^ETj - ^Ml - 

Let us now consider the ^-integration of da^ b m {p a ,p b ). We first use the phase space 
convolution in Eqs. (|5.149|,p!l50 ) in order to factorize the Pi integration. Then we can 
integrate the splitting function V M > C over p^ and we find: 

/ u in f dx e E 1 

Jm JO r , 1 i 



daT(Pa,Pb) 



m+l 



{m+i} i n s (ai)n s (b)$(xp a p b ) 
1 



d(j) m (pi, ..,p m ] Xp a + Pb + Q) "FT 

*J{m+l} 

Fj(pi, ~,Pm] Xp a ,p b ) \M a ^ h ~. h {px, ...,Pm\ Xp a ,p b )\ 2 



a s 1 



47T/i 2 



27rr(l-e) \2 Pa -pbJ T 
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V a ' ai (x;e) + (a «-> b) , (10.9) 



where the functions V a ' at (x;e) are given in Eq. ( |5.155|) . In order to rewrite Eq. (|10.9| ) 
in terms of a sum over m-parton configurations, we have to perform the corresponding 
counting of symmetry factors. However this counting is exactly the same as that already 
considered in Sect. || (see Eqs. ( |8.17| , |8.18| )) for the case in which the spectator is a final-state 
parton k. Thus we obtain: 

/ da^"'(p a ,p b ) = - [ MnYj I dx d(j)m(Pl,",Pm;XPa+Pb + Q) ~J— 
Jm+1 Jm rj. JO J{m} 

■ F<j m \px,..,p m ;xp a ,p b ) 

' — 77^1 V \M c ^ cb {p u ...,p m ;xpa,Pb)\ 2 

V n s (c)n s (b)$(xp a p b ) 



Collecting Eqs. ( gDg , gg77[ , gUg[ , [mD[ ) and adding Eq. QQ5D , we can write the following 
expression 

1 



/ da£ b (p a ,p b )+ I da^ b (p a ,p b ;fi 2 F ) = ^2 I dx f dy f Mi 

Jm+1 Jm cd m 



'' n s (c)n s (d)§(xyp a p b ) 



d(j) m (pi, ..,p m ; xp a + ypb + Q) F { J n \p l , ..,p m ;xp aj yp b ) 

{m} 

mtCd < l,..,m;a,b\r b ' cd (x,y;e)\l,..,m;a,b> mtCd . 



10.11) 



Here J m+1 da^ip^pi,) + J m da^ b (p a ,p b ; ix 2 F ) is obtained from the leading-order contribution 
d a a\{ x Pai yPb) by replacing the corresponding Born-level matrix element squared 



n s (a)n s (b) ' 



ab <C. — — ^m.ab 



by 



V 1 

7^ n s (c)n s (d) 



d <....\I ab > cd (x,y;e)\....> miCd 



(10.12) 



(10.13) 



and performing the x and y integrations. The insertion operator I(x, y; e) depends on the 
colour charges, momenta and flavours of the QCD partons. Its explicit expression can be 
written as follows 

I ab ' cd (p u ...,p m ;p a ,x;p b ,y;e; n F ) = S ac 5 bd 5(1 - x) 5(1 - y) I(p h ...,p m ,p a ,p b ; e) 
+ 5 ac 5 bd \5(l-y) (pi, ...,p m ;p a ,x;e) + 5(1 - x) I {1) (p u ...,p m ;p b ,y;e) 



+ 



5 bd 5(1 - y)!^^, ...,p m ;p a ,x;p b ;e; fi F ) + 5 ac 5(1 - x) I$?(pi, -,Pm;p a ;Pb, V, e; /4) 



■(2) 



b,d, 
(2)> 



(10.14) 



where 



I(p 1 ,...,p m ,p a ,p b ;e) 



as 1 
27rr(l -e) 



E^H(6) 



Zpi'Pk, 
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+^2 V.(e) 



£7VT a [^— ) +T b -T a ' 



+ H(e) 



2pi -Pa y 



47T/1 



2 \ e 



+ T b -T a 



\2Pb 'Pa, 



;io.i5) 



/ { i)(pi,...,p m ;p a ,a;;e) 



V^;e)-5(l-x) H(e) 



1 



10.16) 



^2) (Pi > • •• > Pm', Pa, x; p b ; e; fi%) 



as 



1 



(2)v^ ..., jr ,„, jrtl ,~, jru , ^, r ^ 27rr(l - e) 
V a ' c (x; e) - 5 ac 5(1 - x) V«(e)l + T b • T c ( 



Ann 2 

Zpi'Pa, 



V a ' c (x; e) - 5 ac 5(1 - x) V Q (e) 



+ 



/ 47T/X 2 

4f/i 2, 

/4 



P ac (x) +K™{x) 



10.17) 



All the contributions that are not proportional to 5(1 — x) in Eqs. ( 1016f) and ( 10.17 ) 
respectively come from Eqs. ( |10.7| ) and ( |6.13| , |10.8| , |10.10| ). The operator I in Eq. ( |10.15| ) 
instead contains all the terms coming from Eq. (|10.6|) plus those proportional to 5(1 — x) 
that have been subtracted in Eqs. ( |10.16| ) and ( |10.17| ). 

We see that I(p±, ...,p m ,p a ,p b ; e) in Eq. ( |10.15| ) is exactly like that in Eq. (|7.26|) , apart 
from depending on the additional initial-state partons a and b. Therefore it cancels all the 
singularities in the virtual contribution J m da^ b (p a ,p b ). 

The other two operators J(i) and 7(2) contribute as finite counterterms. As a matter of 
fact, the insertion operator in Eq. (|10.16|) is exactly the same as the insertion operator 
in Eqs. ( |8.26| , |8.29| ) for the case of cross sections with a single incoming parton (note, 
however, that in Eq. ( |10.16j ) colour-charge conservation reads J2i T{ = —{T a + T b ) ). As 
for the operator 1(2), in order to show that it is finite itself for e — > 0, we rewrite Eq. ( [10. 17 ) 
as follows 



■ffin(Pi> ■■■,Pm;p a ,x;p b ;e;/j, 2 F ) 



(2) 

+ T b -T c 

. i 

+ T b -T c 



2nT(l 



47T/i 2 

^Pi ■ p a 



2Pb-Pa) 
f 4lTfi 2 



tyb ■ Pa 



1 




/tt2 

c 

V i 


-J 


rrxl 

C 


1 

c 


V 



V a ' c (x; e) + - P ac {x) - 5 ac 5(1 - x) V a (e) 



=o+T b -T, 



47T/i 2 



2p& • Pa 

V a ' c (x; e)-V a ' c (x;e 



1 

7fa 



'47T/X 2 

, Mf 



- P ac (x) 
e 



10.18) 
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Then, the first term in the curly bracket of Eq. ( 10.18 ) gives: 

2 



K ac (x) - P ac (x) lnx + 5 ac Tl 



2 In- 1 



1 — x 1 — X 



1 — X 



ln(2 - x) 



+ 0(e) ,(10.19) 



where we have used Eq. ( |B.31| ) and charge conservation T% + T b = —T c ). From the 
second term in the curly bracket of Eq. ( |10.18j ) we obtain: 



VT T — In ^ F 



+ Tb -T c \ ln-^— 

C T 2 C 2p b -p a 



P ac (x) + (9(e) 



(10.20) 



where we have used charge conservation and performed the e expansion as in Eq. (|8.36|) . 
Finally, the third term in the curly bracket of Eq. ( |10.18| ) gives: 



T b ■ T c -zK ac (x) + S ac T b ■ T a 



2 1 \ 2 
In + ln(2 - x) 

1 — x 1 — xj + 1 — X 



+ 0(e) 



(10.21) 



where we have used Eq. ( |5.155| ). Collecting these results we find: 

I^( Pl , ...,p m ;p a , x- Pb - e; = {-K a \x) - P ac (x) lnx + K°, 

+ 5 ac T 2 n 



x 



2 , 1 \ 2 , , 

In- +- ln2-x 

1 — X 1 — X J + 1 — X 



T z c 2 Pi ■ Pa 



+ T b -T c ^K ac (x)+5 ac T b -T a 



T c 2 Pb ■ Pa 
2 In 1 



P ac (x) 



. 1 — X 1 — X J + 1 — X 
and, adding Jm, we can write the following final expression 

d(Ta b (p,P)+ dag{p,p;fi 2 F ) = \da* b (p,p) ■ I(e 

Jm Jm L 



ln(2 - x) 



m+l 



(10.22) 
0(e) , 

(10.23) 



Jm 



+ J2 dx K a ' a, (x)-da^P,P)\+J2 dx / P a ' a '(xp,x^ 2 F )-d^, b (xp,p) 

^7 JO Jm 1 J ^ JO Jm 1 

+ J2[ 1 dx I \K b ' b '(x)-dag l (p,x P )}+Y, t dx I \P b ' V (xp,x;^ F )-da^(p,xp) 



Jm 



where /(e) is given by Eq. ( |10.15|) and the insertion operators K a ' (x) and P a ' (xp, x; up) 
are: 



(*) = g \K°* (:r )-/v ;"<(,-) 



(10.24) 



1% 



1 — X 



+ 6{l-x) 



—T b ■ T a ,^-K aa> (x) 
2tt T\, K J 



P a ' a '( Pl , ...,p m ,p b ;xp a ,x;n F ) 



27T 1 ; T a , 



^T,-T a , In ^ +T b -T a , In 

i 2xp a ■ Pi 



2xp a ■ p b 



. (10.25) 
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The operators /(e) and P a,a (xp, x; li f ) are completely analogous to those in Eqs. (|8.25| ) 
and ( |3.39| ) for the case with a single incoming parton, apart from the trivial dependence 
on the additional initial-state parton. Note, instead, a new feature of the operator K. 
While the term in the curly bracket on the right-hand side of Eq. ( |10.24| ) is equal to that 
in Eq. ( |8.38| ), in the present case there is an additional contribution to K, namely K aa (x), 
due to parton-parton correlations in the initial state. 



10.2 Final formulae 

The results of the previous Subsection can be combined into the following final expressions 
for the jet cross sections in hadron-hadron scattering processes. 

The hadron-level cross section is obtained by convoluting the partonic cross sections 
on the right-hand side of Eq. ( 6.10| ) with the parton densities. The LO parton-level cross 



section is given by 

^ b °(Pa,p b )= I da*(p a ,p b ) (10.26) 

J m 

{ 1 

= / d$ {m) (p a ,p b ) — — \M m ,abipl, ...,Pm]Pa,Pb)\ 2 F ( f l \p l , ...,p m ;p a ,Pb) , 

J n c (a)n c (b) 

where a and b denote the flavours of the incoming partons, p a and p b are their momenta 
and n c (a),n c (b) are their number of colours. The matrix element |-M mja 6| 2 is the square 
of the tree-level amplitude to produce m final-state partons and Fj 11 ^ is the function that 
defines the jet observable we want to compute (the properties that Fj has to fulfil are given 
in Eqs. (|7.2| -|7.4|)) and Eqs. ( p.3| , |3.4| )). All the other phase space factors on the right-hand 



side of Eq. ( |10.1|) are collected into the factor d^ m \p a ,p b ). 

According to the general notation in Eq. (|2.15| ), the full NLO partonic cross section is 
obtained by adding three different types of contribution: 

NLO I 2\ NLO{m+l), n NLO{m}, \ / lnn7 \ 

°ab (Pa,Pb;H F ) = °ab (Pa, Pb) + O ab S {Pa,Pb) (10.27) 

f 1 A \ ~NLO{m}, 2\ , ~NLO{m}, 2\~ 

+ / dx a ab J {x;xp a ,p b ,fx F ) + a ab s {x; p a , xp bl bt F ) . 
Jo 1 J 

The first contribution has m + 1-parton kinematics and is given by the following ex- 
pression 



6=0 



<?a b LO{m+1} (Pa,Pb) = / (d0-ab(Pa,Pb)) _ 

- ( J2 da^ b {p a ,p b ) ® ( dV Aipolc + dV^ ipole ) 

dipoles 

= f d$ {m+1 \p a ,p b ) | - 1 \M m+1 ,ab(pi, ■■■,Pm+l,Pa,Pb)\ 2 F! J m+1 \p 1 , ...,p m+1 ;p a ,p b ) 

J { n c {a)n c (b) 

£ (P-F( m ))( Pl ,...,p m+1 ;p a ,p 6 )l , (10.28) 

dipoles J 
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where Ai m +i >a b is the tree- level matrix element with m + 1 partons in the final state and 
Sdipoies {p ' F^ m ^ (pi, ...,p m+ i,p a ,p b ) is the sum of the dipole factors contained into the 
curly bracket on the right-hand side of Eq. ( [LP .4 ). 

The NLO contribution with m-parton kinematics is obtained by adding the virtual 
cross section and the first term on the right-hand side of Eq. (|10.23| ). As in all the other 
scattering processes, its explicit expression is given in terms of the square of the one-loop 
matrix element \M. m ,ab\ 2 (i-i oop ) and of the insertion operator /(e): 



NLO {m} 



ab 



(Pa,Pb) 



+ 



,ab< 1, 



d(7ab(Pa,Pb) + da^ b (p a ,p b ) ® / 

1 



e=0 
2 



n c (a)n c (b) 
m; a, b | /(e) |1. 



\M mi ab(pi, ... j Pm] Pai Pb) | (l-loop) 



(10.29) 



,m:a,b > 



m,ab 



e=0 



1- 



,Pm',Pa,Pb) 



In the present case, the colour-charge operator /(e) is explicitly written down in Eq. ( |10.15| ) 
(see also Appendix C). 

The third term on the right-hand side of Eq. ( |10.27| ) comes from the second and third 
lines on the right-hand side of Eq. ( |10.23|) and contains two contributions that are similar 
to that involved in the processes with a single incoming hadron (cfr. Eq. (|B.4(j|)). Each 



of these contributions is obtained by integrating a cross section with m-parton kinematics 
with respect to the fraction x of the longitudinal momentum carried by one of the incoming 
partons. When this parton is the parton a, we explicitly have: 



dx Va b L ° {m} (x;xp a ,p b ,fi 2 F ) = dx \ \da^, b (xp ai p b ) 

, JO Jm 1 

T / dx d^ m \xp a7 p b )F^\ Pu ... 

i Pmi XPai Pb) 

ha , b < l,....,m;x Pa , Pb \ (K a ' a ' (x) + P a ' a ' {xp a ,x; n 2 F f 



[K + P) a ' a (x) 



e=0 



(10.30) 

, ,fTl,Xp a ,p b > mj a'6 ; 

where the colour-charge operators K and P are respectively defined in Eqs. ( |10.24| ) and 



( |10.25| ) (see also Appendix C). The expression for cr 



NLO {m} 
ab { 



x; p a , xp b , ftp) is completely 



analogous to Eq. (|10.30|) , apart from the replacements xp a — > p a ,Pb %Pb and J2 a ' ~> 12b' 
(as in Eq. ( 10.23| )). Note that the right-hand side of Eq. ( 10.3U| ) has exactly the same 
structure as in Eq. ( [B.49 ) for the case with a single incoming parton. However, we should 
recall that the colour-charge operator K a ' a entering into Eq. ( |10.3Qj ) differs from that 
appearing into Eq. flS.49j) by the additional correlation term K a ' a (see Eq. ( |10.24| )), which 
is due to the presence of the other incoming parton b. 

Equations (|10.26|) , ( |10.28|) and ( |10.30|) are directly evaluated in four space-time dimen- 
sions. As for Eq. ( p.0.29| ), one should first cancel analytically the e poles of the one-loop 
matrix element with those of the insertion operator /, perform the limit e — > and then 
carry out the phase-space integration in four space-time dimensions. 
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11 Multi-particle correlations 



In the case of processes involving multi-particle correlations (see Eq. ( |6.20| )), the partonic 
cross section at the Born level is given by 

d<Ta b , iinc i )ai ,.., an (p,p-, qi, .-, q n ) = A/ '- ns ( a ) ns ( 6 ) nc ( a ) nc ( fo ) $ ( p .p) 

• d 4>m(pi, -,Pm',p+P+Q -qi - ■■■ - q n ) -J— 

{m} ^{ m i 

• \M m +a 1+ a 2 ...,ab{qii -AmPl, • • ■ , Pm] P, P) \ 2 (fel], -[q n ], Pi, ■■■,Pm]P,P) ■ (H.l) 

Here, we denote by a and 6 the flavour indices of the two incoming partons with momenta 
p and p, while a±, a n are the flavour indices of the final-state identified partons with mo- 
menta qi, q n . In addition, the leading-order cross section in Eq. (|1 1 . 1|) has m final-state 
unidentified partons with momenta pi, ...,p m (non-QCD partons are understood). The jet 
defining function Fj has the properties already discussed in Sects. |7HT0], namely, infrared 
and collinear safety and factorizability of initial- and final-state collinear singularities. Re- 
member that, by definition, the momenta p, p, qi, q n are supposed not to be parallel to 
each other. 



11.1 Implementation of the subtraction procedure 



According to our general procedure, we write the NLO partonic cross section fl6.23|) in the 
following form 

NLO / 2 2 2\ 

(T ab,{incl)a l ,..,a n \PiP, <?1> <?n, ^Fi P'li ■■i\ i n) 

{d&ab,{incl) ai ,..,aSP->P'i —J ?n) ~ ^ 7 a6,(inci)ai,..,a„ (P> P\ 9l> In)) 



+ 



m+1 



m+1 



^ CT a&,(wcl)o 1 ,..,o„(P)P) 9l> •••> 9n) + / ^ cr a6,(mc/)oi,..,a„(P) P> 9l> —j 9n) 



m 

+ / da abAincl)a 1 ,..,a n (P->P'-' 9l> "•' 9"! ^F' ^4) 
Jm 

where the subtraction term is defined by 



(11.2) 



da 



ab,(incl)ai,..,a n {P> Pi 9l) •••> (7n) -A/in . . —\ 

A ; 11 ' n s (a)n s (o) ( I ) (p • pj 

^ # w+ i(pi, ...,p m+ i;p + p + Q - gi - ... - 9n) ^— — 

{m+1} »J{m+l} 

E (*> • ^i m) L + EE(^^L + E \ip ■ F J m) r + ( V - F 



pairs 
i.j 



1=1 i 



p (m) 



bi 



(11.3) 



Here we have introduced the shorthand notation T> ■ F 1 ^ to denote the different dipole 
contributions in which the emitter is a final-state unidentified parton, a final-state identi- 
fied parton or an initial-state parton. Their explicit expressions, according to the dipole 
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formulae in Sect. |5|, are respectively the following 
(p.if 1 >) = £ V ij , k {q u ..,q n ,p 1 ,...,p m+1 ;p,p) F^\[q 1 },..,[q n },p 1 ,..p^p k ,..,p m+1 ;p,p) 



+ J2 T> i3Al(<lh ■■> < ln,Pl, -,Pm+i;P,P) (M,--, [ft] , "fan] , Pi, -Pij, -, Pm+l] P, P) 



1=1 



+£&(?i> ••, q n ,pi, -,Pm+i;p,p) Fj m \[qi], ••, [q n },pi, -Pij, -,Pm+i;p,p) 



+ T^iMl, ■; Qn,Pl, -,Pm+l,P,P) F J ([Qi], ■; fan], Pi, -Pij, ■; Pm+1, P, p) 



;n.4) 



(V-F^j m) ) . =J2' D a l i,k{ ( ll>-; ( ln,Pl,--,P m +i;P,p) 

• *J ([?l]> ••> [ft], Pi) - ; Pfe, ■;Pm+i; P ,P) 

n 



r-=l 

■)( n ) b / 



+ V a,i (qi,-,Qn,Pl,-,Pm+l,P, ],Pl,-,Pm+\,P,P) , (11-5) 



P • Fj m) )" + (a <- 6) = {E^f(9i> Pi, ...,Pm + i; P ,P) 



■ Fj^iiQl],-, [Qn},Pl,-,Pk,-,Pm+i;P,P) 



+ 



^V^ a \q 1: ..,q n , Pl , ...,p m+1 ; P ,p) +V^ ai ' b (q h ..,q n , Pl , ..., Pm+1 ;p,p) 



1=1 



f j "* ([?i]> ••> [ftJ,Pi, -,Pm+i;P,p)} + (a <->• 6) 



11.6) 



Comparing Eqs. ( 11.3HTTT6 ) with the form of the subtraction terms introduced in Sects. |7]- 
|I0| , we see that the only new feature of the present case is due to the 'pseudodipoles' 
W^l arJ V^j a ,V^ ) i \v^ ai } V ( -^ ai ' b where both the emitter and the spectator are incoming 
partons or identified final-state partons. These dipoles are defined in Sect. |5l| Note, in 



particular, that no distinction is made between initial- and final-state spectator (i.e. T> i 



(n) a 
ai% 



V%l a in Eq. ([TOD and £>M ai > fe = Vf' ai in Eq. fllOD ) and that the dipole partonic states 



{gi, .., q h ..q n ,p l: ..,p m+ x,p,p}, {q u .., q n ,p u .., Pm+ i,p,p}, {q 1: .., q n ,p u ..,p m+1 ,p,p} depend 
only on the emitter (the corresponding jet functions Fj appear as common factors in 
Eqs. (mjXq) ). 



The subtracted contribution (da R — da A ) in Eq. (|11.2|) is integrable in four dimensions 
by explicit construction. 



In order to evaluate the <i-dimensional integral of da A , we decompose it as follows 

da A = I i (da Al + da A " + da Al " + da A "" 



m+l 



11.7) 



The first three terms on the right-hand side of Eq. ( |1 1 . 7|) respectively correspond to the 
integral of (V-Fj)ij, V aii ^ and T>™. Their treatment has been already considered in Sects. |7|- 
H The fourth term contains the integral of all the pseudodipoles and can be handled as 
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da A "" in Sect. [HJ As a result of this integration procedure we find: 



m+l 



,(incl)ai ,..,a n (p,p;qi,...,q n ) + I doS, (incl)ai,..,a n 



/ daft, 

(inci)ai ,..,a n 

(p,p; gi,...,g n ) -1(e) 

Jm L 



+ E / dx / [(^ a * a '(x) + P a ' a '(xp,x;/i^)) •d^6,(incO« 1 ,.,«n( a ^'P; 

q/ JO Jm L N ' 

+ E / [(K b ' b '(x) + P b > b '{xp,x;n 2 F )) ■da^ !{ind)ai ^ an (j),xp;q 1 ,...,q n ) 

+ EE / — / d ^,(inci)a 1 ,..,a'..a n (P»P;9l I " > 9i/2," > ?n) ' + P a > ai (qi/ Z, Z', 



11.8) 



The factor /(e) comes from the integration of X^fc and, in addition, collects all the e- 
poles due to the other dipole factors. The finite parts of the dipoles P?-, Pf -, D^ 4 , T>^ ai , 
jy(n)ai,b^ j)(n)bi,a con tribute to the initial-state operators K(x) + P(xp, x; fx F ), and those of 
the dipoles V i j an V aii ^,V^} ar ,V^} a J V^ b contribute to the final-state operators H(z) + 
P(qi/z,z)fj.f). 

The insertion operator /(e) in Eq. ( [11. 8 ) is fully symmetric with respect to all the QCD 
partons. Therefore, denoting by I a generic QCD parton (/ = {i,ai,a,b}), we can write 
(the singular factors Vj(e) are given in Eq. ( [7.27] ) ) : 

I(q 1 ,...,q n ,p 1 ,...,p m ,p,p;e) = - ^ — - — - ^ V/(e) E T/ • Tj ( ^ ) . 

27rr(l-e) Y T / J^/ VPi-PjJ 

(11.9) 

The initial-state operator P a,a (xp, x; is instead symmetric with respect to all the 
partons except p. Its explicit expression is: 

P a > a \q 1 ,..,q n ,p u ...,p m ,P;xp,x;fi F ) = ^P aa \x) 7 ^J2 T i- T «> ^7T^~ • ( U - 10 ) 



rrtz 

1 a' Ijia' 



2xp ■ pi 



Similarly, for the final-state operator P a ' >ai (qi/z, z; [if) we find 



Pa'ai(qi,~,qn,Pl,-,Pm,P,P;qi/z,z;i4) = - 1 P a ' ai (Zj ., 



1 ^ _ „ , ZM? 



E T ' ■ T < In 



2gz • p/ 

(11.11) 



The insertion operators K a ' a (x) and H a > m {z) are separately symmetric with respect 
to the sets of the unidentified and identified (or initial-state) partons. They are given by 
the following equations 



_ f — aa 



K a ' a fx) 



K aa (x) - KZ (x) 



1 — x 



+ 5(1 -x) 



^^pT ar T a , + T b .T a ,jK aa '( 



J2T ai -T a , £ a > a '(x;p,q h n) + T b -T a , £ a > a ' {x;p,p, 



n) 



i=i 



11.12) 
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H a' vai {z) 



«S f -j^a',ai 
27T 



K^\z)+W a , ai {z) \nz-HZM 



1-2 



+ 5{l-z)-l 



P a , ai {z) \nz-K^{z) 



1 



T ar ■ T a>l C a >' (z; q h q r , n)+T a - T a[ £ a >> (z; q h p, n) 

r=l 

+T b -T a > C a > l (z-q h p,n) 



11.13) 



Tab 



where the flavour kernels K"", K^, K ab , C a ' b , 
Q5|) , (p|), (|5~T56|) , (|5TT76|) and (jsT7p . 



are defined respectively in Eqs. 



K32 



Note that the operators I and P in Eqs. ( |11.9HITTTl"| ) are completely analogous to 
the corresponding operators defined in Sects. f7HTo|. Some new features instead appear 
( |Tl~T2D and 



in 



Eqs. 



11.13|) . Comparing Eq. (|1 1 . 1 2| ) with Eq. ( |10.24j ), we find a new 



contribution (that in the last square bracket of (|11.12|) ) due to correlations between the 
incoming parton a and the identified final-state partons. As a result, the insertion operator 
K explicitly depends on the momenta of the identified partons unlike the previous cases. 
This contribution is indeed non- vanishing only if there are final-state identified partons. A 
similar term is also present in the expression (|11.13|) for H . However, in this case also the 
term in the second square bracket on the right-hand side of Eq. ( 11.13 ) has no analogue 
in the corresponding Eq. ( |9.19| ). This term is due to final-state parton correlations and 
is similar to that proportional to K aa (x) on the right-hand side of Eq. ( |11.12j ). In this 
respect the difference between initial- and final-state parton correlations simply amounts 
to the replacement K a '' a (z) — > K a '' a (z) — P a '' a (z) hiz. 



11.2 Final formulae 



The results of the previous Subsection can be summarized by the following final formulae. 



The LO parton-level cross section, which enters into the calculation of the multi-particle 
hadronic cross section of Eq. ( |6.21| ), is obtained as follows 



a ab°(inci)a 1) ..,a n {P->P'-' 9i> •••> In) — / da^ b ^ ncl - )ai ^ an (p,p;qi,...,q n ) 



d<& {m \p, 



1 



n c (a)n c (b) 



(qi, ...,q n ,Pi, 



F j m) ([?l], -[<ln\,Pl,-,Pm,P,P) ■ 



,p m ;p,p)\ 2 

(11.14) 



Here n c (a) and n c {b) are the number of colours of the incoming partons, M. m +a 1 +a 2 ...,ab is 
the tree-level matrix element to produce m final-state partons in addition to the identi- 
fied partons, and Fj m ^ is the most general jet defining function (it fulfils Eqs. ( |7.2| - |7.4|) , 
Eqs. (|8.3| , |3.4|) and Eqs. fl9.2|J9.3| ) with respect to the dependence on the momenta of the 
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unidentified partons, of the incoming partons and of the final-state identified partons). 
The term d& m \p,p; qi, q n ) in Eq. ( 11.14|) collects all the other kinematic factors on the 
right-hand side of Eq. (|11.1|) . 

As in the symbolic notation of Eq. ( 2.15Q , the NLO partonic cross section contains a 
contribution with m + 1-parton kinematics, a contribution with m-parton kinematics and 
n + 2 terms obtained by a one-dimensional convolution of cross sections with m-parton 
kinematics. We have: 

NLO ( 2 2 2\ _ NLO{m+l} , x 

a ab,(incl)a 1 ,..,a n \PiP^ 9l, •••) Qn] H'Fj Pi) ••> PrJ ~ a ab,(incl)a 1 ,..,a n \PiP^ HU — , Qn) 

f 1 , r ~NLO{m} i _ 2\ , ~NLO{m} , - 2\" 

J u 

+ Z^/ n a ab,(incl)L..,aS Z '>PiP'i^--^l/ Z i-^n;Vl) • (H-15) 



The NLO contribution with m+ 1-parton kinematics is given by the following expression 

(p,p;qi,...,q n )) _ n 



jVLO{m+l} / - \ 
tT o6,(» n d)a 1 ,.:,« n (P.P; 9l.- 5 9n) 



m+1 



e=0 



(P, Pi 91) •••) <7n) ® ( <^Klipole + ^Kiipole 

dipoles 

/ r TTx |AC+l+ai+a 2 ...,a6(gi, -,9n, Pi, -, Pm+1, P, P) | 2 

n c (a)n c (o) 



ci$( m+1 )(p,P;gi,...,g„ 



•^'(fclr-W.Pl.-.Pm+lSftP) - X! ( P ' F(m) ) (9l>->9n,Pl,-,Pm+i;P,P) f , 

dipoles J 

(11.16) 

where .M m +i+ 0l + a2 ... )0 & is the tree-level matrix element with m + 1 unidentified partons 
in the final state and Sdipoios \P ' F^ m ^ (gi, .., q n ,Pi, — ,Pm+i,P, £>) is the sum of the dipole 
factors in the curly bracket on the right-hand side of Eq. (|11.3| ) . 

The NLO term with m-parton kinematics is obtained by adding the virtual cross section 
and the first contribution on the right-hand side of Eq. (|11.8| ). Its explicit form is given by: 

NLO{m} / - x 

a ab,{incl)a 1 ,..,a n \P->Pi 9l> — > Qn) 



,(incl)ai,..,a n (p,p; ?l,.»,?n) ®! 

1 



e=0 



d$M(p >p; qi ,...,q n ) 



\M m +a (Qi, ...,q n ,Pi, —,Pm',P, P)l (1— loop) 



n c (a)n c (b) 

m+a\+a2...,ab < ^ <^1, ••, <^n, 1; , Wlj <^> & | -^"(^) |^1, ••; 1) •••5 ^, ^ ^m+ai+a2...,ab 

• *J -[?n]. Pi,-, Pm,P,P) , 



e=0 
11.17) 



where |A / i m+ai+a2 ... ia b|^ 1 _ ioop x is the one-loop matrix element squared and the colour-charge 
operator /(e) is defined in Eq. ( |11.9| ) (see also Appendix C). 
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The last three terms on the right-hand side of Eq. ( 11.1 5|) are in one-to-one correspon- 
dence with the last three terms on the right-hand side of Eq. ( |11.8| ). Their expressions as 
a function of the QCD matrix elements are the following 



, * NLO {m} i - 2\ 



= J2 dx da a'b,(incl)a 1 ,..,a n ( X ^ X P^P'^^--^n)®(K + P) a ' a ( y a 

a , JO Jm L 

= E / dx I d ^ {m) ( x P,PiQU-,Qn) -F/ m) (M, -[ffn], Pi, -,Pm;xp,p) 
- Jo J 

■ m+a 1 ..,a'b< Ql, ; Qn, h 1™', ^P, P I (-K"°' a ' (x) 

+ P a,a '(xp,£;/4)) \qi,;qnA,-,m; x P,P>m+ai..,a'b , 



e=0 



(11.18) 



1 dz 



a 



NLO {in} i I 2\ 

ab,(incl)L-.,aS Z 'i V,V\ Ql, Ql/Z, q n \ ft ) 



z" 



E 



1 dz 



2r Jm 



e=0 



• m+ .. +a ;+.. )ab < gi, .,®/z, ..,g„, 1, ..,m;p,p| (if^O) 

+ Pa'^im/Z,*)^)) \qi,;Ql/z,..,qnA,-,m;p,p> m +.. + a' l +..,ab ■ 



The contribution ^^inciWi a„( x ' P' X P' 9i> 5™ A*! 1 ) * s obtained from Eq. (|11.18|) by the 



>,(inc()a_ 

replacements xp p, p — > xp, X) tt ' ~^ Sfe' ( see Eq. flU-D)- Note that each of the 
contributions ^ab^ncVflx a n m Eqs. ( |1 1 • 1 5| , [TT . 1 8| , |1 1 . 1 9|) depends on a single factorization 
scale. 



The colour-charge operators K and P of Eq. ( |11.18| ) are defined in Eqs. (|11.12|) and 
C P- 1.1 01") respectively. The colour-charge operators H and P of Eq. ( T1.19Q are respectively 
given in Eqs. ( |11.13|) and ( ll.llj) . The definitions of the related flavour kernels P aa (x), 

K (x), K™ (x), K aa (x), H^ a {x) and of the functions £ a ' a are also recalled in Appendix 
C. 



The actual evaluation of Eqs. (|11.14j) , ( |11.16|) , (jll,18[) and ( |11.19|) is directly performed 
in four space-time dimensions. As for Eq. ( |11.17 ), one should first cancel analytically the 
e poles of the one-loop matrix element with those of the insertion operator J, perform the 
limit e — * and then carry out the phase-space integration in four space-time dimensions. 

Note that the formulae presented in this Subsection can be applied also to the simplified 
cases of multi-particle correlations with a single incoming hadron or with no hadrons in the 
initial state. For this purpose, it is sufficient to remove one or both of the contributions in 
the square bracket on the right-hand side of Eq. ( |11.15| ) and to set equal to zero the colour 
charge of the removed incoming parton both in the dipole factors of Eq. (|11.16|) and in the 
colour-charge operators of Eqs. ( |1 1 . 1 7| — P~l . 1 y| ) . 
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The results obtained in this Section show that, in the most general case of multi-particle 
correlations, the NLO partonic cross section in Eq. ( |1 1 . 1 5 ) is finite, i.e. free from soft and 
collinear singularities. We should point out that its finiteness directly follows from the 
primitive definition in Eq. ( |6.23| ), which is nothing other than the formal restatement, to the 
lowest non-trivial order, of the factorization theorem of mass singularities 0]. Doubts have 
been raised in the past about the validity of the theorem for processes involving more than 
one identified parton, for instance, the Drell-Yan process. Although those doubts have been 
proved to be unfounded, a Cartesian (and fully general) proof of the factorization theorem 
in the context of QCD is still missing. Our check, by means of the explicit calculations in 
this Section, that the partonic cross section in Eq. ( |6.23D is finite is, to our knowledge, the 
first proof (at NLO) of the factorization theorem in the most general case of QCD jet cross 
sections. 



8N 



12 Summary and discussion 



12.1 Summary 



In this paper, we have introduced a new general algorithm for calculating arbitrary jet cross 
sections to NLO in arbitrary scattering processes, based on the subtraction method. The 
key ingredients are the dipole factorization formulae, which implement both the usual soft 
and collinear factorization formulae, smoothly interpolating the two. The corresponding 
dipole phase space obeys exact factorization, so that the dipole contributions to the cross 
section can be exactly integrated analytically over the whole of phase space. 

The steps that are necessary to set up a general method for evaluating NLO QCD 
cross sections are recalled in Sect. [l|. In order to compute a jet quantity to NLO, one 
must calculate the contributions from one-loop corrections to the Born-level cross section 
and from the cross section for processes in which one additional parton is present in the 
final state. Each leads to singularities, which are regularized by working in a number 
of dimensions d = 4 — 2e other than four. Analytical integration must then be used to 
extract the singular terms as poles in e and combine them with one another to yield a 
finite result. However, almost all experimentally important jet quantities are sufficiently 
complicated that analytical integration is impossible and one must resort to numerical 
techniques. Thus one must somehow extract the singular parts of the cross section in a 
way that is independent of the exact details of the observable and treat them analytically. 
This leaves a remainder that depends on the full complications of the jet quantity, but 
which is finite so can be treated numerically. 



One way of doing this is provided by the subtraction method, described in Sect. 
This works by introducing an approximate cross section, which is defined in such a way 
as to match all the singularities of the real cross section to produce one additional parton. 
Thus the difference between the two is guaranteed to be finite and numerical integration 
can be used in 4 dimensions. However, the approximate cross section is required to be 
simple enough to be integrated analytically without knowing the details of the definition 
of the jet observable. 

This is achieved by considering, for every m + 1-parton configuration, a mapping to 
a 'similar' m-parton configuration (or in general, several of them), defined such that in 
the singular regions of phase space the two configurations become indistinguishable. The 
approximate cross section is then defined to be proportional to the jet quantity calculated 
from this similar m-parton configuration. Thus, the analytical integration can be performed 
without knowing the definition of the jet quantity, by holding the m-parton configuration 
fixed and integrating over all m + 1-parton configurations that map onto it. This gives 
rise to e poles, which cancel those from the one-loop cross section to produce that same 
m-parton configuration. The result is a finite m-parton cross section that can then be 
integrated numerically. 

Note that, using the subtraction method, no approximation is actually performed in the 
evaluation of the NLO cross section. The approximate cross section is subtracted from the 
real cross section and then added back to the one-loop cross section. Moreover, and most 



89 



importantly, no approximation is performed in the analytical integration of the approx- 
imate cross section. In that respect, the adjective 'approximate' may sound misleading. 
Rather than approximating the true cross section, the subtracted contribution defines a 
fake cross section that has the same dynamical singularities as the real cross section and 
whose kinematics are sufficiently simple to allow its analytic integration. 

To define this fake cross section, one must find an approximation (in d dimensions) 
to the matrix element (or its square) that matches the real matrix element in all the 
singular regions of phase space. While this can be done by calculating the full matrix 
element in d dimensions and taking all the relevant limits, this is an extremely laborious 
and ungeneralizable procedure. Instead, one can use an approach based on the soft and 
collinear factorization theorems (see Sect. |j), which guarantee that the singular terms are 
process independent. Thus one can define the subtraction cross section as a sum of soft and 
collinear pieces in a completely process-independent way and integrate it once-and-for-all. 

In fact, we go one step further even than this. By introducing a dipole factorization 
theorem (see Sect. |5[), we are able to construct a single approximation to the matrix element 
that matches all of its soft and collinear singularities. This approximation is given in terms 
of the Born-level matrix element times universal (process independent) dipole factors. Thus 
we can provide a single formula that approximates the real matrix element (squared) for 
an arbitrary process, in all of its singular limits. This avoids problems that can arise 
from using the separate soft and collinear approximations, where one must either define 
an explicit arbitrary cutoff to separate the two regions, or add back on another term to 
compensate for the double- counting of the region of overlapping singularities that gives rise 
to double poles. 

While the soft and collinear factorization theorems dictate the form the dipole factor 
must take in the exactly singular limits, we are free to choose the extrapolation away from 
those limits arbitrarily. Clearly we should use this freedom to make the necessary analytical 
integrals as straightforward as possible. Physical processes in which some external momenta 
are fixed (by incoming hadrons or by measurement of outgoing hadrons) impose additional 
constraints on the phase space. In order to ensure that the phase-space integrals are still 
analytically tractable, we make different choices for the extrapolation away from these 
limits, depending on which of the participating parton momenta are fixed by external 
hadrons. Thus we end up with a finite set of different dipole formulae (see Sects. |Q|-|5.6|), 
applicable to different physical processes. 

Closely related to the choice of dipole formulae is the definition of the 'similar' m-parton 
configuration. Each dipole factor smoothly describes the merging of three partons into two 
new partons (emitter and spectator) while one of the three partons approaches a singular 
region. Thus, for each of the different constraints on the m + 1-parton phase space, we are 
able to define a one-to-one mapping to a set of m momenta plus a single-parton momentum 
in such a way that the single-parton subspace obeys exact phase space factorization. That 
is, the m + 1-parton phase space can be written as the product of a physical m-parton 
phase space times the dipole phase space. By physical, we mean that all partons are on- 
shell and energy-momentum conservation is implemented exactly, as are all the phase-space 
constraints. The fact that each dipole factor depends on more than two parton momenta 
is essential in order to implement these kinematic features. 
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Owing to the exact phase-space factorization and the convenient choice of dipole for- 
mulae, it is possible to integrate all of the dipole contributions analytically over the full 
dipole phase space in d dimensions. These result in a set of e poles that cancel those in the 
one-loop cross section, as well as a set of finite corrections (see Sects. |772| , |7.3| ). It is worth 



noting that because of the simple definitions of the dipole formulae, the origin of all the 
finite terms that arise can be simply traced, and all are well-known constants, allowing a 
powerful check that the integrals have been performed correctly. 

In addition to the poles that cancel with the one-loop contribution, when there are 
identified external hadrons the integration of the dipole factors leads to additional poles that 
must be subtracted into the process-independent parton distribution functions. Practically, 
this means cancellation against universal (but factorization scheme- and scale-dependent) 
collinear counter-terms (see Sect. |]). The scheme- and scale-dependences resurface in the 
finite remainder left over after the cancellation (see Sect. [2.3| and Sects. |8HTT|). Once again, 
the finite terms can be easily checked, because they have simple physical origins related to 
different integrals and projections of the Altarelli-Parisi splitting functions. 

One feature of our algorithm is that it does not require the convolution with the parton 
distribution function to be made during Monte Carlo integration. One is free to choose 
either to calculate a hadron-level cross section, including the convolution, or a parton- 
level cross section as a function of the partonic momentum fraction. The latter can then 
be convoluted with the distribution function after Monte Carlo integration. This can be 
extremely useful in many respects. One can produce cross sections with a wide variety 
of parton distribution functions, or study the scheme- and scale-dependence of the re- 
sults without having to reintegrate for each new scheme or scale. Moreover one can check 
whether, in extreme phase space regions, the NLO partonic cross section contains en- 
hanced (typically, logarithmically-enhanced) contributions that may spoil the convergence 
of the fixed-order perturbative expansion, thus requiring all-order summations. Finally, the 
parton-level calculation can also be important for comparing different computations, since 
all the dependence on non-perturbative input can be removed, and the result is a purely 
perturbative, parameter- free, quantity. 

Starting from the dipole formulae and the integrals of the dipole factors, in Sects. 0-1 1 
we have explicitly carried out all the d-dimensional analytic work that is necessary for a 
straightforward numerical implementation of NLO calculations in any scattering process. 
The results are collected in effective final formulae, which are recalled in the last Subsection 
of Sects. f7HTl| for each different class of scattering processes. Using these final formulae, any 



NLO calculation requires only the corresponding matrix elements as input. More precisely, 
the only additional ingredients needed to construct a numerical program to calculate the 
NLO corrections to arbitrary jet quantities in a given process are: 

• a set of independent colour projections of the matrix element squared at the Born 
level, summed over parton polarizations, in d dimensions (if the total number of 
QCD partons involved in the LO matrix element is less than or equal to three this is 
unnecessary, because the colour structure exactly factorizes, see Appendix A); 

• the one-loop matrix element in d dimensions; 

• an additional projection of the Born level matrix element over the helicity of each 
external gluon in four dimensions; 
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• the tree-level NLO matrix element in four dimensions. 



We should emphasize that, independently of the actual set-up of our algorithm, the 
dipole formalism is fully general and, hence, highly flexible. The main point is that it 
provides an explicit, universal and simple (in d dimensions) expression for the approximate 
cross section da A . Starting from it and having the process-independent dipole splitting 
functions and their integrals to hand, by direct inspection one can try to modify the sub- 
traction term (for instance, introducing finite weighting factors or cut-offs for the dipole 
terms) to simplify its treatment in any particular scattering process. This may be useful 
for improving the convergence of the numerical program and can always be done at the 
expense of some extra analytic work in four dimensions on the finite difference between the 
two approximate cross sections. In other words, using the dipole formalism one can set up 
ones own sub-algorithm. 

Generalizing the procedure for constructing NLO Monte Carlo programs for arbitrary 
quantities has several advantages. These are principally because of the reduction in the 
number and complexity of ingredients that have to be calculated for each new process, 
and because the d- dimensional integrals only need be done once and can be easily checked 
independently, rather than being buried inside a specific calculation. 

Perhaps the single biggest advantage is the fact that the NLO matrix element for addi- 
tional real emission can be calculated in four dimensions. In calculations involving several 
partons this can result in great savings in computation time and in size of the final ex- 
pressions, because helicity amplitude techniques can be used. On the other hand, specific 
calculations that construct the approximate cross section directly from the real one must 
work in d dimensions, producing very cumbersome formulae at intermediate stages even 
though the final result is simple, since it is just the soft and collinear limits, which factorize. 

It is often said that the bottleneck in producing new NLO calculations is the calculation 
of the necessary one-loop amplitudes. While this is partially true, one can cite many 
examples where the relevant matrix elements have been available for a long time, yet no 
Monte Carlo programs for arbitrary jet quantities have been available. Most notorious has 
been the case of jet production in deep inelastic scattering, where the matrix elements can 
be simply obtained by crossing the e + e~ ones, which have been known for many years, 
but only recently have any NLO jet calculations been produced pi , pl| . This has severely 



hampered jet studies by the HERA experiments, which have been forced to compare data 
with parton shower models, or partial calculations. Other examples include the longitudinal 



fragmentation function in e + e annihilation, which was, until very recently ||32|| , perhaps 
the single simplest uncalculated QCD quantity. 

It is thus clear that the numerical implementation of NLO calculations forms a second 
bottleneck. General algorithms such as ours will certainly help to reduce the amount of 
work needed to make these implementations, so will help to reduce this bottleneck and 
hence increase the number of processes in which perturbative QCD can be compared with 
data in a quantitative way. 

As a final comment on the Monte Carlo implementation of the subtraction method, 
we should note that the cancellation of soft and collinear singularities does not completely 
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solve all the numerical problems. Although the subtracted cross section is finite, integrable 
square-root singularities may arise in the m+1 parton integral. Square-root singularities 
cannot be integrated by naive Monte Carlo methods because they have infinite variance 
and the integration procedure would never converge. The presence of these singularities 
can depend on the behaviour of the particular jet observable, so no subtraction procedure 
can universally overcome them. However, it is straightforward to control the square-root 
singularities in a completely general way, without any knowledge of the jet quantity being 



calculated, using the standard technique of multi-channel Monte Carlo integration |33 



during the phase-space generation. In the Monte Carlo implementations of our algorithm a 
user-supplied routine analyses the generated momenta, and any infrared safe jet observable 



can be used |2Q,21 



12.2 Comparison with other general approaches 

The first general algorithm for constructing NLO Monte Carlo programs for jet cross 
sections in arbitrary processes was proposed in Refs. MM, using the phase-space slicing 
method. Several subsequent calculations have been based on it. In the simple case in 
which there are no identified hadrons, discussed in Ref. |J, colour-ordered amplitudes are 
used to derive the soft and collinear factorization formulae. An arbitrary, but small, pa- 
rameter is introduced to separate the phase space into 'resolved' and 'unresolved' regions. 
In the latter the cross section and jet quantity are really approximated by their soft and 
collinear limits. The resulting formulae can then be integrated analytically over the unre- 
solved regions, to cancel the poles in the virtual cross section. 

Although this is an approximate calculation, it becomes exact in the limit that the 
cutoff parameter goes to zero. Unfortunately the errors on the numerical integral over the 
resolved region of phase space then diverge and one must always make a compromise in the 
choice of cutoff. While rules of thumb have developed to give a rough idea of how small it 
needs to be, there is no substitute for explicitly checking that there is no dependence on it. 
However, this is not always feasible owing to constraints on computer time. A particularly 
poignant example is that of the energy-energy correlation in e + e~ annihilation. Since 
discrepancies exist between dedicated analytical calculations and general-purpose Monte 
Carlo programs, a very high-statistics comparison between three different Monte Carlo 
programs was made ||34|| , to check that they really were in agreement. The result was a 
small, but statistically significant, discrepancy between the program of Ref. || and those 
of Refs. 0,^3). While this was attributed to a residual dependence on the cutoff, it was 
not possible within the available computer time to make a sufficiently accurate calculation 
to confirm this, or to extrapolate to zero cutoff. 

These problems become increasingly severe as one approaches the edges of phase space 
because the approximation performed within the slicing procedure can strongly inter- 
fere with the actual definition of certain jet observables. Indeed, at the edges of phase 
space there are physical parameters, namely ratios of physical scales, that become large. 
These large parameters can produce logarithmic and even power-like (!) enhancement of 
the cutoff-dependence, thus requiring such small cutoffs that numerical stability is never 
reached, in practice. 
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However, we should stress that these problems do not imply that, in some case, the 
phase-space slicing method is unable to provide reliable QCD predictions. On one hand, in 
the evaluation of jet quantities far from the phase-space boundary, these problems are all 
obviated if proper care is taken to explicitly check the cutoff-dependence for every result 
quoted and to work in a region in which it is negligible. On the other hand, in the case 
of quantities at the edges of phase space, the kinematic region where one can lose control 
of the numerical stability is that in which the perturbative expansion is affected by very 
large coefficients: in this region the convergence of the fixed-order expansion is spoiled, the 
NLO calculation is insufficient to provide reliable QCD predictions and the latter require, 
anyhow, analytic summation techniques to all orders. 

A possible disadvantage of the phase-space slicing method regards precisely this last 
point. For most jet observables the only analytical NLO results available are for the co- 
efficients of logarithmically enhanced terms near the kinematic limits. Thus, it is a useful 
test of a complete NLO numerical calculation to compare it with these analytical results. 
At the same time, the numerical calculation can be used to check these analytical results, 
which are the first step in the resummation procedures. The approximation embodied in 
the phase-space slicing method, and the ensuing problems of numerical stability in ex- 
treme phase-space regions, may reduce the amount of information that, otherwise, can be 
provided by a NLO calculation. 

The slicing method of Ref. || was extended to include identified hadrons in Ref. |J. 
The approach is to first consider the case in which all partons are outgoing and then cross 
some to the initial state as necessary. As we have seen in our approach, the kinematic 
constraints imposed by identifying external hadrons do not act symmetrically on initial- 
and final-state partons. Therefore 'crossing kernels' have to be introduced. 

These crossing kernels are similar to the insertion operators K and P (and H , in the 
case of fragmentation processes) that, in our approach, arise as finite remainders left over 
after the subtraction of the universal collinear counter-terms. The main difference between 
the crossing kernels of Ref. || and our kernels K and P is that the latter do depend on 
the colour charges of the partons involved in the scattering process. In order to compute 
the full parton- level cross sections in NLO, the crossing kernels as well as the insertion 
operators K and P must be convoluted with m-parton cross sections (see Sect. |2.3| ). 

In Ref. @ , the crossing kernels are preconvoluted with the parton density set to provide 
a new effective set of 'crossing functions', which are used in the main integration stage of 
the Monte Carlo program. This procedure can be a considerable problem in one of the 
main applications of NLO jet calculations, the extraction of parton distribution functions 
by fitting to the data, since this process is usually iterative, with the distribution functions 
gradually converging on the best fit. The use of the crossing functions would require a 
new Monte Carlo integration for each iteration. In our method on the other hand, it is 
straightforward to calculate parton-level cross sections. These can then be convoluted with 
any parton distribution functions after the Monte Carlo integration is completed. Thus 
each iteration would then only require a simple one- dimensional numerical integral. 

This disadvantage of the usual implementation of the method of Ref. can obviously be 
avoided if the crossing kernels are first convoluted with partonic cross section contributions 
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rather than with parton densities. 



The general properties of soft and collinear radiation were used to construct a subtrac- 
tion algorithm for the first time in Ref. 0, for one- and two-jet production in hadron-hadron 



collisions. This has recently been modified to deal with three-jet production [16, T7] and, 
in general, with n-jet cross sections. Although this formalism is based on the subtraction 
method, so looks superficially similar to our algorithm, there are in fact a great many 
differences. 



The subtraction procedures used in Refs. [lj| and [|17| differ in many details but share 



some common features. Firstly, they select energy and angle variables by working in a 
definite reference frame and thus breaking Lorentz invariance at intermediate steps of the 
calculation (although of course it is restored in the final results). The definite choice of 
a reference frame then unambiguously specify the integration variables that can lead to 
singularities in the integration of the real matrix element. Thus, one can introduce a 
partition of the m + 1-parton phase space in such a way that in each subregion only one 
energy or angle variable can kinematically vanish. Having done that, the approximate cross 
section da A is defined by means of a double (soft and collinear) subtraction procedure 
as the product of the singular variable times its residue. The residue is evaluated by 
exactly performing the soft and collinear limits and thus it can be computed in a process 
independent manner. 

The two main differences between this 'residue approach' and our formalism regard 
the treatment of the soft and collinear limits and the related definition of the subtracted 
m-parton configuration. While the dipole formulae provide a single and smooth approx- 
imation of the real matrix element in all of its singular limits, in the residue approach 
the soft and collinear regions are treated separately. Correspondingly, while in the dipole 
formalism the subtracted parton configuration is obtained by a one-to-one mapping from 
the original m + 1-parton configuration, in the residue approach the mapping is achieved 
by a projection procedure. This projection of the m + 1-parton phase space onto m-parton 
phase space does not allow for an exact phase space factorization. For this reason, one 
has to introduce arbitrary cutoffs to define the upper limits on the integration variables, 
although the dependence on these cutoffs should cancel, to within the numerical accuracy, 
in the final result. 

The residue approach can be recovered as a particular set-up of the dipole formalism. 
Starting from our general expression for the subtracted cross section da A one can project 
it onto any properly defined soft and collinear subspace. 



12.3 Future outlook 



Looking to the future, there are several avenues along which the present work could be 
continued. The first is obviously to apply it to as many processes as possible, in particular 
those for which no other calculations exist. We have already constructed Monte Carlo 
programs for jets in e + e~ annihilation p0[ and deep inelastic scattering [fJTJ and several 
more applications are in progress. 
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In terms of extending the algorithm itself, the most obvious missing feature at present is 
the treatment of heavy quark effects. It is straightforward to extend the dipole formalism to 
incorporate massive partons, either in jet calculations or with fragmentation functions [|l9|j . 



This is clearly an important extension, as many heavy quark processes are good probes of 
perturbative QCD, but many of these that are being experimentally measured are not yet 
predicted to NLO. 

Another potentially important extension is the generalization to polarized partons, 
which is also straightforward in the dipole formalism. 

Looking further ahead, it is to be hoped that at some stage NNLO calculations of 
jet observables will be attempted. Even once the necessary two-loop matrix elements are 
calculated, the amount of work needed to provide a numerical implementation will be enor- 
mous. Clearly any progress that can be made in the meantime to set up a general-purpose 
NNLO subtraction algorithm will speed up the process of bringing the new calculation to 
the marketplace. The dipole formalism seems particularly suited to this task. 

When starting this research project, we had in mind a main final goal: a method for 
carrying out NNLO QCD calculations. Having set up a completely general NLO algorithm, 
we are confident that this main goal can indeed be achieved. 
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Appendix A: Colour Algebra 



In order to practice with the colour algebra we can consider some simple examples. In the 
simplest cases with two or three partons (regardless of whether they are final- or initial-state 
partons), the colour algebra can be performed in factorized (closed) form. 



For the case with two partons, using colour conservation, we have: 

T 1 -T 2 1 1, 2 >= -T 1 -T 1 1 1, 2 >= —T\ 1 1, 2 >= -T\ |1,2 > , 



(A.l) 



so that all the charge operators {T 2 , T 2 ,, —T\ -T 2 } are factorizable in terms of the (scalar) 
Casimir operator C± — C 2 . 

Using colour conservation for the three-parton case we have: 

= (E T *) |1,2,3> 

= (T 2 1 +T 2 2 + Tl + 2T 1 -T 2 + 2T 1 -T 3 + 2T 2 -T 3 ) |1,2,3> , 



and 



(T 1 -T 2 + T 1 -T 3 ) |1,2,3>=-T? |1,2,3> . 
Combining these two equations we get: 

2T 2 -T 3 1 1, 2, 3 >= (T\ -T\ -if) |1,2,3> 



(A.2) 
(A.3) 

(A.4) 



and similarly for T 1 ■ T 3 and T\ • T 2 . Therefore, all the charge operators are factorizable 
in terms of linear combinations of the Casimir invariants Ci, C 2 , C 3 of the three partons. 

When the total number n of partons is n > 4 the colour algebra does not factorize any 
longer. For instance, if n = 4 we have four trivial relations, namely 



7t|l,2,3,4>= Ci|l,2,3,4> , 



1,...,4 



(A.5) 



As for the remaining six charge operators Tj • Tj(i ^ j), we can use the following four 
identities (due to charge conservation) 



T,-^T J |1,2,3,4>=0 , 



1,...,4 



(A.6) 



in order to single out two independent charge operators. For instance we can write: 



T 3 -T 4 |l,2,3,4> 

7VT 4 |1,2,3,4> 

7VT 3 |1,2,3,4> 
7VT 4 |1,2,3,4> 



l - (d + C 2 - C 3 - C 4 ) + T 1 ■ T 2 
l - (Ci + C 3 - C 2 - C 4 ) + T 1 ■ T 3 



- (C 4 — C\ — C 2 — C 3 ) — T\ ■ T 2 — T\ ■ T 3 
-(C 1 + T 1 -T 2 + T 1 -T 3 )|l,2,3,4> , 



1,2,3,4> , 
1,2,3,4> , 

1,2,3,4> 



(A.7) 
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and express all the charge operators in terms of Casimir invariants and Ti-T 2} Ti-T 3 . The 
actual values of the latter depend on the detailed colour configuration of the four-parton 
state, namelyQ 



< 1,2,3,4|2VT 2 |1,2,3,4> 

< l,2,3,4|Ti.r 3 |l,2,3,4> 



M 



b\b2 Q3 124 



(Pl,P2,P3,P4) 



T b c iai T b c 2a2 Mr 2a3a4 (PnP2,P3,P4) 



^^ 2bsa4 (pi,p 2 ,P3,P4)]*T b % l T b C 3a3 ^^ a2 ^(p 1 ,p 2 ,p3,P4) 

(AA 



In the general case with n > 4 partons, one should consider n(n + l)/2 charge operators 
Ti-Tj. Among them there are n trivial factorizable contributions, i.e. T?\... > n = Cj|... > n . 
In addition, one can use n charge conservation constraints (i.e. • 2? =1 Tj|... >„= 0) to 
end up with n(n — 3)/2 independent charge operators whose action onto the matrix elements 
has to be evaluated explicitly. 



** To b e precise, if some of the partons in |1, 2,3, 4 > are initial-st ate pa rtons, the right-hand sides of 
Eq. (A.S) should contain an additional normalization factor as in Eq. ( |3.13| ). 
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Appendix B: Soft Integrals 



In this Appendix we explicitly perform the integration of the soft terms Vi )Cl b in Eqs. Q5.174I ) 
and (|5.191|) over the phase space volumes in Eqs. ( |5.171| ) and (|5.188|) respectively. 

Let us first discuss the space-like case (see Eq. ( |5.173| )). We have to carry out the 
following integration 

1 



Isoft(x) = 8na s p. e / [dpi(n,p a ,x)] 



Vi,ab 



(b.i; 



Since the vector n M in Eq. ( 5.157 ) is time-like, we can work in its rest frame and consider 
the integration variables Ei,8,(p defined as follows 



rr 



Pi 



e(i. 



E a (l,. 



Pa+Pb = E ab (l, ■ ■ . ,vsinx, vcosx) 



Ei(l, ..'angles'.., sin#cos0,cos( 



(B.2) 



In Eq. ( [B.2| ) the dots stand for vanishing components, while the notation 'angles' denotes 
the dependence of Pi on the d — 4 angular variables that can be trivially integrated in 
Eq. O- 

Using Eqs. ( |5.161| , [5~171|) for the phase space [dpi(n,p a ,x)] and Eq. ( |5.169|) for and 
introducing the variables in Eq. ( |B.2| ), the integral ( |B.1| ) becomes: 



Isoft(x) = S(x) 9(1 - x) 



a s T(l - e) f°° dEi ( irfi 
2^ r(l - 2e) 



E< El 



d0(sin0) 



l-2e 



sin . 



\-2e 



c f E a — Ei 
1 — v cos x 



(B.3) 



/o 7r ' (1 — cos 0) (1 — v sin x sin 6 cos <p — v cos x cos ( 

The integration over the energy E^ in Eq. ( B.3|) is straightforward: 



00 dEi ( vr/i 2 



o Ei \E 



5 x- 



Eg — Ej 

E a 



1 



7T/X 



l-x \{l -x) 2 E 2 at 



(B.4) 



while for the angular integrals we can use the result of Ref . |35| , namely 



d9{sm9f- 2e 
1-v 2 



sin i 



-2c 



1 — v cos x 



7T 



' 1 — cos 6) ( 1 — v sin x sin cos <fi — v cos x cos i 



+ ln 2 



(1 — vcosx) 2 
1-v 1 
2 



2Li 2 



v(l — cosx) 



1-v 



2Li, 



v(l + cosx) 
1 — v cos x 



vcosx 



1 + v 1 
~ 2 



_2 l n 2±^_^ ln 2 



+ 0(e 2 



(B.5) 



— v cos x) 2 

which, by means of the identity Li 2 (— v(l + cosx)/(l — vcosx)) = ~ Li 2 (v (1 + cosx)/(l + 
v)) - (l/2)ln 2 ((l - vcosx) /(I + v)) (i.e. Li 2 (l - 1/x) = -Li 2 (l - x) - (1/2) In 2 x, for 
< x < 1), can be rewritten as follows 

1 — v cos x 



define) 1 - 2 ' 

1-v 2 



sin . 



-2c 



7T 



cos 0) (1 — v sin x sin 9 cos — v cos x cos 6) 



vcosx)' 



+ e 



2Li 2 | 



v(l — cos x) 
1-v 



+ 2Li 2 



' v(l + cosx) 



+ 0(e 2 



(B.6) 
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Thus, we find 



I soft (x) = 0(x) 0(1 - x) ^ ^4 (! - ^'" 1_2 ' 



1 

e 

e 



2Li- 



2tt T(l - 2e) 
/ d(1 — cosx) 



2Lio 



(l-t; 2 )7r/i 2 
-E 2 (l — f cosx) 2 
f (1 + cosx) 



l — v J \ 1 + v 

We can then perform the e-expansion of the factor (1 — x)~ l ~ 2<L according to Eq. 

<* ~ I)_I " 2< = -h* 1 ~ x) + (rb) + + 2£ (r 1 ^ ln t^) + + ° (e!) • 

and, relating the variables t>, cosx, i n Eq. (|B.2|) to the Lorentz invariants: 

4£ 2 (1 - v cosx) 2 n 2(1 -v cos x) (p a 

: 2p a ■ p fe , 



) .(B.7) 
(B.8) 



1 



1-^ 



Pa ■ ™ 



n 2 (p a +Pb) 2 



\ [(Pa+Pb) -n] 2 

we end up with the final result (note, T 2 (l - e)/T(l - 2e) = 1 - e 2 7r 2 /6 + 0(e 3 )) 



J so/t (x) = 0(x)e(i-x)^ T ^ T l- T (^) (i^l -x) 



1 / 1 



2 In 1 



2tt T(l - e) \2p a p b J [2e 2 '" e \1 - x 

(l + v) (Pa+Pb) ■ n 



6(1 -x) 



2Li 2 1 



+ 2Li 2 1 



1 — X 1 — X. 

1 (1 —v) (Pa + Pb) -n\ 7T 2 



Pa ■ n 



Pa ■ n 



12 



0(e) 



(B.9) 



(B.10) 



Let us now consider the time-like case (see Eq. ( |5.190|) ) and compute the integral 

1 Vi ab 



Jsoft(z) = 87ra S Ai 6 / [dpi(n;p a , z)\ 



ZPaPi z 



(B.11) 



As in the previous case we work in the rest frame of n M . Using Eqs. ( |5.182|]5T88|) for 
[dpi(n;p a , z)] and Eq. ( p.!69| ) for t>j ab and introducing the kinematic variables in Eq. ( |B.2p , 
we obtain an expression for J so ft(z) that is equal to that in Eq. ( |B.3| ) apart from the 
replacement: 

E a — Ei\ 1 _ 2e r ( E a 



8 x- 



E a 



z L ~' € 8 z- 



E a + Ei 



Therefore the Ei integration in J so ft(z), namely 



00 dEi ( nji 2 



z l ~ 2t 5 ( z 



E n 



7T/i 



/o Ei \Ef J "V E a + Ej l-z\(l-z) 2 E 2 aj 
gives exactly the same factor as in Eq. (p.4| ) and we immediately obtain 

Jsoft(z) = hoft(z) ■ 



(B.12) 



(B.13) 



(B.14) 
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Appendix C: Collection of Main Formulae 



In this Appendix we collect together the main formulae that are needed to implement our 
algorithm for calculating jet cross sections. 

According to our method, the final expressions for the NLO cross sections are given in 
terms of contributions with m + 1-parton and m-parton kinematics, denoted as u NL °{ m + 1 } ; 
a NLO{m} ^ g.NLO{m}^ cross sec tion a NL °{ m + 1 } j s obtained by subtracting the counter- 
term da A from the real contribution da R . The counter-term da A is constructed by using 
the dipole factors introduced in Sect. |. These dipole factors are collected in Table 1, 
where we list the Equations with the dipole definition, the corresponding kinematics and 
the related splitting functions. Remember that, when computing a NL °{ m + 1 } ; a n the dipole 
factors have to be directly evaluated in four space-time dimensions. 

The cross sections a NLO ^ and a NLO ^ are obtained by adding the virtual contribu- 
tion da v and the collinear counter-term da° to the integral of the subtraction counter-term 
da A . In the rest of this Appendix we concentrate on the formulae needed to construct the 
final integral of the subtraction counter-term. 



dipole 
factor 


definition 


dipole 
kinematics 


phase 
space 


dipole splitting 
functions V 


Vf 
V 

u b i 

ai,bi ai 


o 

$.61|) 

$.94) 

$.114) 
$.136) 
(5.159) 
(5.180) 


(U, o, (H) 

(P7|), (^38p, O 
C|5.62-5.64j), (|539D 
(15.95-5.97), 

(|5.115|), (15.116), (5.1211) 
(15.137-5.140), 

(|5.157|), (p. 160-5. 164)), (|5.169) 

(|5.157|), (|5.181), (|5.182|), 
(5.162H5.164|), (|5.169|) 


$.20) 

(psp 

(5.72) 

$.103) 

$.124) 

$.151) 

(5.172) 

$.189) 


(ET7-53) 

(5.39-5.41) 

(|535|-5.68) 

(5.98-5.100) 

(5.117-5.120) 

(5.145-5.148) 

(5.165-5.168) 

(5.183-5.186) 



Table 1: List of the dipole factorization formulae used to construct the universal subtraction 
term da A . The numbering refers to the Equations in Sect. 
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Flavour kernels 



The majority of the functions and constants we encounter are derived from the spin- 
averaged Altarelli-Parisi splitting functions in d = 4 — 2e dimensions: 



< P qq (x; e ) > = < P qq (x; e) > = C F 

< P qg {x; e) > = < Pg g {x; e) > = C F 

< P gq (x;e) > = < P gq (x;e) > = T R 

< P gg {x; e) > = 2C A 



l + x 2 



1 — x 

l + (l-:r) 2 
x 

2x(l - x) 



e(l-x) 



— ex 



1-e 

x 1 — x 



l-x 



+ 



X 



+ x(l — x) 



(C.l) 

(C.2) 

(C.3) 
(C.4) 



and < P q q(x; e) > = < Pq q (x; e) > = 0. It is sometimes useful to express these as a power 
series in e, 

< P ab {x- e)> = < P ab {x- e = 0) > - e P^x) + 0{t 2 ) . (C.5) 

The regularized Altarelli-Parisi probabilities are the sum of these and the virtual correc- 
tions, evaluated in four dimensions, 

l + (l-x) 2 



P« J (i) = P q9 {x) = C F 

PM(X) = P 9 %X) = T R 

pii(x) = pm(x) = C F 



X 



x 2 + {l-xf 
(l + x 2 \ 



1 — X 



+ 



P"(x) = 2C A 



l-x 



H l + x(l-x) 

X 



+ S(l-x)(^C A -^N f T R 



(C.6) 
(C.7) 

(C.8) 
(C.9) 



Note that the final-state Altarelli-Parisi probabilities P a b{ x ) are identical to these initial- 
state ones, P a b( x ) = P ab (x), so we do not usually make any distinction. Their regular parts 
are denoted as follows 



P±(x) = P a \x) 



if a 7^ b , 



P£(x) = ~Cf (1 + x) , P r %{x) = 2 C 



l-x 



X 



-l + x(l-x) 



(CIO) 



After integrating the dipole formulae for final state emission, we obtain the constants 

7, = 7* = |c> , lg = l ±C A - 2 -T R N f , (C.11) 

and 



K q = K- q 



IT 



C F , K g ^(^-^)c A -^T R N f , (C.12) 



which are related to various integrals of the Altarelli-Parisi splitting functions. As a matter 
of fact, we have 



\Y,f Q d* W-z)r <P ab (z;e)>=2Tl±+ la +(K a -^Tl) e + 0(e 2 ) . (C.13) 
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After integrating the expressions for emission from a dipole in which one identified 
hadron participates, we obtain the functions: 



K ab (x)=Pj b (x)+P£(x) In — 



x 



+ 5' 



ah 



■ -) , 2 , 1 - x 
T„ In 



a VI — x x 



-S(l-x) (^ a + Ka -K 2 T 2 a 



(C.14) 



or, explicitly, 



K qq (x) = K m (x) = C 



K q9 (x) = K q9 (x) = P q9 (x)\n- — - + C F x , 

x 

K gq (x) = K gq (x) = P 9q (x)\n^^ + T R 2x(l-x) , 

x 



2 1 1 — x\ 1 — x . 

In -(l+x)ln \-(l-x) 

1 — x x ) + x 



5(1 -x)(5- 7T 2 ) C F , 



K yy (x) = 2C A 



In- -\ +(- — - - 1 + x(l - x) ) In 



6(1 -x) 



1 — X X 
50 



x 



(C.15) 
(C.16) 



(C.17) 
1 — x 



X 



(C.18) 



and = = 0. 

Finally, when two identified hadrons participate, we obtain the following additional 
functions 



K ah (x) = P r fJx) ln(l -x)+ 5 ab T 2 a 



rb mi -*>)-£*<!-*) 



(C.19) 



Although it is not as closely related to the Altarelli-Parisi splitting functions, we in- 
clude in this Section the integral of the pseudodipole splitting function encountered in 
multiparticle correlations, 



£ a > b (x;p a ,p b ,n) = 5 ab 5(l-x)2T 2 a 



ah 



Li 2 1 



(1 +v) (p a +Pb) -n 



+ Li 2 1- 



[1 - v) (Pa+Pb) • n 
2 p a ■ n 



Pa ■ n 



2(p a ■ ny 



1 - 



n 2 (pa+Pbf 



\ [iPa + Pb) 

x dz 



n\ 



f x (17 

Li 2 (x) = - — ln(l - z) 
Jo z 



(C.21) 
(C.22) 



Remember that the four-momentum n M is defined as follows 



E P"a , 
aefinal state 



(C.23) 
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where p? n is the total incoming momentum in the scattering process and the second term on 
the right-hand side is the sum of all the momenta of the identified partons in the final-state. 



In addition, partonic cross sections for processes involving identified hadrons also depend 
on the scheme-dependent flavour functions K£ & (x) and H^(x). In the MS scheme, all are 
zero, 

^ 0. 



K— = H MS 

MS ~~ ba 



(C.24) 



In the DIS scheme, the initial-state functions are given by (see [^6| and Appendix D), 



K 



n 

DIS 



K" i 



X 



x) 



K 



11 

DIS 



x 



1 + x 2 ( 1 

I m 

1 — X 



(x 2 + (1 -x) 2 )lo£ 



- + -(9 + 5x) 



+ 8x(l - x) 



X 



X 



and K 



19 
DIS 



X 



K 



19 
DIS 



X 
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J - fi DIS 



X 



-2N f K°? ls (x),K™ s (x) = K 



ii 

DIS 



X 



(C.25) 
(C.26) 
= 0. 



In the case of fragmentation processes, a factorization scheme conceptually analogous to 



the DIS scheme has been introduced in Ref. 37 



Insertion Operators 

The final result for singular part of the integral of the dipole splitting functions is the same 
in all processes, and depends on the universal insertion operator I. In order to write the 
result in a uniform way for all dipole contributions, we use the notation {p} to denote a set 
of parton momenta, without specifying which are identified and whether they are in the 
initial or final state. / and J are indices over all these momenta, and we obtain 

-£f(ibj? h VM £ T, ' Tj (S)' • (a27) 

Note that the scalar products pi ■ pj are always positive. If either momentum is crossed 
between the initial and final states, our uniform notation ensures that pi ■ pj retains the 
same sign. The universal singular function V/(e) in Eq. (|C27j ) depends only on the parton 
flavour and has the following e-expansion 

V/(e) = Tf - + 7/ \ + 11 + Kj + 0(e) , (C.28) 



where the constants 77- and Kj are given in Eqs. ( |C.11 , C.12 ). The insertion operator / 



enters into the calculation of the cross sections a NLO ^ with m-parton kinematics. 

When there are identified external hadrons, we obtain finite remainders from the in- 
finite subtraction into the parton distribution functions. These finite remainders, which 
contribute to the one-dimensional convolution of the cross sections a NLO ^ with m-parton 
kinematics, contain two terms which respectively depend on the factorization scale and on 
the factorization scheme. 
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The factorization-scale dependent term is proportional to the insertion operator P, 
which is again universal. Using the {p} and I notation again, it is given by 



P a > b ({p}; x Pa , x; 11%) = -± P ab (x) — 2 £ T l ■ T b In 

Z7T * " — 



b I^b 



2xp a ■ pi 



(C.29) 



for initial-state partons and 



pU{p};p.A. *; /40 = £ ^ E ^ • T b in 

Z7r b Ijtb P a ' P 1 



(C.30) 



for final-state partons. Here P ab = P ab are the Altarelli-Parisi probabilities in Eqs. ( |C.6| - 
U.9[) . Note that if we denote by p b the rescaled momentum in the insertion operator 
P (i.e. pb = xp a in Eq. (|C.29|) and pb = p a /z in Eq. ( |C.30|) ), the right-hand sides of 
Eqs. ( U.29| ) and ( |U.30|) only differ by the transposition ab — > 6a of the flavour indices in 
the Altarelli-Parisi probabilities. Nonetheless, as discussed at the end of Sect. |9.1| , we have 
P a,b ({p};Pb,x;/j, 2 F ) ^ P a , b ({p};PbiX;iJ, 2 F ). 



As for the term that depends on the factorization scheme, it is proportional to the 
initial-state insertion operator K and/or to the final-state insertion operator H. When 
there is only one initial-state hadron, the insertion operator K is: 



K a >\x) = ^ K a \x) - K*{x) + 5 ab £ Ti ■ Tb^ 



2tt 



1 — x 



+ S(l-x) 



. (C.31) 



Similarly, when there is only one final-state identified hadron, the insertion operator H is: 



H b , a (z) = ^ ( K ba (z) + 3P ba (z) \nz- H™(z) 



2tt 



+ 5 ab ■ T b ^ 



1-z 



5(1 -z) - 1 



(C.32) 



Likewise, when there are two initial-state hadrons, we obtain 



K a > a x) 



2tt 



K" (x) - K™ (x) 



+ 5™' £2VT a ^ 



1 



1 — X 



+ 5(1 -x) 



2tt T, 



(C.33) 



Finally, in the most general case of multiparton correlations, we obtain 



K a > a '(x) = ^{K^ (x) - KTAx) 



2tt 



+ 5 aa 'Y,T l -T a 



li 



1 — x 



+ 5(l-x) 



1 a' \l=l 



Y,T ar T a , + T b -T a ,)K aa '(x) 



J2T ar T a > £ a > a '(x-p,q h n) + T b -T a , C a ' a '(x;p,p, 



n 



i=i 



(C.34) 
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and 



H a'aX Z ) 



as 
2tt 



+ ^a[ai X/ r ^- ' i ' ^ a 'i ji2 



+ S(l-z)-l 



+ 



Tl 



In 



P a[ai {z) \nz-K a '^{z) 



T 2 , 



J2 T ar • T a>l C a > 1 (z; q h q r , n)+T a - T a[ C a > 1 (z; q h p, n) 

r=l 

+T b -T a[ C a > l (z; qi ,p,n)^ . 



(C.35) 



Note that setting T b = in Eqs. ( C.34 ) and (|C.35 ) we respectively obtain the operators 
K and H for the case of multiparton correlations with a single incoming parton. Likewise, 
setting T a = T b = in Eq. ( |C.35|) we obtain the operator H for multiparton correlations 
in processes with no hadrons in the initial state. 



The definition of the flavour kernels K ab (x), K ab (x), C a ' b , K^(x) and H b ^(x) that 



appear in Eqs. ( |a3THa3l ) is recalled in Eqs. ( gig ), (KHSD , (|O20| ) and (|a24H026|) . 



In order to evaluate the NLO cross sections with m-parton kinematics, the colour-charge 
operators /, P, K and H have to be inserted into the tree-level matrix elements. This 
leads to the computation of colour-correlated tree-amplitudes. We conclude this Appendix 
by recalling their definition. As above we denote by {p} a generic set of iV parton momenta. 
The square |A / 1 / ' J | 2 of the colour-correlated amplitude has the following expression in terms 
of the coloured tree-level amplitude Ji4 ai " aN ({p}): 



l-M 7 ' J (M)P 



<{p}\T I -T J \{p}> 



(C.36) 



n c (a)n c (b) 



M 



ai..b I ..b J ..a N 



(M) 



T b I a I T bja J -M ai - ai - aj - aN ({p}) . 



Here the labels a and b refer to the initial-state partons in \{p} > and n c (a) and n c (b) is 
their number of colours. The factor 1 / (n c (a)n c (b)) on the right-hand side of Eq. ( |C.36| ) 
comes from the definition of the state vector \{p} > (see Eqs. ( |3.3| , |3.iiD ) used throughout 
this paper. If there is only one initial state parton, or none, then this factor becomes 
l/n c (a) or 1 respectively, as in Eqs. ( [3.13|) and (|3.9|) . 
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Appendix D: Examples 



In this appendix we give some simple applications of our method. 



e + e 2 jets 

We begin with a trivial example: two-jet observables in e + e~ annihilation. We use cus- 
tomary notation for the kinematic variables: Q 2 is the square of the centre-of-mass energy, 
Uij = 2pi ■ Pj/Q 2 and Xj = 2pi ■ Q/Q 2 , where Pi is the momentum of any QCD parton in 
the final state. 

The LO contribution is the parton model process e + e~ — > q(pi) + q(p2), with matrix 
element M. 2 . We average over event orientation, so M. 2 has no dependence on the parton 
momenta. Moreover, we choose the overall normalization of Ai 2 such that the two-parton 
phase space is: 

= dy 12 5(l-y 12 ) , (D.l) 
and the LO cross section in Eq. ( |7.33|) is given by 



\M 2 \ 2 J dy 12 6(1 -y 12 ) Ff\p uP2 ) . (D.2) 



The NLO real-emission process is e + e — > q(pi) + q{p 2 ) + gips), with matrix element 
■Ms(pi,p 2 ,p3). In four dimensions, the matrix element is: 



\M 3 ( Pl ,p 2 , P3 )\ 2 = C/-^ %l+ f 2 \M 2 \ 2 , (D.3) 

Q 2 (1-X 1 )(l-X2) 

and the phase-space is given by 

O 2 

rf$ (3) = -^d Xl dx 2 0(1 - x 1 )e(l - x 2 )Q(x 1 +x 2 - 1). (D.4) 

l07T Z 

The calculation of the subtracted cross section involves the evaluation of two dipole 
contributions: ©13,2 and ^23,1 • Their definition is given in Eqs. ( |5.2| , |5.7|) . The associated 
colour algebra is trivial, as shown in Eq. (|A.1|) , and we find 

V 13>2 (px,p 2 ,p 3 ) = —-V^ \M 2 \ 2 , (D.5) 

<ZPlP3 

with the following dipole kinematics 

Pf^-Pt, p& = qt--a- (d.6) 

x 2 x 2 

The dipole contribution Z>23,i is obtained from Eqs. ( p.5| , p.6|) by the replacement p± <-> p 2 . 
Inserting the (four-dimensional) definitions of V qig3t2 and Xi, we obtain the final expression 
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for the three-parton cross section in Eq. (|7.35 ): 



a 



NLO {3} 



daf =0 - daf =0 



Z7T JO 



,71 2 I ,-v^2 



(1 - Xi)(l - x 2 ) 



Fj 3 \pi,P2,P3) 



1 — £2 \2 — xi — £2 
1 / 2 



[1 + ari) + 



:i + ^ 2 ; 



x 2 

1 - x 2 



1 — X\ \2 — X\ — X 2 J X\ 

Since the three-parton matrix element can be written as follows 

1 / 2 



(l-x l )(l-x 2 ) 



X 2 \2 — Xi — x 2 



l+^l) + \X\ X 2 ) 



it is straightforward to see that for any infrared safe observable (implying that F 
as Xi — > 1), Eq. ( |D.7|) is finite. 



(3) 

j 



(D.7) 



(D.8) 



,(2) 



Next we have to evaluate the insertion operator /(e) and combine it with the virtual 
cross section. The one-loop matrix element in the MS renormalization scheme is given by 



\M 2 \l 



2 Cfois 

(l-loop) - W vl 2\ 2 



while the insertion operator in Eqs. (|C27| , |(J.28| ) gives 



2 < 1, 21 1(e) |1,2> 2 = \M 2 \ 2 (Vns 



1 /47ru 2 V f 2 3 „ , 



2tt r 

(D.10) 

Note that in Eqs. ( D.S , D.1Q ) the two-parton matrix element J\A 2 is consistently evaluated 
in d = 4 — 2e dimensions. Combining these two contributions as in Eq. (|7.36|) , we obtain a 
finite (for e — > 0) final expression for the two-parton cross section: 



a 



NLO {2} 



J da v + J^da J 



e=0 



71 



J dy l2 5(l-y l2 )Ff\p uP2 ) . (D.ll) 



It is straightforward to check that the total cross section (Fj 3) = Ff> = 1) agrees with 
the well-known result, = \C F ^a£°. 



,(2) 



3 jets 



In Ref. [2(| we presented the simplest non-trivial case: three-jet production in e + e~ anni- 
hilation. For completeness we repeat it here. We again average over event orientation, so 
our formalism can be directly compared with that in Ref. ||. 

The LO partonic process to be considered is e + e~ — ■> q + q + g, with matrix element 
and kinematic variables as defined above for the case of e + e~ — > 2 jets. 
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At NLO, two different real-emission subprocesses contribute: a) e + e~ — > q(pi) +q{j>2) + 
g{p 3 ) + g{pi)\ b) e + e~ — > q(pi) + q{p 2 ) + q{pz) + qiPi)- The calculation of the subtracted 
cross section ( |7.35|) for the subprocess a) involves the evaluation of the following dipole 
contributions: X>i 3i2 , £>i 3 ,4, Pu,2, £>i4,3, ^23,1, £>23,4, XV i> X> 24 ,3, X> 3 4,i, £>3 4i 2- The associated 
colour algebra can again be easily performed because the several colour projections of the 



three-parton matrix element fully factorize (see Eq. (|A.4[) ). Thus we do not need to compute 
any colour-correlated tree amplitudes and we find 



^13,2(Pl,P2,P3,P4) 
^13,4(Pl,P2,P3,P4) 
£>34,lG°l,P2,P3,P4) 



1 - 



V qigz ,2 |-M 3 (Pl3,P2,P4) 



2 Pl p 3 V 2C F 

1 c 

^ l£V«*A I^3(P13,P2,P4) 
2^2^LlV(Pl>^,P34) • 



(D.12) 



The dipole contributions ^23,1, ^23,4, ^34,2 are obtained respectively from ^13,2,^13,4,^34,1 
by the replacement p\ <-> £>2- Analogously, one can obtain X^^ and XV3 respectively from 
Pi 3i 2 and Pi 3i 4 by the replacement p 3 <-> ^4, and XV, 1 and ©24,3 respectively from X>i 3; 2 
and X> 13i4 by the replacement p\ <-> p 2 , £>3 ^ £>4- 

In the case of the subprocess b) we have to consider the following dipole contributions: 
^12,3, £>i2,4, ^14,2, £>i4,3, ^23,1, X> 23i4) X> 34 ,i, "^34,2- Performing the colour algebra we get 



^34,l(Pl,P2,P3,P4; 



1 1 



2p 3 p 4 2 



(D.13) 



and all the other dipoles are obtained by the corresponding permutation of the parton 
momenta. 

The splitting functions of Eqs. ( p,12| , [DTT3|) are explicitly given in Eqs. ( |5.7| -|5.9|). 

The tensor T^ v is the squared amplitude for the LO process e + e~ — > qqg not summed over 
the gluon polarizations (/x and v are the gluon spin indices and —g^T^ = I.M3I 2 ). This 
can be easily calculated. In the case of jet observables averaged over the directions of the 
incoming leptons (un-oriented events) we find (in d = 4 dimensions) 



where 



T^(p h p 2 ,p 3 ) 



ry 2 I ry 2 



\M 3 ( Pl ,p 2 ,p 3 )\ 2 T^ 



M v 
1 — Xi 



2^1-2 1 



t 1 V 

X\ PiPi 



Q 2 

x 2 + x\ 



1 — ^2 



l-x 2 Q 2 
P1P3 + PsPi 



-2- 



x 2 p 2 p 2 



Q 2 



+ (1 + \x\ + \x\ — X\ — x 2 ) g 



Q 2 



x l Q 2 
1 — x 2 — X\ + x\ 

1 — X\ 



P2P3 



(D.14) 



P3P2 



Q 2 Q 2 



(D.15) 



The few ingredients listed in Eqs. (|D.12HD.15| ) have to be combined with the four-parton 
matrix elements A^4 for evaluating the four-parton cross section a NL °^ in Eq. (|7.35|) . 
Obviously, due to the very long expressions for the matrix elements 0], we do not report 
here the explicit formula for a NL °^ . 
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To complete the NLO calculation we also need the virtual cross section. In the case of 
un-oriented events, we take the one-loop matrix element in the MS renormalization scheme 
from Ref. |J (we use slightly different notation): 



\M 3 (pi,p 2 ,P3)\li-ioo P ) = \M 3 (p l ,p 2 ,p^ 1 ' 2 



1 



1 r 



27iT(l-e) V Q 2 
11 



(2C F - C A )yg + C A (y£ + y 23 t )} - - [3C F + -C A - -T R N 



71 



Q's 



+ -{2C F + C A ) - 8C F + [F(y 12 , y 13 ,y 23 ) + 0{e)\ 



2tt 



(D.16) 



where F(yi 2 , yi 3 , 2/23) is defined in Eq. (2.21) of Ref. [fj. 

The explicit evaluation of the insertion operator /(e) in Eqs. (|C27| , |C728D gives: 



3 < 1,2,3| /(e) |1, 2, 3 > 3 = \M 3 (p uP2 ,p 3 )\ 2 



as 



2vrr 1 



47T/! 2 



1 r 



(2C F - C A )yn + C A (yi£ + y 23 e )] + ~ (2 7g + l g ) 



71 



C F 



{2C F - C A ) lny 12 + ~C A \n(y X3 y 23 ) 



Ipg Infers) 



{2C F + C A ) + 2( 7g + K q ) +lg + K g + (9(e) 



(D.17) 



Combining the one- loop matrix element (p,16|) with the result ( p,17|) , and using the explicit 
expressions ( |C. 1 1|JCTT2]) for 77- and Kj, all the pole terms cancel. Note that as well as the 
pole terms, the closely related tt 2 and In terms cancel: 

|A*3(pi,P2,P3)|(i- lo op) + 3< l,2,3|/(e) |1,2,3> 3 = |A^ 3 (Pi,P2,P3)| 2 



as 
2tt 



-l(2C F - C A ) lnyi 2 - ^(5C A - T R N f ) HtoiVza) 



+2C F + ^-C A - ^-T R N f 
9 9 



+ -±[F(y 12 ,y l3 ,y 23 ) + 0(e)) 

Z7T 



(D.18) 



The integration of the expression ( D.18|) (with the matrix element A4 3 given in Eq. ( [D . 3|) ) 
over the phase space in Eq. ( p.4j) provides the three-parton cross section a NL °^ in 
Eq. (pH). 



The results presented here form the basis for a Monte Carlo program |20| that can 
calculate the NLO corrections to arbitrary three-jet-like observables in e + e~ annihilation: 
it will be described in more detail elsewhere. 



1 + 1 jets in DIS and the Structure Function F 2 

We next discuss the simplest case in which there is an incoming hadron: 1+1-jet observables 
in deep inelastic lepton-hadron scattering (DIS). We limit ourselves to considering the 
case of virtual-photon exchange and, in particular, we compute a = ctt + o~l, &t and <tl 
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respectively being the scattering cross sections off transversely and longitudinally polarized 
photons. In the fully inclusive limit (i.e. when no final-state jet observable is defined), 
our cross section is simply proportional to the customary structure function F 2 . We use 
standard notation for the kinematic variables: p is the incoming momentum, q is the 
off-shell photon momentum (q 2 = —Q 2 < 0), x = Q 2 /(2pq) is the Bjorken variable and 
Z% — PiP/PQ, Pi being any parton (hadron) momentum in the final state. The relevant 
matrix elements are obtained from the hadronic tensor W^ u by applying the following 
(d-dimensional) projection operator 

2pq 



pr 



(d-l)Q'< 

(pq) 2 



p^p h 



(D.19) 



Note that to implement the final formulae of our algorithm, the <i-dimensional definition 
of P 2 V is relevant only for a consistent calculation of the one- loop matrix element. 



The hadronic cross section a(p) in Eq. Q6.3Q is obtained by convoluting partonic cross 
sections a a = a q , a q , a g with parton densities f a . Since we are considering only photon 
exchange, by charge conjugation invariance we have a q = a q . Thus, we explicitly compute 
only a q and cr g . At LO, <j^° = while is obtained by the parton model process 
q(p) + 7*(q0 — > q(pi)- By momentum conservation we have p\ = q + p, so the LO matrix 
element A4i >q (q + P',p) has, in practice, no dependence on final-state parton momenta. We 
choose the overall normalization of Aii, q (q + p;p) such that the single-parton phase space 
is: 



5(1 



X 



and the LO cross section in Eq. ( p.45 ) is given by 



ri°(p) = — \M 1 , q (q + p;p)\' 2 5(i 



N, 



x) 



+ p;p) 



(D.20) 



(D.21) 



At NLO, we have to consider real emission processes with two final-state partons with 
momenta p\ and pi- By momentum conservation we have 2pip 2 = (1 — x)Q 2 /x and the 
two-parton phase space is given by 

Q 2 x 



d& 2 \p) 



-dz\ dz 2 



e(zi)6(z2) 8(1 - z x - z 2 ) 2tt 8(tt + fa - fa) . (D.22) 



16tt 2 " l 2vr 2vr 

In the following, we consider jet quantities that are averaged over the azimuthal angles fa t 2, 
so they can be trivially integrated and / d(j)\dfa 5(n+fa—fa) — > 2n. The NLO cross sections 
u q LO (p)i <y^ LO (p) respectively receive contributions from the real emission processes q(p) + 
7*(<?) ~~ > q(Pi) + g(P2) and g(p) +7*(g) — > q(pi) +q(p2)- The corresponding matrix elements 
in d = 4 dimensions are: 



^r\M 2 , q (pi,p 2 ;p)f 



N, 



CU 



87ras 



x 2 + z\ 



[l-x)(l- Zl ] 



2(1 + 3^; 



-jy\Mi, q (q + xp;xp) | 2 



N 2 

c 



■\M 2 , g (pi,P2;p)\' 



T, 



R 



was 



(z 2 + (l- Zl ) 2 )(x 2 + (l-x) 2 ) 



zi(l - Zl ) 



j^\Mi, q (q + xp;xp)\ 2 



+ 8x(l - x) 

(D.23) 
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a 



According to Eq. ( (3.470 , the calculation of the two-parton cross sections a^ L °^ and 
is performed by subtracting, in the first case, the final-state dipole T>\ 2 and the 



NLO{2} 



initial-state dipole Vf 2 and, in the second case, the two initial-state dipoles Df\ D 9 \ ■ The 
associated colour algebra is trivial (as in the example of e + e~^ 2 jets) and we find 



T )g 2 1 (p 1 ,P2;p) 



^-^IAM* + p;p)I 2 , 



(D.24) 



with Vf 2 being obtained from T> 9 2 by the replacement p\ «-> p 2 . Note that the dipole 
kinematics turns out to be the same for all these dipole contributions: the incoming and 
outgoing parton momenta in each dipole respectively are p^ = xp^ and (by momentum 
conservation) + xp 11 . 

Inserting into Eq. (p.24j) the definitions (|5.39| , pl65| , ^.66|) of the splitting functions V and 
combining Eqs. (p.22jJTX2^ , p.24|) as in Eq. ( |S.47| ), we obtain the following expressions for 
the NLO cross sections contributions a^ L °^ 



a 



NLO {2} 



(P) 



f 1 1 

/ dz x dz 2 5(1 -z x - z 2 ) — \M 1<q (q + xp; xp)\'- 
Jo iV c 



Cpas 
2tt 



X 



x 2 + z\ 



_(l-x)(l- Zl ) 

x 2 + z\ 
{l-x){l- Zl ) 



+ 2(1 + 3x^1 
Fj(q + xp; xp) 



Fj 2 \pi,P2;p) 



(D.25) 



a 



NLO {2} 



(P) 



f 1 1 

/ dzi dz 2 5(1 -z x - z 2 ) — \Mx, g (q + xp; xp)\ 
Jo iv c 

(z 2 + (l- Zl ) 2 )(x 2 + (l-x) 2 ) 



Zx{l-Zx) 

x 2 + (1 - x) 2 



ixd — x , 



2 T R a s 

x 

2tt 

Fj 2 \pi,P2;p) 



Fj 1 ^ (q + xp; xp) 



Clearly, for any jet observable (implying that F 
integrals in Eqs. (|D . 25| , pT26[ ) are finite. 



(2) 
J 



F^ as zx 



1, or x 



(D.26) 
1) the 



In order to compute the NLO cross sections with 1 — > 1 parton kinematics, we have to 
evaluate the insertion operators I, P and K in Eqs. ( |C . 2 7| JC72^| , |C . 3 1| ) . We find 

1 



i, q < i;p| 1(e) |1;p >i,. 



\Mx, q (q + p;p)\' 



2tt r(l-e) \ Q 2 ' 
1 



(D.27) 



2 3 
I2- + 7 



+ 10-vr 2 + O( e )| 



£ i,6< 1; zp\P q \zp, z; /4) |1; zp > lj6 = ~ \Mx, q (q + zp; zp)\ 2 ^P qq (z) In ■ 

h c 



Xflp 



(D.28) 
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1 ■ • • ' v, 2 _ as n n QI \ i X ^F 



J2 i,b< 1; *p|P fl 'V, z; /4) |1; sp >i, 6 = |A*i, g (g + zp; zp)\ 2 2 ^P OT (*) In zQ2 



(D.29) 



1 



£ li6 <l;zp|l^> 6 (z) |l;zp> 1>6 = \Mi, q (q + zp; zp)\ 2 ^ (D.30) 

b iV c 27T 



F^)-i^i(*)--CV 



1 



1 - 2 



+ 5(1-^) 



E i,6< 1; *Pl*'' 5 (*) |1; zp > lf6 = ^ |JWi, ff (g + zp)| 2 2 g - 



(D.31) 



According to Eq. (|8.48|) , the result in Eq. (|D.27| ) has to be combined with the following 



one-loop matrix element 
1 



1 



\M 1)q {q + p;p)\ 2 {l _ loap) = j^\M hq (q + p;p)\ 2 



2tt r(i -e) V Q 2 

Obviously, the e poles cancel and we obtain 



C F a s I f^ 2 Yf 2 _3_ g + 0(eU(D32) 



e 



™«(p) = 1 \ Ml , q ( q + p;p)\ 2 %^ (2-vr 2 ) 5(1 - x) (q + p; p) . (D.33) 



Nr. 



2tt 



From Eqs. (p.28| - [D.31|) we can evaluate the NLO cross section contributions in Eq. Q8.49Q . 
Using the phase space factor in Eq. ( p.20| ), we find 

f 1 dz &^ L °^{z; zp, i4) = f 1 dz S(z/x - 1) F^\q + zp; zp) \M 1>q {q + zp; zp)\ 2 
Jo Jo iv c 



■^\K q \z)-K^{z)-\c F 



1 - z 



+ 6{l-z) 



P qq {z) In 



X JJi jp 

zQ 2 



— \Mx, q (q + xp; xp)\ 2 F { p(q + xp; xp) ^ C F x 

iv c Z7T 



(D.34) 



l + x 2 / (l-x)Q 2 3 s 

1 — X \ X fMp 4 1 



1 - 7x 



+ S(l-x) 



39\ 1 



Ti 



8 J C, 



K™{x) 



jf* dz a^ LO{1} (z; zp, /4) = £ dz 5(z/x - 1) F?\q + zp; zp) ^ \M 1>q (q + zp; zp)f 



.g 2 |^(,)-^i(,)-P-(,)ln^} 

= — \Mx >q {q + xp; xp)\ 2 F { p(q + xp; xp) ^ T R 2x 
iv c 271 



[x 2 + (l-x) 2 ]ln (1 X)Q1 



1 



+ 2x{l-x)-—K«{x) 
x n F l R 



(D.35) 
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where, in the last expression on the right-hand side of Eqs. (p.34|]D~3"5| ) we have performed 



the ^-integration and introduced the explicit formulae ( |C.16| , |CTT7|) for the kernels K 
( |C.7|JC.8|) for the splitting probabilities P ab . 



and 



We can explicitly check that our calculation correctly reproduces the known NLO re- 
sults |3(J for the structure function F2. The partonic coefficient functions FJ^ LO for F2 are 
obtained from our cross sections by integrating over all possible final states (i.e. by setting 
1) and by fixing the overall normalization with \Ai\ yq (q + xp;xp)\ 2 = N c 



,(2) 



Gathering together all the terms in Eqs. (|EX25| , pT33ip^) and those in Eqs. (P-26|JD.35|) , 
we obtain 



T?NLO 
h 2q 



In 



In ^ p- - - + - 9 + 5x 

1 - x \ x n 2 F 4/ 4 V 



rpNLO 



(p,fi 2 F ) = ^T R 2x{[x 2 + (l-x) 2 }ln il x)Q * 



2 



+ 8x(l - x) 



C F 
1 



K«(x) 



(D.36) 



m(x) 



rp - l ES 

J-R 



(D.37) 



The NLO expressions in Eqs. ( p.36| , |D~3"TD are factorization-scheme dependent and this 
dependence is accounted for by the kernels and The DIS factorization scheme 

is defined in such a way that F^ q LC \p,Q 2 ) = F 2 N g LO (p,Q 2 ) = 0: thus Eqs. flCTSjlCTD 
follow. The definition of the gluon kernels K^ ls and K^ ls is chosen in order to define 



parton densities that fulfil momentum conservation [36 |. 



The next simplest process involving an incoming hadron is the 2+1-jet rate in DIS. 
Conceptually, there are no additional problems that are not dealt with either in the example 
of 1+1 jets in DIS or in e + e~ — > 3 jets, and it is straightforward to implement a general 
purpose Monte Carlo program for arbitrary 2+1-jet-like quantities in DIS. Specific details 
of the algorithm and numerical results will be presented elsewhere. 
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Note added 



All the detailed calculations in this paper have been performed using the conventional 
dimensional- regularization scheme (see Sect. H). As discussed in Sect. jOj the unphysical 
dependence on the regularization scheme can be parametrized in terms of simple coefficients 
7i (see Eq. ( |3.15|) ) that enter in the one- loop contribution. However, it is worth emphasizing 
some peculiar features^ of one particular regularization scheme, namely, the dimensional 
reduction (or four-dimensional helicity) scheme^]. 



All the final results of our algorithm, summarised in Sects. [7.4| , [872], |9.2| , |10.2| and pL 1.2 



can be directly translated into the dimensional-reduction scheme by simply modifying the 
explicit expression of the cross section component a NLO * m ^ (p) . This is the contribution to 
the NLO cross section that involves m-parton kinematics and no additional convolutions 
with respect to longitudinal momentum fractions. It is obtained by integrating the following 
combination of the one-loop matrix element and the insertion operator I 



77, 



-TT 777 {•M m + ai +a,2~,ab(c[l, •••> Qn,Pl, ■■■,Pm',P, P)\(l-loop) C^-l, 
.yCLjTl c yuj 



+ m+ai+a 2 ...,ab< a l; a ni 1; ••••) 777; (7, b \ 1(e) |di, .., d n , 1, Ui; (7, b >m+a 1 +a,2...,ab \ 



e=0 



The corresponding expression in the dimensional-reduction scheme is obtained by the fol- 
lowing replacements 

\M\l_ loop) -> \M n %_ loop) , (E.2) 
< ...| /(e) |... > - < ...| I DR (e) |... > D . R , (E.3) 

where \-M nK \%-i oop ) and |... > DR respectively denote the one-loop and tree-level matrix ele- 
ments evaluated in the dimensional-reduction regularization and the colour-charge operator 
I DR (e) has the same general expression as in Eq. ( |C.27 ) apart from the replacement: 



Vj(e) - VT (e) = V/(e) - 7/ + 0(e) , (E.4) 

lq = lq = \ C F , % = g C A . (E.5) 

The dimensional reduction scheme is particularly relevant because most of the one-loop 
matrix elements |§ have been computed in precisely this scheme. 

Moreover in dimensional reduction, just because of its definition, the cZ-dimensional 
and four- dimensional tree-level matrix elements exactly coincide. Thus, no calculation 
of ci-dimensional Born-level matrix elements is necessary. Actually, owing to the general 
structure of the one- loop corrections discussed in Sect. |77J| and Ref. p9| , it follows that from 
the e-expansion (see Eq. ( |7.32| )) of |^ nR "|(i_j oop ) one can directly extract all the independent 
colour projections of the matrix element squared at the Born level. 



^We thank Z. Trocsanyi for pointing out these features to us. 

M See: W. Siegel, Phys. Lett. 84B (1979) 193; D.M. Capper, D.R.T. Jones and P. van Nieuwenhuizcn, 
Nucl. Phys. B167 (1980) 479; Z. Bern and D.A. Kosower, Nucl. Phys. B379 (1992) 451; Ref. |1. 
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Thus, for our purposes, the main feature of dimensional reduction is that, provided the 
one-loop matrix elements are computed in this regularization scheme, the shopping list (see 



Sects. [T2] and |12|) to construct a numerical program to calculate the NLO corrections to 



arbitrary jet quantities in a given process can be shortened as follows 

• the one-loop matrix element in the dimensional-reduction scheme in d dimensions; 

• an additional projection of the Born level matrix element over the helicity of each 
external gluon in four dimensions; 

• the tree-level NLO matrix element in four dimensions. 

The computation of one-loop matrix elements in the dimensional-reduction scheme is 
greatly simplified by directly evaluating helicity amplitudes Using these calculations, 
the above shopping list can be further shortened by eliminating the second item. 
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